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ALICE and its pp physics 
programme 



M. Monteno for the ALICE Collaboration 

1.1 Introduction 

ALICE (A Large Ion Collider Experiment) is the dedicated heavy-ion experiment designed 
to measure the properties of the strongly interacting matter created in nucleus-nucleus 
interactions at the LHC energies {^/s = 5.5 TeV per nucleon pair for Pb~Pb collisions) 

mm- 

In addition, with its system of detectors, ALICE will also allow to perform interesting 
measurements during the proton-proton LHC runs at y/s = 14 TeV [3]. Special strength of 
ALICE is the low pr cut-off (~100 MeV/c) due to the low magnetic field, the small amount 
of material in the tracking detectors, and the excellent capabilities in particle identification 
over a large momentum range (up to ~100 GeV/c). The above features of its conceptual 
design as soft-particle (low p^) tracker make ALICE suitable to explore very effectively 
the global properties of minimum-bias proton-proton collisions (such as the distributions 
of charged tracks in multiplicity, pseudorapidity and transverse momentum) in the new 
domain of the LHC energies. 

In addition, these measurements will provide also an indispensable complement to 
those performed in the other pp experiments, ATLAS and CMS, where the superposition 
of minimum-bias collisions at the highest LHC luminosity will be the main source of 
background to the search for rare signals (Higgs boson, SUSY particles, 'new physics^). On 
the other hand also the Underlying Event (i.e. the softer component accompanying a hard 
QCD process) must be carefully understood, since it accounts for a large fraction of the 
event activity in terms of the observed transverse energy or charged particle multiplicity 
and momenta. 

Furthermore, as it will be shown in the following, the ALICE proton-proton pro- 
gramme will include also cross-section measurements of strange particles, baryons, reso- 
nances, heavy-flavoured mesons, heavy quarkonia, photons, and also jet studies. 

Another motivation for studying pp events with ALICE is the necessity to provide a 
reference, in the same detector, to measurements performed with nucleus-nucleus (and 
proton-nucleus) collisions. The latter could be done via interpolation to ^/s = 5.5 TeV 
(the centre-of-mass energy for Pb-Pb runs) between the Tevatron and the maximum LHC 
energy. However, since this interpolation will be affected by rather large uncertainties, 
additional dedicated runs at the same centre-of-mass energy as measured in heavy-ion 
collisions could be necessary to obtain a more reliable reference. Indeed, as it was shown 
by the past experiments at the SPS and at the RHIC, such comparison is important in order 
to disentangle genuine collective phenomena and to be more sensitive to any signatures of 
critical behaviour at the largest energy densities reached in head-on heavy-ion collisions. 

Last, a more technical motivation of pp studies is that pp collisions are optimal for the 
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commissioning of the detector, since most of the cahbration and ahgnment tasks can be 
performed most efficiently with low multiplicity events. 

For all the above reasons, from the Technical Proposal onwards the proton-proton pro- 
gramme has been considered an integral part of the ALICE experiment. This programme 
is going to be started at the commissioning of the LHC, which will happen with proton 
beams at low luminosities . 

This review paper is organized as follows. After a section describing the ALICE detec- 
tor, we will review the features of ALICE operations with pp collisions at the LHC, and 
the statistics and triggers required to accomplish its physics programme. Several physics 
topics will be addressed, but special emphasis will be given to the soft physics programme 
(event characterization, strange particle and resonance production, particle correlations, 
event-by-event fluctuations and baryon asymmetries measurements). Then, the ALICE 
capabilities in measuring some diffractive processes and its potentialities in the study of 
hard processes (jet and photon physics) will be presented, to conclude with some hints 
of possible studies of exotic processes (like mini black holes eventually produced by large 
extra dimensions). 

The physics programme will include also measurements of heavy-flavoured mesons 
(open charm and beauty) and of quarkonia states, both in the central detector and in the 
forward muon spectrometer. However these studies will not be discussed here, since they 
are already reported in other contributions included in these proceedings ([HIS]). 

1.2 ALICE detector overview 

ALICE, whose setup is shown in Fig. is a general-purpose experiment whose detectors 
measure and identify mid-rapidity hadrons, leptons and photons produced in an inter- 
action. A unique design, with very different optimisation than the one selected for the 
dedicated pp experiments at LHC, has been adopted for ALICE. 




Figure 1.1: The ALICE experiment in its final layout. 

This results from the requirements to track and identify particles from very low 
(~100 MeV/c) up to fairly high (~100 GeV/c) Pt, to reconstruct short-lived particles 
such as hyperons, D and B mesons, and to perform these tasks even in a heavy-ion colli- 
sion environment, with large charged-particle multiplicities. 
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Theoretically founded predictions for the multiplicity in central Pb-Pb collisions at 
the LHC range at present from 1000 to 4000 charged particles per rapidity unit at mid- 
rapidity, while extrapolations from RHIC data point at values of about 1500. The ALICE 
detectors are designed to cope with multiplicities up to 8000 charged particles per rapidity 
unit, a value which ensures a comfortable safety margin. 

The detection and identification of muons are performed with a dedicated spectrometer, 
including a large warm dipole magnet and covering a domain of large rapiditie^ (—4.0 < 
r? < -2.4). 

Hadrons, electrons and photons are detected and identified inside the central barrel, a 
complex system of detectors immersed in a moderate (0.5 T) magnetic field provided by 
the solenoid of the former L3 experiment. 

Tracking of charged particles is performed by a set of four concentric detectors: the 
Inner Tracking System (ITS), consisting of six cylindrical layers of silicon detectors, a large- 
volume Time-Projection Chamber (TPC), a high-granularity Transition-Radiation Detec- 
tor (TRD), and a high-resolution array of Multi-gap Resistive Plate Chambers (TOF). 
These detectors allow global reconstruction of particle momenta in the central pseudo- 
rapidity range \rj\ < 0.9 (with good momentum resolution up to Pt ~ 100 GeV/c), and 
particle identification is performed by measuring energy loss in the ITS and in the TPC, 
transition radiation in the TRD, and time of flight with the TOF. 

However, in the case of pp collisions, the lower particle density allows to increase the 
TPC acceptance by considering also tracks with only a partial path through the TPC, i.e. 
ending in the readout chambers; in that case the pseudorapidity coverage can be enlarged 
up to \r]\ < 1.5, with a lower momentum resolution. 

Two additional detectors provide particle identification at central rapidity over a lim- 
ited acceptance: the High-Momentum Particle Identification Detector (HMPID), that is 
an array of Ring Imaging Cherenkov counters dedicated to the identification of hadrons 
with Pt > I GeV/c, and a crystal Photon Spectrometer (PHOS) to detect electromagnetic 
particles and provide photon and neutral meson identification. 

Additional detectors located at large rapidities, on both sides of the central barrel, 
complete the central detection system to characterise the event on a wider rapidity range 
or to provide interaction triggers. The measurement of charged particle and photon mul- 
tiplicity is performed respectively by the Forward Multiplicity Detector (FMD) (over the 
intervals —3.4 < rj < —1.7 and 1.7 < i] < 5.1) and by the Photon Multiplicity Detector 
(PMD) (over the range 2.3 < r] < 3.5). The VO and TO detectors, designed for triggering 
purposes, have an acceptance covering a rather narrow domain at large rapidities, whereas 
a set of four Zero-Degree Calorimeters (ZDC) will measure spectator nucleons in heavy-ion 
collisions and leading particles in pp collisions around beams' rapidity. 

Finally, in order to complete the ALICE capabilities in jet studies, a large lead- 
scintillator electromagnetic calorimeter (EMCal) [6, 7[ will be located between the TOF 
and the L3 magnetic coils, adjacent to HMPID and opposite to PHOS. In its final config- 
uration, the EMCal will have a central acceptance in pseudorapidity of \ri\ < 0.7, with a 
coverage of 180° in azimuth, and an energy resolution of AE/E = IOVo/^/E. It will be 
optimized for the detection of high-p^ photons, neutral pions and electrons and, together 
with the central tracking detectors, it will improve the jet energy resolution. 

The charged-particle multiplicity and the dN^h/d-V distribution will constitute the first 
basic observable which will be measured in ALICE, both for pp and Pb-Pb collisions at the 
LHC. In the central region the best performance in these measurements will be obtained 
with the Silicon Pixel Detector (SPD), the first two layers of the ITS, with approximate 
radii of 3.9 and 7.6 cm. A simple algorithm can be used in the SPD to measure multiplicity 

^In ALICE the z-axis is parallel to the mean beam direction, pointing in the direction opposite to the 
muon spectrometer 
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in a robust way by using 'tracklets', defined by the association of clusters of hits in two 
different SPD layers through a straight line pointing to the primary interaction vertex, 
assumed to be known |8J. The limits of the geometrical acceptance for an event with 
primary vertex at the center of the detector are \ri\ < 2 for clusters measured on the 
first SPD layer and \r]\ < 1.5 for tracklets measured with both SPD layers. However, 
the effective acceptance is larger, due to the longitudinal spread of the interaction vertex 
position, and its limits extend up to about \r}\ < 2 for the multiplicity estimate with 
tracklets [9]. 

Therefore, by considering the partial overlap between the rj ranges covered by the 
SPD (-2 < r? < 2) and by the FMD (-3.4 < tj < -1.7 and 1.7 < 7/ < 5.1), it follows 
that the pseudorapidity range covered by the ALICE experiment for the charged-particle 
multiplicity and the dNch/drj measurements spans over about 8 pseudorapidity units. 

1.3 ALICE operation with pp collisions 

The proton-proton programme of ALICE will start already during the phase of commis- 
sioning of the LHC, when the luminosity will be low (L < 10^^ cm-'^s"^). This time will 
be a privileged period for ALICE to measure pp collisions, because there will be only a 
small pile-up in its slowest detectors and a low level of beam background [3j . 

However, when higher luminosities will be delivered by the LHC, a limiting factor for 
ALICE will be given by the readout of its detectors, essentially by the TPC, which is the 
slowest detector with its drift time of 88 /xs, during which additional collisions may occur, 
causing several superimposed events (pile- up). 

From the point of view of track reconstruction this would not be a problem, since the 
piled-up interactions in the TPC will keep a regular pattern with virtual vertices shifted 
along the drift direction. This can be tolerated, although at the price of heavier tracking 
and larger data volume for the same physics information, at least up to L = 3 x 10^^ 

— 2 —1 

cm s . 

At this luminosity the interaction rate amounts to about 200 kHz, assuming that the 
total inelastic pp cross section is 70 mb. The TPC records tracks from interactions which 
have occurred during the time interval 88 fj,s before and after the triggered bunch crossing. 
Hence on average 40 events will pile-up during the drift time of the TPC, before and after 
the trigger. However, on average only half of the tracks will be recorded, due to the fact 
that the other half will be emitted outside the acceptance. Therefore the total data volume 
will correspond only to the equivalent of 20 complete events. The charged-particle density 
at mid-rapidity in pp collisions at the nominal LHC centre-of-mass energy of ^/s = 14 
TeV is expected to be about 7 particles per unit of pseudorapidity, resulting in a total of 
~ 250(400) charged tracks within the TPC acceptance \r]\ < 0.9 (or within the extended 
acceptance \rj\ < 1.5, when including also tracks with only a short path through the TPC). 
Clearly, tracking under such pile-up conditions is still feasible, since the occupancy is more 
than an order of magnitude below the design value of the TPC. 

For higher luminosity pile-up becomes progressively more difficult to handle, since 
events start to pile-up also in other detectors (silicon drift and silicon strip detectors, 
and then the HMPID). Therefore the luminosity L = 3 x lO^'' cm-^s"^ is the maximum 
that can be tolerated, and in the following we will consider it as a benchmark for ALICE. 
When the LHC will reach its design luminosity {L = 10^^ cm^^s^^) some strategies will be 
needed to record meaningful pp data by reducing the luminosity at the ALICE interaction 
point (e.g. beam defocussing and displacement). 

For the benchmark luminosity L = 3 x lO^'' cm^^s^^ the total pp event size (including 
pile-up and possible electronics noise) is estimated to be of the order of 2.5 MB, without 
any data compression. Thus, running at the foreseen maximum TPC rate of 1 kHz would 
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lead to a total data rate of 2.5 GB/s. However, the online tracking of the High Level Trigger 
will select only tracks belonging to the interesting interaction. This pile-up suppression 
will reduce the event size by at least a factor 10. 

According to current estimates of event sizes and trigger rates, a maximum data rate 
(bandwidth) of the Data Acquistion (DAQ) system of 1.25 GB/s to mass storage, consis- 
tent with the constraints imposed by technology, cost and storage capacity, would provide 
adequate statistics for the full physics programme. This will be possible by using a com- 
bination of increased trigger selectivity, data compression and partial readout. 

The above needs of ALICE for data acquisition are well within the limits of band- 
width to mass storage provided by the central computing facility (TIER-0) of the LHC 
Computing GRID project, that will be installed at CERN. 

1.4 Required statistics and triggers 

An extensive soft hadronic physics programme will be feasible in ALICE using LHC pro- 
ton beams since the machine commissioning phase, when the low luminosity will limit 
the experimental programme to the measurement of large-cross-section processes. This 
programme will include: 

• event characterization, with the measurement of charged particle multiplicity, pseu- 
dorapidity and momentum spectra, and of the (pT)-multiplicity correlation; 

• particle production measurements, i.e. yields and spectra of various identified par- 
ticles, like strange particles (A, S, CI, etc) and resonances (i.e. p, K* and $), and 
baryon-antibaryon asymmetries; 

• particle correlations (i.e HBT interferometry and forward-backward correlations) 
and event-by-event fluctuations. 

The soft hadronic physics programme will rely on data samples of minimum-bias trig- 
gered events. 

The statistics needed depends on the observable under study and spans the range from 
a few 10^ to a few 10® events. For a multiplicity measurement, a few 10^ events will give 
a meaningful data sample; an order of magnitude more is needed for particle spectra; 
to study rare hadronic observables (e.g. CI production) we will need a few times 10® pp 
events. Therefore, a statistics of 10^ minimum-bias triggered events will fulfill the whole 
soft hadronic physics programme. 

Since the readout rate of ALICE is limited to 1 kHz by the TPC gating frequency, the 
requirement is to be able to collect the data at the maximum possible rate: 1000 events/s, 
at an average of 100 events/s. In this way, at an average acquisition rate of 100 Hz, the 
required statistics can be collected during one typical year of operation (10^ s) 

However, at the same acquisition rate, a reasonable statistics for different physics topics 
can be collected already in the first few hours, days, or weeks of data taking. For example 
a few minutes will be sufficient to measure pseudorapidity density with 10^ events, while 
a few hours will allow to collect sufficient event statistics for multiplicity studies. 

For all the above outlined soft physics programme ALICE will need a simple minimum- 
bias trigger for inelastic interactions, that will be provided by two of its sub-detectors: the 
Silicon Pixel Detector (the two innermost layers of the ITS), and the VO. 

The basic building blocks of the Silicon Pixel Detector (SPD) are ladders, arranged in 
two concentric layers covering the central pseudorapidity region, and consisting of a 200 
/xm thick silicon sensor bump-bonded to 5 front-end chips. The signals produced by each 
chip are logically combined to form the global fast-OR (GLOB.FO) trigger element. 
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The VO detector is composed of two independent arrays of fast scintillator counters 
located along the beam pipe on each side of the nominal interaction point and at for- 
ward/backward rapidities. Two different trigger elements are built with the logical com- 
bination of the signals from counters on the two sides: VO . OR requires at least one hit in 
one counter on one side, while VO.AND requires at least one hit in one counter on both 
sides. 

The main background to minimum-bias events are beam-gas and beam-halo inter- 
actions. The rate of beam-gas collisions is expected to be much smaller than the rate 
of beam-halo collisions, whose magnitude should be of the same order as proton-proton 
collisions. It has been shown that the structure of beam-halo events is similar to that of 
beam-gas events, the difference being that beam-halo events happen at greater distances 
to the nominal interaction point (more than 20 m). 

The proposed proton-proton minimum bias triggers, that use logical combinations of 
the above outlined trigger elements, result to be sensitive to interactions corresponding 
to ~ 90 % of the total inelastic cross section (and ~ 99 % of the non-diffractive cross 
section), and still reject the majority of beam-gas interactions |10j . 

On the other hand the SPD global fast-OR (GLOB . FO) trigger element can also be used 
to provide a high multiplicity trigger, that will allow to collect enriched statistics in the 
tail of multiplicity distributions. 

As regards the other physics topics (open heavy flavour mesons and quarkonia pro- 
duction; diffractive processes studies; jet and photon physics) they require separate high 
statistics data samples that would need high rates and bandwidth. Dedicated trigger and 
HLT algorithms will significantly improve the event selection and data reduction, and will 
allow to collect data samples of adequate statistics already in one year of data taking. 

Some details on triggers for diffractive processes and jets will be given in following 
dedicated sections. 

1.5 Event characterization 

For the first physics measurements, shortly after the LHC start-up, in order to minimize 
the uncertainty stemming from non-optimal alignment and calibration, a few detectors 
systems will be sufficient: the two inner layers of the ITS (the Silicon Pixel Detector), the 
TPC, and the minimum-bias trigger detectors (VO and TO). Indeed, particle identification 
will have a limited scope during initial runs, since it requires a precise calibration and a 
very good understanding of the detectors. Four measurements of soft hadronic physics 
which can be addressed during the first days of data taking will be outlined in this sec- 
tion: 1) the pseudorapidity density of primary charged particles, 2) the charged particle 
multiplicity distribution, 3) transverse momentum spectra and 4) the correlation of {pt) 
with multiplicity. 

These measurements of global event properties will be discussed in the context of 
previous collider measurements at lower energies and of their theoretical interpretations. 

1.5.1 Pseudorapidity density 

The pseudorapidity density of primary charged particles at mid-rapidity dA^^ch/dr^l^a^g = 
■^l^vnei . (^({(j^j^l({rq^^^Q has been traditionally among the first measurements performed by 
experiments exploring a new energy domain. Indeed this measurement is important since 
it gives general indications on the interplay between hard and soft processes in the overall 
particle production mechanisms, and furthermore it brings important information for the 
tuning of Monte Carlo models. A simple scaling law (~ In s) for the energy dependence 
of particle production at mid-rapidity was predicted by Feynman [TT], but it appeared 
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clearly broken in collisions at the SPS \1'2\ I13| . Indeed, the best fit to the pp and pp 
data, including that from SPS and Tevatron colliders [2], follows a In'^s dependence, 
whose extrapolation at ^/s=14 TeV gives about 6 particles per rapidity unit for non-single 
diffr active interactions. 

A reasonable description of the energy dependence of the charged particle density 
is obtained within the framework of the Quark Gluon String Model (QGSM) [T5], a 
phenomenological model that makes use of very few parameters to describe high-energy 
hadronic interactions. In this model, based on the ideas of Regge theory, the inclusive cross 
sections dach/df] increase at very high energies and at r/ ~ as a power-law ~ (s/sq)^, 
where so = 1 GeV^ and A = ap — 1 is related to the intercept ap of a Pomeron (Regge) 
trajectory. Indeed, with the value ap = 0.12 ±0.02 found from the analysis of a^"^(s) [16] 
it results that the QGSM model reproduces successfully the observed growth of pseudo- 
rapidity distributions with energy |17] . 

Furthermore, the increase with energy of the charged particle density as well as the bulk 
properties of minimum bias events and of underlying event in hard processes are success- 
fully reproduced (up to Tevatron energy) by models assuming the occurrence of multiple 
parton interactions in the same pp collision |18| [T9| I20j . Examples of such models, ex- 
tending the QCD perturbative picture to the soft regime, are implemented in the general 
purpose Monte Carlo programmes PYTHIA j2l], JIMMY j22j, SHERPA l23j and HER- 
WIG++ [23], all of them containing several parameters that must be tuned by comparison 
against available experimental data. On the other hand, another successful description of 
the available data is provided by the Monte Carlo model PHOJET [251 which is based on 
both perturbative QCD and Dual Parton Model. However, the growth of particle density 
predicted by PHOJET is slower than in multiple parton interaction models, and so the 
charged particle density at LHC energy results to be ~ 30% smaller. 

ALICE will measure the dNch/d.r] distribution around mid-rapidity by counting corre- 
lated clusters (tracklets) in the two layers of the SPD (\ri\ < 2), and/or by counting tracks 
in the TPC (up to |?7| = 1.5). At the low multiplicity typical for proton-proton events, 
the occupancy in the highly segmented detectors will be very low, and corrections for geo- 
metrical acceptance, detector inefficiency and background contamination (from secondary 
interactions and feed-down decays) will be applied on track level. A second correction, 
taking into account the bias introduced by the vertex reconstruction inefficiency, will be 
applied on a event-by-event level (see Ref.[26] for more details). 

The measurement can be done with very few events (10^ events will give a statistical 
error of ~ 2 % for bins of Ar] = 0.2, assuming dNch/dr]\^^Q = 6). 

In addition, the measurement of the pseudorapidity distribution can also be performed 
in the forward region (on the pseudorapidity intervals —3.4 < rj < —1.7 and 1.7 < rj < 
5.1), with the Forward Multiplicity Detector, but a complete understanding of secondary 
processes, which are dominant at low angles, is required. 

1.5.2 Multiplicity distribution 

The multiplicity distribution is the probability P„ to produce n primary charged particles 
in a collision. 

At energies below y/s=63 GeV (up to the ISR domain), the multiplicity distribu- 
tions still scale with the mean multiplicity |27j, following an universal function (P„ = 
(n)~"^$(n/ < n >))[28]. For higher energies, starting from the SPS, the KNO-scaling ap- 
pears clearly broken [29] . The peculiarities of the measured multiplicity distributions (as 
the shoulder structure in their shape) have been explained in a multi-component scenario, 
by assuming an increased contribution to particle production from hard processes (jets 
and minijets). Multiplicity distributions are fitted to a weighted superposition of negative 
binomial distributions corresponding to different classes of events (soft and semi-hard) 
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|3m EU E2]- III alternative approaches, the violation of the KNO-scaling is understood 
as an effect of the occurrence of multiple parton interactions [33], or in terms of multi- 
Pomeron exchanges |34j. 




However, the general behaviour of multiplicity distributions in pp collisions in full 
phase space is quite uncertain. For example data at y^=546 GeV from E735 and UA5 
differ by more than a factor of two above Nch ~ 80 [35] . Therefore extrapolations to higher 
energies or to full phase space of distributions measured within limited rapidity intervals 
are affected by rather big inaccuracies. 

Experimentally the multiplicity distribution is not straightforward to extract. The 
detector response matrix, i.e. the probability that a certain true multiplicity gives a 
certain measured multiplicity, can be obtained from detector simulation studies. Using 
this, the true multiplicity spectrum can be estimated from the measured spectrum using 
different unfolding techniques [361 133 EE] • The procedure of measuring the multiplicity 
distribution with the ALICE detector (using the Silicon Pixel Detector of the ITS, as well 
as the full tracking based on the TPC), is thoroughly described in Ref.[_39j 

ALICE reach in multiplicity with the statistics foreseen for the first physics run (> lO'^ 
minimum-bias triggered events) is about 125 (|r/| < 0.9). However, a large statistics of 
high-multiplicity events, with charged-particle rapidity densities at mid-rapidity in the 
range 60-70 (i.e. ten times the mean multiplicity) can be collected by using a high- 
multiplicity trigger based on the SPD Fast- OR trigger circuit. This class of events may 
give access to initial states where new physics such as high-density effects and saturation 
phenomena set in. 

Also, local fluctuations of multiplicity distributions in momentum space and related 
scaling properties (intermittent behaviour) might be a possible signature of a phase tran- 
sition to QGP[30]. This makes it interesting to study such multiplicity fluctuations in pp 
collisions. 

1.5.3 Transverse momentum spectra 

Collider data on charged-particle px spectra have shown that the high yield rises dra- 
matically with the collision energy, due to the increase of the hard processes cross sections 



At high the transverse momentum spectra are well described by LO or NLO pQCD 
calculations, but involving several phenomenological parameters and functions (K-factor, 
parton distribution functions and fragmentation functions) which need an experimental 
input to be determined. At lower p^, where perturbative QCD calculations cannot be 
performed, theoretical foundations of different models are even more insecure. Therefore, 
early measurements of p^ spectrum are important for the tuning of the model parameters 
and for the understanding of the background in the experimental study of rare processes. 
Also, the measurement of the Pt spectrum is important to perform high-p^ hadron suppres- 
sion studies in in heavy-ion collisions, where the proton-proton data is used as reference. 

In ALICE the track reconstruction is performed within the pseudorapidity interval 
|r/| < 0.9 through several steps (see section 5.1.2 of Ref. [2] for a detailed description of the 
procedure). Firstly, track finding and fitting in the TPC are performed from outside inward 
by means of a Kalman filtering algorithm [3?. In the next step, tracks reconstructed in 
the TPC are matched to the outermost ITS layer and followed in the ITS down to the 
innermost pixel layer. As a last step, reconstructed tracks can be back-propagated outward 
in the ITS and in the TPC up to the TRD innermost layer and then followed in the six 
TRD layers, in order to improve the momentum resolution. 

As it was already said before, the TPC acceptance covers the pseudorapidity region 
\rj\ < 0.9, but this range can be extended up to \r]\ ~ 1.5 when analyzing tracks with 
reduced track length and momentum resolution. 
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The Pt spectrum is measured by counting the number of tracks in each bin and 
then correcting for the detector and reconstruction inefficiencies (as a function of r/ and 
Pt). Finally, the Pt distribution is normalized to the number of collisions and corrected 
for the effect of vertex reconstruction inefficiency and trigger bias. 

With an event sample of > lO''' event that could be collected in the first runs ALICE 
could reach p^ > 40 GeV/c. 

1.5.4 Mean transverse momentum versus multiplicity 

The correlation between charged-track (p^) and multiplicity, describing the balance be- 
tween particle production and transverse energy, is known since its first observation by 
UAl [33] , and it has been successively studied at the ISR jH] and Tevatron jJSl |l6] ener- 
gies. The increase of {px) as a function of multiplicity has been also suggested by cosmic 
ray measurements |47j. 




This correlation between {pt) and multiplicity is generally attributed to the onset 
of gluon radiation, and explained in terms of the jet and minijet production increasing 
with energy ||48|. However, CDF data [46j have shown that the rise of the mean pj< with 
multiplicity is also present in events with no jets (soft events). This behaviour is not yet 
satisfactorily explained by any models or Monte Carlo generators (as PYTHIA [21] and 



In ALICE it will be relatively straightforward to obtain the correlation between the 
mean p^ and the charged particle multiplicity, once the multiplicity distribution and the 
Pt spectra have been measured. The Pt cut-off imposed by the detector (~ 100 MeV/c 
for pions and ~ 300 MeV/c for protons) introduces a rather large systematic uncertainty 
on the {pt) estimate. 

Detailed measurements of the {pt) versus multiplicity (eventually in different regions 
in r]-(j) relative to leading-jet direction, as in CDF analyses [50]) will give an insight to jet 
fragmentation processes and to the general Underlying Event structure. 

Another interesting subject for ALICE, due to its powerful particle identification sys- 
tem at low and high p^, will be the correlation between {pj) and multiplicity studied 
separately for pions, kaons and proton/antiprotons. The data collected at Tevatron by 
the E735 experiment [l5] indicate that the correlation has rather different behaviour for 
the three types of particles, especially as regards the proton and antiproton {pr)-, that 
do not appear to saturate at high multiplicity as pions (and maybe also kaons, within 
experimental uncertainties). This is not yet understood in terms of the available hadronic 
models. 

1.6 Strange particle measurements 

There are basically two main motivations for ALICE to measure strange particle produc- 
tion in pp collisions at the LHC centre-of-mass energy of 14 TeV: 1) extending the range 
where strange quark production has been probed in pp collisions; 2) providing a refer- 
ence for the measurement of strangeness production in heavy ion collisions, in view of the 
strangeness enhancement which was observed to set in at the SPS centre-of-mass energy 



Strange and light quark production rates are usually compared by means of several 
observables. The most simple ones are measured particle ratios, like the widespreadly 
used K/vr, that features a slight and remarkably stable increase in pp and pp collisions, 
between = 27 and 1800 GeV. 

However, ideally one would like to extract directly from data the ratio between newly 
produced s and tt, d quarks at hadronization, before hadronic decays take place. A useful 
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way to measure such strangeness content is the so-called "Wroblewski ratio", defined as: 
Xs = — The earliest attempts to determine Xs were done in [51] , on the basis of the 



uu+dd ' 

models in Refs. [52]. In pp collisions and in the ^/s range from 10 GeV to 900 GeV a fairly 
constant value of As 



0.2 has been extracted from the data (see Fig. 1.2) and there is 
no evidence of a rise [53] with increasing energy. This figure also summarizes the results 
from heavy-ion collisions. 
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Figure 1.2: The Wroblewski factor As as a function of ^/s. 



On the other hand, more recent analyses based on statistical models of hadronization 



have had great success in describing experimental data. This is shown in Fig. 1.2 



by the dashed line, which comes from a canonical description using a correlation volume 
of two protons. This correlation volume causes a strangeness reduction as compared to 
heavy-ion collisions, which have a As around 0.43. 

In the case of heavy-ion collisions the parameter of the statistical models are interpreted 
thermodynamically, ascribing a "temperature" and some "thermodynamic potentials" to 
the system. However, it remains unclear as to how such models can successfully describe 
particle production in systems of small volume like those occurring during pp collisions. 
On the other hand it must be remarked that a pp system does not have to be thermal on a 
macroscopic scale to follow statistical emission. The apparently statistical nature of par- 
ticle production observed in pp data could be simply a reflection of the statistical features 
of underlying jet fragmentation, or a result of phase space dominance considerations. 

From Fig. 1.2 one might conclude that the strangeness content in elementary collisions 
will hardly increase with incident energy. However, this is far from being clear. The 
number of produced particles in pp collisions will increase up to values similar to those 
observed in heavy-ion collisions at the SPS. Hence, the volume parameter in the canonical 
description may have to be increased to account for the higher multiplicity. This would 
result in an increasing strangeness content. 

Furthermore, at the LHC energies, the dominance of jet and minijet production will 
raise new questions since the final hadronic yields will originate from two different sources: 
a source reflecting the equilibrated (grand) canonical ensemble (soft physics) and on the 
other hand the fragmentation of jets (hard physics), that differs from the behavior of an 
equilibrated ensemble. By triggering on events with one, two or more jets, a 'chemical 
analysis' of these collisions will be possible. This very new opportunity would allow us to 
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study whether the occurrence of hard processes influences the Underlying Event distribu- 
tions. Particularly interesting in this context is the behaviour of strange and multi-strange 
particles, e.g. the K/vr or O/vr ratio, in combination with extremely hard processes. 

Possible effects of strangeness enhancement might be amplified when selecting events 
with high multiplicity. In this respect ALICE, thanks to the high-multiplicity trigger 
provided by the Silicon Pixel Detector can collect samples enriched in high multiplicity 
events, and so it will reach multiplicities 10 times the mean multiplicity. This study has 
gained special interest recently, when arguments for 'deconfinement' have been advocated 
in pp collisions at Tevatron energies |55j . 

Moreover, several kinematical properties of strange particle production, like the multi- 
plicity density dependence of their yield and their pj. spectra, have been measured in the 
past, up to Tevatron energies, and still await a full theoretical explanation. For example 
the and A spectra recently measured with high statistics and rather large pr cov- 
erage by the STAR experiment [56* in ^/s = 200 GeV pp collisions at the RHIC collider 
have been compared to NLO pQCD calculations with varied factorization scales and frag- 
mentation functions (taken from [57] for and from |58] for A). Although for the 
a reasonable agreement is achieved between the STAR data and the NLO pQCD calcu- 
lations, the comparison is much less favorable for the A. A better agreement is obtained 
by Albino, Kniehl and Kramer in |59] when using a new set of fragmentation functions 
constrained by light-quark flavour-tagged e^e~ data from the OPAL experiment |60| . It 
will prove helpful to perform similar comparisons at LHC energies. 

It was shown by the E735 Collaboration at Tevatron f6TJ that kaons {pt) has stronger 
correlations with the charge multiplicity per unit rapidity than the pions {pt)'- while the 
latter shows a saturation at {dN/drf) ~ 10, the former continues to grow, although slightly 
decreasing its slope; the same behaviour is seen for antiprotons. Since {dN/drf) can be 
related to the energy density or entropy density [BT, this behaviour is certainly relevant 
for quark-gluon plasma searches, besides providing constraints on models attempting to 
describe hadron production processes. ALICE can test this behaviour at much higher 
multiplicity densities and for other identified mesons and baryons carrying strangeness 
quantum numbers. 

Finally, the measurement of higher resonances in pp will be important to obtain the 
respective population. This can be useful as input for the statistical models, but also for 
comparison with what is found in heavy-ion collisions, though in the latter case the yields 
are likely to be changed by the destruction of the resonances following the rescattering in 
the medium. 

For all the reasons discussed above it appears very important to measure strange 
particles over a broad range of transverse momentum in the new regime of LHC energy. 
The ALICE experiment will face this challenge, for both pp and Pb-Pb collisions, thanks to 
the large acceptance and high precision of its tracking apparatus and particle identification 
methods. 

Strange particle can be identified over a wide range in px from the topology of their 
decays ("kinks" for charged kaons and secondary vertices for X^, A, H and O decays) or 
otherwise from invariant mass analyses (for resonance decays). 

The decay pattern of charged kaons into the muonic channel, with one charged daugh- 
ter track (a muon) and one neutral daughter (a f^) which is not observable in the tracking 
detectors, is known as a "kink", as the track of the charged parent (the candidate) 
appears to have a discontinuity at the point of the parent decay. The kink-finding soft- 
ware loops on all charged tracks by applying to them some cuts to look for pairs of tracks 
compatible with the kink topology described above. The reconstruction of the kink topol- 
ogy is a key technique for identifying charged kaons over a momentum range much wider 
than that achieved by combining signals from different detectors (ITS, TPC, TOF and 
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HMPID). Simulation studies have shown that for a total sample of 10^ pp events, in a 
full year of pp data taking at the LHC, a usable statistics of kaons can be obtained up to 
14 GeV/c. However, when exploiting the relativistic rise of the energy loss signal in the 
TPC, the momentum reach can be further on enlarged up to 50 GeV/c. 

Strange particles as ivT^, A, H and 17 decay via weak interactions a few centimeters 
away from the primary vertex, and therefore they can be identified by using topological 
selections. 

In the case of and A the dominant decay channels are — > vr^vr^ and A piT~ . 
The charged tracks of the daughter particles form a characteristic V-shaped pattern known 
as a "VO" , whose identification is performed by pairing oppositely charged particle tracks 
to form VO candidates. Then, a set of geometrical cuts is applied, for example to the 
distance of closest approach (DCA) between the daughter track candidates and to the VO 
pointing angle, in order to reduce the background and to maximize the signal-to-noise 
ratio. More efficient algorithms for VO reconstruction, named "on-the-fly" , i.e. performed 
during track finding, are also under study. 

The identification of the so-called "cascades" (H~ Avr" and 0,^ AK^), goes 
through pairing VO candidates with a single charged track, referred to as the "bachelor" , 
and then using selections on the VO mass and impact parameter, the DCA between the 
VO and the bachelor, the bachelor impact parameter and the cascade pointing angle. 

The reconstruction of secondary vertices (that has been thoroughly investigated in 
section 6.2 of [2]), relies on the primary vertex reconstruction. ALICE shows good per- 
formances for the identification of secondary vertices from strange hadrons decays, both 
in Pb-Pb and pp collisions. In the latter case particle multiplicities are low, and com- 
binatorial background is even lower, so that topological cuts can be loosened in order 
to gather more signal. However, the reconstruction of the primary vertex position in 
the low-multiplicity events produced by pp collisions is affected by a large error, which 
substantially alter the reconstruction efficiency. 



In any case a clear signal is obtained, as can bee seen in Fig, 1.3 that shows in the 
left panel, for the case of the A reconstruction, the signal and background invariant mass 
spectra obtained for pp events generated with PYTHIA 6.214. The right panel shows the 
estimated distribution of reconstructed A versus Pt in the central rapidity range \y\ < 0.8, 
for one year of LHC running (10^ events). 
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Figure 1.3: Invariant mass distributions of A reconstructed in pp events for tight selections 
(left); distribution of reconstructed A as a function of for 10^ pp events at ^/s = 14 TeV 
(right). 

Therefore the simulation (presented in j2] and in [63]) shows that transverse mass 
spectra for A should be measurable up to ~ 11 GeV/c in a full year of pp data taking 
at the LHC. On the other hand the maximum pr reachable for and the multi-strange 
hyperons H and O are 12, 8 and 7 GeV/c respectively. 
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It must be remarked that in all the topological methods described above the single 
particle identification is not required, which makes them especially efficient at intermediate 
and high pj-, where particle identification based on energy loss (in ITS or TPC) and time- 
of-flight (in TOF) measurement fail. Therefore these studies can be performed by using 
only the basic ALICE tracking devices (the TPC and the ITS). 

As regards strange resonance identification, as for example the i^*(892) and the 
<I>(1020) since they decay very early, their daughters are not discernible from other primary 
particles. 

Their main decay modes are K* K^'k~ and <I> — > K'^K~ . Therefore these reso- 
nances are identified via invariant mass reconstruction methods that combine all possible 
pairs of primary daughter candidates. The background is very high since no selection other 
than particle identification or track quality is applied, and can be accurately estimated 
by means of 'like-sign' or 'event mixing' procedures. The kaon and pion identification is 
obtained from the energy loss in the TPC and from the time-of-flight measured in the 
TOF. 

It has been found by preliminary analyses that for both resonances reaching as high 
as 4 GeV/c is not problematic. However, the identification of higher-p^^ resonances should 
be done without using particle identification. 

All the tools to identify strange secondary vertices and resonances will provide first- 
physics observables, and a rather large statistics can be detected within the very first hours 
of LHC run. However, within a larger time scale (like the first full year of pp LHC run), 
the statistics of strange particles reconstructed with ALICE will by far overstep that of 
the previous pp and pp experiments, and will allow several new studies that were barely 
achievable up to now because of statistics, such as the properties of Pt spectra in a range 
of Pt covering soft, intermediate and hard regimes, as wide as possible to understand 
the underlying QCD processes; and to be compared with the phenomena observed in 
nucleus-nucleus collisions at comparable centre-of-mass energy. 

1.7 Baryon measurements 

Studies of baryon production in the central rapidity region of high energy pp collisions 
provide a crucial possibility to test the baryon structure and to establish how the baryon 
number is distributed among the baryon constituents: valence quarks and sea quarks and 
gluons. 

Hadronic processes are described by several models (DPM[6T, QGSM [15], PYTHIA 
|21] ) in terms of color strings stretched between the constituents of the colliding hadrons. 
In the framework of such models the dominant contribution to particle production in pp 
collisions involves diquark-quark string excitations followed by string breaking. The un- 
broken diquark system, playing the role of carrier of baryon number, will take large part of 
the original proton momentum and subsequently fragment into leading baryons, concen- 
trated in the fragmentation region of the colliding protons. Such approach, where baryon 
number transfer over wide rapidity intervals is strongly suppressed, describes successfully 
the bulk of data on leading baryon production. 

However, the observed high yield of protons in central rapidity region observed in 
experiments at the ISR pp collider (in the energy interval y/s = 23 — 63 GeV), cannot 
be explained in the framework of such models, that assume an indivisible diquark. These 
measurements indicate that the baryon number can be transported with high probability 
over a rather large rapidity gap (Ay < 4) [65] . An appreciable baryon stopping is observed, 
with baryons exceeding antibaryons, and in association with higher hadron multiplicities. 

To describe such data other mechanisms of baryon number transfer have been sug- 
gested, following the approach originally introduced by Rossi and Veneziano |66) . They 
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have shown how it is possible to generahze to baryons the successful schemes employed to 
unify gauge, dual and Regge-Gribov theories of mesons. Their results on the topological 
structure of diagrams of processes involving baryons can be rephrased in a dual string 
picture in which the baryon (for Nc = 3) is a Y-shaped object with valence quarks sitting 
at the ends and with a string junction in the middle. Then, it can be assumed that the 
baryon number is carried by valence quarks, or otherwise by the string junction itself, 
which is a non-perturbative configuration of gluon fields. 

In a first approach [67] the baryon number of the incident proton is assumed to be 
transferred to a more central rapidity region through a mechanism by which a valence 
quark is slowed down to the central rapidity region, while a fast spectator diquark is de- 
stroyed. The cross section of the baryon number flow has been estimated using perturba- 
tive QCD calculations: it has been found to depend on the rapidity gap Ay approximately 
as exp~^^ and nicely agrees with the data at ISR energies. Another estimate [68], based on 
the topological approach and Regge phenomenology, and considering also the stopping of 
string junctions in the central rapidity region, finds a similar dependence of single baryon 
stopping cross section on energy and rapidity, in agreement with the ISR data. 

In an alternative approach [69] the baryon number is assumed to be transferred dom- 
inantly by gluons. This mechanism does not attenuate baryon number transfer over large 
rapidity gaps, since the transfer probability is independent of rapidity. 

The HERA data on high-energy photon-proton collisions have offered a unique oppor- 
tunity to study the mechanisms of baryon number transfer. The asymmetry in the e-p 
beam energies made it possible to study baryon production in the photon hemisphere up 
to 8 units of rapidity distance from the leading baryon production region. It has been 
shown in [69] that at such large rapidity intervals the gluonic mechanism give a dominant 
contribution to the baryon number transfer. An experimental observable that is useful to 
distinguish between different baryon production models is the proton to anti-proton yield 
asymmetry Ap = 2 , where Np and Np are the number of protons and anti-protons 

produced in a given rapidity interval. The calculations made in ^9] predicted the asym- 
metry to be as big as about 7 %, which appeared to be in reach with the statistics collected 
by the experiments at the HERA ep collider [70]. However, both the gluonic and valence 
quark exchange mechanisms were estimated in [69] to give about the same asymmetry 
at ?7 = 0, and appeared to explain the HERA data within experimental and theoretical 
uncertainties. On the other hand, it was shown in [71^ that the two mechanisms can be 
discriminated by studying the dependence of the baryon asymmetry on the multiplicity 
of the produced hadrons. Comparison with HERA data from [70J strongly supports the 
assumption that the baryon asymmetry is dominated by the gluonic mechanism, and ex- 
cludes a large contribution of baryon number transfer by valence quarks. Such asymmetry 
reflects the baryon asymmetry of the sea partons in the proton at the very low x values, 
that are reached (down to x ~ 10~^) at HERA. 

More recently, the R = p/p ratio has been measured at the RHIC collider by the 
BRAHMS experiment [72] in pp collisions at y/s = 200 GeV. The introduction of a string 
junction scheme appears to provide a good description of their data over the full coverage 
of < y < 2.9. 

The ALICE detector at the LHC, with its particle identification capabilities and abun- 
dant baryon statistics in the central-rapidity region, (7 • 10^ p, lO'^ A, 2 • 10^ H and 10^ O 
will be recorded with 10^ minimum bias events), is ideally suited to perform baryon flow 
studies. 

Experimental observables that are useful to distinguish between such different models 
are the proton to anti- proton yield ratio R = p/p and their asymmetry Ap = 2^-^^. 
Similar observables can be defined also for A and other identified hyperons, and can be 
studied as a function of particle multiplicity. 
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At the LHC energies, the rapidity gap between incoming protons and central rapidity 
win be 9.6. That would allow the contribution from valence quarks to be probably negligi- 
ble in comparison to that from gluons. On the other hand, within the limited acceptance 
of the ALICE central detectors (|r/| < 0.9) the proton-antiproton asymmetry predicted by 
different baryon flow models (being on the order of 5 % at the LHC), would differ only 
shghtly (a few %). 

A detailed studied [73] has been performed of the systematic errors affecting the asym- 
metry measurement, coming from transport code and material uncertainties, contamina- 
tion from beam-gas events, and from the different quality cuts imposed at the event and 
track level. The estimated upper limit of the systematic error in the anti-proton to proton 
ratio and in the asymmetry Ap is below 1%, that is sufficient to keep these measurements 
at the level of accuracy required at the LHC. 

Such measurements will be also relevant for comparison to heavy ion collisions where 
baryon stopping should be dramatically enhanced. 

Finally, ALICE can also study heavy flavour baryons (A;,, H^, O,},...) which are poorly 
known. Due to the branching ratio (4.7 it 2.8) • 10~^ of the decay channel Af, J/ip A, 
10^ events triggered on J ftp using the TRD detector should produce a few thousands A^. 

1.8 Correlations and fluctuations 

The study of the correlations among particles emitted in hadronic collisions is important 
in order to unveil the properties of the underlying production mechanisms. 

First, the analysis of the two-hadron momentum correlations provides valuable infor- 
mation to constrain the space-time description of the particle production processes. These 
measurements are of great interest both for nucleus-nucleus collisions, where collective ef- 
fects in nuclear matter are studied, and for pp collisions, where they provide clues about 
the nature of hadronization. 

Momentum correlations can be analysed using an interferometric technique that ex- 
tracts space-time information on the particle emitting source by means of a Fourier trans- 
formation of the measured two-particle correlation function. Such technique was initially 
developed in astronomy by Hanbury-Brown and Twiss (HBT) to infer star radii from the 
measurement of a two-photon correlation function. 

Particle correlations arise mainly from quantum statistics effects for identical particles 
and from final state interactions (Coulomb interactions for charged particles and strong 
interactions for all hadrons). The two-particle correlation function C(pi,p2) is defined 
as the ratio of the differential two-particle production cross section to a reference cross 
section which would be observed in the absence of the effects of quantum statistics and 
final state interactions. Therefore, experimentally the two-particle correlation function 
can be obtained from the ratio C{q,K)= A{q,K) /Ji(q,K) , normalized to unity at large q, 
where q is the relative momentum of a pair, and K is the average pair momentum: the 
numerator A{q,K) is the the distribution of the relative momentum for pairs of particles 
in the same event, whereas the denominator is the same distribution for pairs of particles 
in different events. 

In order to extract information from the measured correlation function about the space- 
time geometry of the particle emitting source, it is generally assumed that the source 
distribution can be parameterised as a Gaussian. Simple analyses generally reconstruct 
source size in one dimension, thus providing a correlation function of the form: 

C = l + Aexp(-i?2g2) 

where qt is the component of g = pi — p2 normal to pi + p2, and R and A (chaoticity) 
are the parameters related to the source size and to the strength of the correlation effect. 



20 



respectively. 

Most recent analyses have done a major effort in reconstructing the 3-dimensional 
source shape using the so-called Pratt-Bertsch cartesian parameterization to decompose 
the relative momentum vector of a pair q into a longitudinal direction qi along the beam 
axis, an outward direction qo transverse to the pair direction, and a sideward direction 
perpendicular to those two. Then, according to some given assumptions, the correlation 
function takes the simple form: 



Thus, three HBT parameters {Ro, Rs and Ri) are extracted from the data, containing 
information about the space-time extent of the particle emitting source in the out, side 
and long directions. 

A pronounced dependence of HBT parameters on charged particle multiplicity in 
hadron-hadron collisions has been observed by several experiments: UAl [71] and E735 
|75) in pp collisions at respectively \/s =630 GeV and 1.8 TeV, and more recently STAR 
|76] in pp collisions at ^/s =200 GeV. Furthermore pion HBT results from the STAR 
experiment |,76j have shown a transverse mass dependence (my = + rn^) of the HBT 
radii which is surprisingly independent of collision system (pp or nucleus-nucleus colli- 
sions), and very similar to the rriT dependence measured by NA22 [77] in hadron-hadron 
reactions at the lower CERN SPS energies. Since the vtit dependence of the HBT radii 
in heavy-ion collisions is usually attributed to the collective flow of a bulk system, results 
observed for hadron-hadron collisions could suggest that also in this case a thermalized 
bulk system undergoing hydrodynamical expansion is generated j78| . However, alternative 
scenarios have been proposed to explain the observed rriT dependence, and the question 
is still open. 

As shown in sect. 6.3 of Ref.[2], all such studies of particle interferometry can be 
performed with good accuracy also in ALICE, thanks to its accurate tracking devices and 
its low Pj. cutoff, in order to test different theoretical models of particle production in pp 
collisions in the TeV region. 

Since the expected source sizes in pp collisions are of the order of 1-2 fm, two-particle 
correlation functions are much wider than those obtained with nucleus-nucleus collisions. 
This, together with the smaller track density, makes in principle the momentum correlation 
analysis easier in pp than in heavy-ion collisions. On the other hand in pp collisions 
additional correlations come from the fact that at LHC energies a substantial fraction of 
the particles is produced inside jets. Therefore, additional analysis cuts are needed to 
prevent the merging of close track pairs. Predictions for two pion correlations in y/s =14 
TeV collisions are provided in Ref.[79j, where it is shown how it might be possible to 
obtain information on the hadronization time in these collisions. 

Besides momentum correlations, other kinds of correlations among final state particles 
are important, in order to reveal the properties of the underlying production mechanisms. 

First, we can consider two-particle correlations in rapidity: if Cn{rji,r]2) = P2 (^i' V2) — 
Piirji) Pi{r]2) is the semi-inclusive two-particle correlation function for events with a fixed 
multiplicity n, written in terms of the single and two-particle densities, then we can define 
an inclusive correlation function Cs(7?i,?72) in terms of the Cn(^i,f?2) as: 



with Pn the probability to find an event with the multiplicity n. As it was shown by the 
UA5 data at ^/s = 200 and 900 GeV, Cs{r]i, 7/2) is sharply peaked at 771 = 772, and for this 
reason it is usually referred to as a "short-range" correlation function. The qualitative 
shape of such correlation is well reproduced by a model [SO] where the equation of the 



C{q,K) = 1 + XeM-Rl{K)ql - Rl{K)ql - Rj{K)qh 




n 
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perturbative Pomeron results from the summation the of all orders of pQCD in the Leading 
Log Approximation (LLA). 

On the other hand in hadron-hadron collisions clear evidence exists for strong long- 
range correlations between the charged particles produced into opposite (forward and 
backward) c.m.s. hemispheres of a collision, and also between the particles produced 
in two rapidity bins separated by a wide rapidity gap Ar/. For pp and pp collisions the 
forward-backward multiplicity correlation coefficient increases logarithmically with energy 
over a large energy interval (from ISR to Tevatron energies). On the other hand, such 
dynamical correlations are absent or quite small in e+e^ collisions, up to LEP energies. 
Several attempts have been made to explain such correlations within the framework of 
hadronic string models, or by assuming that particles are produced through the decay of 
ancestors bodies named clusters (or clans) [HU [HH ESI EH ES] , but the exact dynamical 
origin of such correlations still seems unclear. Therefore the study of the forward-backward 
multiplicity correlations represents a useful tool to test any model of hadron production, 
also in the LHC energy domain [86]. On such respect the ALICE experiment is well 
designed for such studies, since its Forward Multiplicity Detector extends the charged 
particle multiplicity measurement from the pseudorapidity interval — 2 < r/ < 2 covered 
by the SPD to the range —3.4 < rj < 5.1, thus allowing to study the multiplicity correlation 
between largely separated rapidity bins. 

Another interesting subject is the study of two-particle correlations in azimuthal angle 
(f), initially proposed by Wang [87] as a method to understand the role of minijets in 
high energy hadronic interactions. It was argued that calculating C{(f)i, ^2) for samples of 
particles with above a given p^*, the influence of the underlying soft processes could 
be reduced: the higher the Pj^^, the more the correlation should look like the profile of 
high-pr jets. 

New analysis approaches have been developed recently by STAR collaboration |88| 
[89] to study two-particle correlations in 200 GeV pp (and nucleus-nucleus) collisions at 
RHIC. By looking at the two-particle correlations on transverse rapidity yx = ln[(mj' -|- 
PT)/'m'TT]i pseudorapidity ??, azimuth and on the angular difference variables ?7a = i]i~V2 
and = — </>2, they found that low-Q^ parton fragments (minijets) dominate the 
correlation structure observed both in pp and in nucleus-nucleus collisions. In particular 
they found that at low the fragmentation process in pp differs markedly from the pQCD 
factorization picture, the 'jet cone' being strongly elongated in the azimuth direction. 

Additional valuable information on the collision dynamics may be obtained in the 
event-by-event studies of the correlations between various observables measured in sepa- 
rated rapidity intervals. Model-independent detailed experimental information on long- 
range correlations between such observables as charge, net charge, strangeness, multiplicity 
and transverse momentum of specific type particles could be a powerful tool to discrimi- 
nate theoretical reaction mechanisms. 

On the other hand the experimental studies of the correlations in small domains of 
the phase space have to cope with the problem of the local fluctuation of the produced 
hadrons and, more generally, of the experimental observables. Indeed, large concentrations 
of particles in small pseudorapidity intervals for single events have been seen in JACEE 
cosmic ray experiment |90j, and in the fixed-target experiment NA22 [91 . A possible 
explanation of these spikes was related to an underlying intermittent behaviour, i.e. to 
the guess that there exists a correlation at all scales which implies a power-law dependence 
of the so-called "normalized factorial moments" of the multiplicity distribution on the size 
of the phase-space bins. If the above mentioned scaling law will be confirmed in the LHC 
energy domain a new horizon will be opened on self-similar cascading structure and fractal 
properties of hadron-hadron collisions. 

The ALICE experiment is well designed for correlation studies, as well for the event-by- 
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event measurement of several observables. Charged particle measurement in the central 
region is given by the combination of the ITS, TPC and TOF detectors, that provides 
momenta and particle identification of hadrons. The charged particle multiplicity mea- 
surement in the pixel-detector of the ITS can be measured up to rj = ±2 , and the FMD 
extends this range to —3.4 < rj < 5.1. In the central rapidity region the calorimeter 
PHOS with a rather limited coverage provides photon multiplicity and photon momenta, 
whereas FMD is designed for photon multiplicity in the high particle density region of 
forward rapidity (2.3 < r] < 3.5). Therefore the combination of the information coming 
from these detectors provides an excellent opportunity to study particle correlations as 
well event-by-event physics and fluctuation phenomena at the LHC energies. More details 
can be found in sect. 6.5 of Ref. [2i. 

1.9 Diffractive physics 

Diffractive reactions in proton-proton collisions are characterised by the presence of ra- 
pidity gaps and by forward scattered protons. Experimentally, a diffractive trigger can 
therefore be defined by the tagging of the forward proton or by the detection of rapidity 
gaps. 

In ALICE, in absence of Roman pot detectors for proton tagging, a diffractive double- 
gap Level-0 trigger can be defined by requiring little or no activity in the forward detectors 
(as the VO), and a low multiplicity in the Silicon Fixel Detector (SFD) of the central barrel 
|92] . However, in defining a LO diffractive trigger, also the signals of other fast detectors 
of the central barrel must be used, especially the TRD, that is put in sleep-mode after the 
readout of an event. Therefore, since the SFD signal would not be in time for the wake-up 
call of the TRD, the VO signals are firstly trasferred to the TRD pre-trigger system, where 
a wake-up call signal is generated by using the information provided by the time-of-flight 
(TOF) array. The output of such a trigger unit is fast enough to reach the ALICE central 
trigger processor well before the time limit for LO decision. 

The acceptance and segmentation in pseudorapidity of the VO detectors allow to select 
a gap width of approximately 3 and 4 pseudorapidity units beyond \ri\ = 2 on the two 
sides, in steps of half a unit. Then, the high-level software trigger (HLT), having access to 
the full information coming from the central tracking detectors, can enlarge the rapidity 
gap to the range —3.7 < r] < —0.9 and 0.9 < r] < 5. 

Furthermore, the information of the zero-degree calorimeter (ZDC) can be used in 
the high-level trigger to identify different diffractive event classes. Events of the type 
PP ^ pN*X (where X denotes a centrally produced diffractive state), are characterised 
by a signal in either the two ZDC calorimeters, whereas events pp — > N*N*X present a 
signal in the calorimeters of both sides. 

Therefore the geometry of the ALICE experiment is suited for measuring a centrally 
produced diffractive state with a rapidity gap on either side. Such topology results from 
double-Fomeron exchange with subsequent hadronization of the central state. It is ex- 
pected that such events show markedly different characteristics as compared to inelastic 
minimum bias events. For example mean transverse momenta of secondary particles are 
expected to be larger, and also the K/tt ratio is expected to be enhanced. 

A soft/hard scale can also be defined according to whether the pr of the secondaries 
is smaller or larger than some threshold value pthr- The invariant-mass differential cross 
section is thought to follow a power law: ~ A study of the exponent A as a 

function of the threshold value pthr can reveal the contribution from soft/hard exchanges. 
Such analysis can be carried out as a function of rapidity gap width. 

Signatures of Odderon exchanges can be searched for in exclusive reactions where, 
besides a photon, an Odderon (a color singlet with negative C-parity), can alternatively 
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be exchanged. For example, diffractively produced C-odd states such as vector mesons 
(j), J T can result from photon-Pomeron or Odderon-Pomeron exchanges. Any excess 
beyond the photon contribution would be an indication of Odderon exchange. Estimates 
of cross sections for diffractively produced J/ ^ in pp collisions at LHC energies [93] result 
to be at a level that in 10^ s of ALICE data taking the J/^p could be measured in its e^e~ 
decay channel at a level of 4% statistical uncertainty (see Section 6.7.5 of p] for more 
information on quarkonia detection in the dielectron channel in the ALICE central barrel) . 
Furthermore, a transverse momentum analysis of the J/ip might allow to disentagle the 
Odderon and photon contributions, following their different t-dependence. 

Finally, diffractive heavy quark photoproduction, characterised by two rapidity gaps 
in the final state, represents an interesting probe to look for gluon saturation effects at 
the LHC [21], where the cross sections for diffractive charm and bottom photoproduction 
amount respectively to 6 nb and 0.014 nb |95j . Heavy quarks with two rapidity gaps 
in the final state can, however, also be produced by central exclusive production, i.e. 
two-Pomeron fusion. However, since the two production mechanisms have a different 
t-dependence, a careful analysis of the Pt dependence of the QQ pair might allow to 
disentangle the two contributions. 

1.10 Jet physics 

The measurement of jet production in pp collisions is an important benchmark for under- 
standing the same phenomenon in nucleus-nucleus collisions. The energy loss experienced 
by fast partons in the nuclear medium (through both radiative [96] [971 IMl [99] and col- 
lisional |100| I101| I102| I103| 1104] mechanisms) is expected to induce modifications of the 
properties of the produced jets. This so-called jet quenching has been suggested to be- 
have very differently in cold nuclear matter and in QGP, and has been postulated as a 
tool to probe the properties of this new state of the matter. The strategy is to identify 
these medium-induced modifications that characterise the hot and dense matter in the 
initial stage of a nucleus-nucleus collision, by comparing the cross sections for some jet 
observables in benchmark pp collisions at the same centre-of-mass energy. 

An accurate understanding of jet and individual hadron inclusive production in pp 
collisions is therefore quite important in order that this strategy be successful. In this 
respect, the LHC will open a new kinematic regime, in which the pp collisions involve fea- 
tures which are not well understood yet. Therefore, the ALICE experimental programme 
will also involve specific studies on jet and high-pj. particle production in pp collisions. 

In its original design ALICE can only study charged-particle jets by using the tracking 
detectors of the central barrel part of the experiment, covering the region \ri\ < 0.9. Their 
high-pr capabilities, with a momentum resolution better than 10% at pr = 100 GeV, are 
sufficient for jet identification and reconstruction up to Et — 200 GeV. 

However, the strength of ALICE consists in the possibility of combining these features 
with Iow-Pt tracking and particle identification capabilities, to perform detailed studies of 
jet-structure observables over a wide range of momenta and particle species |105] . 

Furthermore, since the charged-jet energy resolution is severely limited by the amount 
of charged to neutral fluctuations (~ 30%), an electromagnetic calorimeter (EMCal) has 
been designed [Sl[7] to complete the ALICE capabilities at high Ej^. The EMCal covers the 
region \ri\ < 0.7, 60° < (p < 180° and has a design energy resolution of AE/E = 10%/^/E. 
The EMCal will improve the jet energy resolution, increase the selection efficiency and 
further reduce the bias on the jet fragmentation through the measurement of the neutral 
portion of the jet energy. Furthermore, it will add the jet trigger capabilities which are 
needed to record jet enriched data at high Et- 

The low and high transverse momentum tracking capabilities combined with electro- 
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magnetic calorimetry will represent an ideal tool for jet structure studies at the LHC over 
a wide kinematic region of jet energy and associated particle momenta, from the hardest 
down to very soft hadronic fragments. A similar strategy has also been used by the STAR 
collaboration at the RHIC collider to reconstruct jets with an electromagnetic calorimeter 
and a TPC, and then to perform systematic studies of fragmentation functions in inclusive 
jets from pp collisions at ^/s = 200 GeV |106j . 

ALICE will study jet production on a large range, from minijet region (E^ >2 GeV) 
up to high-£^r jets of several hundred GeV. However, the event-by-event jet reconstruction 
will be restricted to relatively high-energy jets, approximately Ej^ >30-40 GeV, whereas 
leading-particle correlation studies will play an important role at Iow-Et- 

Observables of interest for jet studies will include: 1) the semi-hard cross sections, 
measured by counting all events with at least one jet produced above some given Ej'j 2j 
the relative rates of production of 1, 2 and 3 jets as a function of the lower Ej- cutoff; 
3) the double-parton collision cross-section and their distinction from the leading QCD 
2^4 process; 4) the properties of the Underlying Event (UE) in jet events, as it has 
been done extensively by the CDF Collaboration at the Tevatron [50] by examining the 
multiplicity and the spectra of charged tracks in the "transverse" region in ij-cj) space 
with respect to the direction of the leading charged particle jet. 

The jet yield that can be measured with ALICE in a running year (10^ s) has been 
estimated by using the hadronic cross sections calculated at NLO |107j for a cone algorithm 
with R=0.7, and usi ng C TEQ5M p.d.f. and factorization and renormalization scales equal 
to 11 = Et/2. Fig. 1.4 shows the annual jet yield for inclusive jets with Ej^ > E'j,™™ 
produced within the ALICE central barrel fiducial region \r]\ < 0.5 for minimum bias pp 



collisions at the nominal luminosity in the ALICE interaction point L 
(or Lint = 50 pb~^ per year). 



5 X 10'^° cm ^s 




Figure 1.4: Number of jets with E^ > produced per year within \rj\ < 0.5 (ALICE 

central barrel fiducial region) in minimum bias pp collisions at y^=14 TeV with luminosity 
5 X 10^° cm-V^. 



However, the rates estimated as above are production rates, which could only be 
exploited by fast dedicated hardware triggers. The EMCal will provide 7, vr*^ and electrons 
triggers, that can be considered to be jet triggers of a sort, but the resulting sample will 
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be dominated by relatively low-i^^ jets that fragment hard. A more refined selection of 
high-E'T jets requires a jet trigger which sums energy over a finite area of phase space and 
finds the location of the patch with the highest integrated EMCal energy. The expected 
enhancement in statistics due to the EMCal trigger can be estimated by comparing the 
rates to tape of EMCal-triggered observables and equivalent observables using only charged 
tracks in the TPC and simple interaction ('minimum bias') triggers. 

It has been estimated (see Sect. 6.8 in Ref. [2j and Sect. 7.1 in Ref. [6J that the 
EMCal trigger will significantly increase the statistics, by a factor ~ 70 for the vr*^ trigger 
(relative to untriggered charged pion measurements), and by a factor ~ 50 for finite-area 
jets of trigger patch-size At/ x = 0.4 x 0.4. The enhancement will be limited by the 
EMCal acceptance (~ 25% of the TPC acceptance) and by its reduced effective value 
for jet triggers of finite extent in phase space relatively to small-area triggers (7, and 
electrons). However, jet measurements incorporating both the EMCal and tracking have 
significantly better resolution and less bias than jet measurements based solely on charged 
particles. Thus, the EMCal-triggered jets provide more robust measurements even for 
modest trigger enhancements. 

Depending on the setup of the Level- 1 (LI) triggering detectors, the software High- 
Level- Trigger (HLT) will be used to either verify the LI hypothesis or to solely inspect 
events at L2. A very simple online algorithm can run on the nodes of the HLT system to 
online search for jets using the full event information from the central tracking detectors. 
This algorithm is supposed to trigger if it finds at least one charged particle jet with more 
than m GeV in a cone with R=0.7. 

Trigger simulations |108] show that for m=30 GeV data rates in pp can be reduced 
by a factor of 100 relative to LI rates, while keeping 1/5 of the events where E^^™" > 50 
GeV and shghtly more than half of the events with i^j,™" > 100 GeV. 

In case of a HLT running without the help of a jet trigger at LI (as in the running 
scenario before the installation of the EMCal, and neglecting the possibility of a trigger 
provided by the TRD), the yields will drop by a factor of ~ 350. The inspection rate of 
the HLT will be limited to the TPC maximum gating frequency of 1 kHz. The expected 
jet yield accumulated in one year for £'y'"*"= 100 GeV when ALICE is running in this 
configuration, is on the level of 10^ events, that is at the statistical limit for the analysis 
of jet fragmentation function at high-z. 

Therefore a trigger with EMCal will be necessary to collect jet enriched data at 
E-r > 100 GeV and extend the kinematic reach for inclusive jets to above 200 GeV. 
For di-jets, with a trigger jet in the EMCal and the recoiling jet in the TPC acceptance, 
the kinematic reach will be about 170 GeV. 

1.11 Photons 

The study of prompt photons processes, in which a real photon is created in the hard 
scattering of partons, offers possibilities for quantitative and clean tests of perturbative 
QCD (pQCD). 

Lowest-order QCD predicts that prompt photons can be produced directly at a parton 
interaction vertex mainly by two processes: quark-antiquark annihilation (q -|- q ^ 7 -|- g) 
and quark-gluon Compton scattering (g-|-q ^ 7-l-q). Because of the latter process, which 
dominates the photon production in pp collisions, the measurement of prompt photons 
provides a sensitive means to extract information on the gluon momentum distribution 
inside the proton. 

However, an additional source of high-p^ prompt photons is due to the hard bremm- 
strahlung of final state partons (fragmentation photons). The latter is a long-distance pro- 
cess which is not perturbatively calculable, since it emerges from the collinear singularities 
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occurring when a high-pj. parton undergoes a cascade of successive splittings ending up 
with a photon. These singularities can be factorised and absorbed into a parton-to-photon 
fragmentation function which has to be determined experimentally and then included in 
the theoretical calculations. 

The calculations of the production cross section of prompt photons at large px have 
been carried out in the framework of perturbative QCD up to next-to-leading order (NLO) 
accuracy in Og. Their results have been found to describe rather well, within experimental 
errors and theoretical uncertainties, all the prompt photon data collected in pp and pp 
collisions over the last 25 years, both at fixed-target experiments (-^=20-40 GeV) and 
at colliders (^5=63-1800 GeV) [Tmi [TTnl ITTT] . 

Corrections for bremmstrahlung processes, for higher-order QCD diagrams and for 
higher-twist processes have been applied to the theory in recent years. Furthermore, full 
QCD calculations have been implemented up to NLO accuracy in more flexible Monte 
Carlo programmes at the partonic level, that allow to account easily for any kind of ex- 
perimental cut jll2| 1113] . This is particularly important for the analyses of data collected 
in collider experiments, that require isolation criteria on photon candidates in order to 
suppress the huge background of secondary photons coming from hadron decays (mainly 
TT^, 7]) and to reduce the fragmentation component of prompt photon production. 

It has been found that NLO pQCD predictions agree very well also with the most 
recent data collected by DO [HI at the Tevatron Run II (^=1.96 TeV) and by PHENIX 
|115j at the RHIC (-y/s=200 GeV). DO has measured isolated prompt photons in the range 
23 < Pt < 300 GeV/c, whereas PHENIX has collected both inclusive and isolated photon 
data in the range 4 < < 16 GeV/c. 

In a recent phenomenological analysis |116] all available prompt photon cross section 
data, including the most recent data from DO and PHENIX, have been compared with 
NLO pQCD theoretical predictions evaluated at the common scale h=Pt/2. The data 
span two orders of magnitude in energy and there is an agreement over nine orders of 
magnitude in the cross sections between theory and experimental data. 

An exception still comes from the fixed-target experiment E706 |117j , at the Fermilab, 
that measured cross sections several times above theoretical predictions based on NLO 
pQCD calculations, with data and theory differing both in magnitude and shape. Al- 
though resummed calculations |118| accounting for recoil effects due to soft gluon radiation 
have reduced the theoretical uncertainties, E706 data still suggest large non perturbative 
parameters (i.e. an intrinsic |119] ) not required by any other data sets. Data in the 
small-x domain probed by LHC may contribute to clarify this issue, and its relation to 
the recoil resummation. 

Predictions of production rates at LHC, obtained from calculations performed at next- 
to-leading order, still suffer from rather large uncertainties. These uncertainties are associ- 
ated with the choice of renormalization, factorization and fragmentation scales |110l[T20j (of 
the order of 30%). As for the uncertainties associated to the structure functions they are 
expected to be relatively small (almost 10% in the lowest Pt region). 

However, for photon production at the LHC energies, a new kinematic region of small 
X values (x ~ xt = ^p-rj \fs will be explored, especially at low transverse momenta (for 

= 2 GeV/c, this corresponds to x = 3 x 10~^ at = 14 TeV). Therefore one may 
question the reliability of straightforward NLO pQCD calculations in a kinematic domain 
where they have never been tested before, and where recoil corrections and resummed 
calculations may be required. 

Furthermore, a specific feature of photon production at very high energy is related to 
the fact that the bremmstrahlung component becomes large and dominant at small x. The 
bremmstrahlung from a gluon is dominant at not too high pj. (up to 20 GeV), whereas the 
production of photons from final-state quarks still remains important at higher p-^ (about 
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40% of the prompt photon yield at ~ 50 GeV/c jl21] ). However, the bremmstrahlung 
component is not weh under controh in particular the gluon fragmentation function into 
a photon is not sufficiently constrained by previous data, and that results in a factor 
two uncertainty in the prompt photon rate for < 20 GeV. On the other hand the 
uncertainties in the quark fragmentation component might introduce difficulties in the 
calibration of the energy of a jet through the measurement of the energy of the recoiling 
photon in 7-jet events. 



Figure 1.5 shows the predictions at the LHC energies of prompt and decay photon 
spectra |122] . Prompt photon spectra are calculated with NLO pQCD, while for decay 
photons the NLO pQCD estimates, extracted from calculated vr'^ spectra, are compared 
to the predictions of PHOJET/DPMJET yjSJ, an implementation of Dual Parton Model 
which includes soft physics (pomeron exchanges) and semi-hard dynamics. 

The results show that the prompt photon spectrum is dominated by more than an order 
of magnitude by the decay photon spectrum, for which there is an excellent agreeement 
between the predictions of DPMJET and NLO pQCD. 
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Figure 1.5: 
Vs=14 TeV. 



Comparison of prompt and decay photon spectra in pp collisions at 



NLO pQCD predictions for the ratio "/prompt / are in the range 5 x W-^-W-^ for 
Pt < 10 GeV/c, and in the range 10~^-10~^ for 10 < < 100 GeV/c. The ratios, slowly 
increasing with pj., are rather small at LHC energies, and isolation cuts will certainly be 
necessary to reduce the amount of the huge background due to decay photons. 

At low Pt the uncertainty in the photon spectrum is largely due to the choice of the 
gluon fragmentation function into photons, which is hardly constrained by the present 
data |124] . and it goes from a factor up to 2.5 at 3 GeV/c to a factor 10% or less for 
Pt > 20 GeV/c. The uncertainties associated to the structure functions are much smaller, 
of the order of 10% and 2.5% at low and high Pt, respectively. 

Photons will be detected in ALICE by the Photon Spectrometer PHOS (see |125] and 
Sect. 3.9 in Ref. [T| for details), an electromagnetic calorimeter with high resolution but 
limited acceptance (r/ < 0.12, A<I> = 120°). The identification power of prompt photons 
in PHOS is limited by the background created by decay photons (mainly, n^, ^ 7 + 7), 
and is optimal for photons with energy larger than 20 GeV. Below this value, decay 
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and prompt photons cannot be efficiently distinguislied on an event-by-event basis, and 
statistical methods are needed. 

Photon measurements can be performed also with the EMCal, the second electro- 
magnetic calorimeter in ALICE, which has coarser granularity compared to PHOS, but 
a factor eight larger phase space coverage, and so it may allow to extend ALICE photon 
measurements to higher pj.. However, the EMCal capabilities in the photon measurements 
will not be discussed here. 

Two different procedures to select (with PHOS) prompt photons among inclusive pho- 
tons are possible: Shower Shape Analysis (SSA), and Isolation Cut Method (ICM) (see 
|125j and Section 6.9 of Ref. [2j). The former identifies photons by analysing the shape of 
the shower in PHOS, and the latter tags and identifies a photon as prompt if it appears 
isolated, i.e., without charged particles emitted in the same direction. 

To optimize the number of parameters conveying the maximal information about the 
shower topology, a principal component analysis was performed. The values of the two 
principal components, corresponding to the largest eigenvalues, have a Gaussian distribu- 
tion. Low, medium and high purity photon samples can be defined by cutting at three, 
two and one standard deviations, respectively. For medium purity photon samples prompt 
photon identification efficiency is about 85% for pp collisions. The contamination from 
misidentified neutral pions ranges from 0% at E-y = 40 GeV to 40% at = 100 GeV 
and the contamination from misidentified charged hadrons and neutrons ranges from 5% 
at E^ = 20 GeV to 15% at E'-y = 100 GeV. Requiring higher purity photons the rejection 
improves at the cost of an important identification efficiency reduction. To improve the 
situation, additional identification procedures are required. 

The main source of background to the prompt-photon spectrum at high p-^ is due to tt^ 
which at Pt > 40 GeV/c produce single clusters in PHOS. Therefore isolation algorithms 
can be devised, that search for hadrons inside a cone centred around the direction (ryo, V'o) 
of high-pr photon candidates (p^ > 20 GeV/c) identified by PHOS with the SSA method. 

In the case of pp collisions a prompt photon identification probability of 100% and a 
TT^ misidentification probability of 3% were estimated from simulations with cone radius 
R = 0.2. 

Finally, our Monte Carlo simulations indicate that with PHOS the photon spectrum 
in pp (and also in Pb-Pb) collisions can be measured with the statistics of one standard 
year up to about 80-100 GeV/c, with a total systematic error of the order of 20%. 

As already explained in another section of the present report, because of the large cross 
sections available for hard processes at LHC, exclusive jet measurements will be within 
reach. In particular, the measurement of jet topology (jet shape, jet heating, fragmentation 
functions, etc.) will require the identification of jets and the measurement of the parton 
(or jet) energy. A very attractive method of performing these studies is to tag jets with 
prompt photons emitted in the opposite direction to the jet direction. 

A 7-tagging algorithm was developed [126j in ALICE to identify 7-jet events and 
to reconstruct the hadronic jet features. The algorithm was tuned for two experimental 
configurations of ALICE: (i) Charged particles are detected in the central tracking system 
and neutral particles in EMCal (this configuration is labelled as 'TPC-t-EMCal'); (ii) Only 
the central tracking system is available and consequently only charged particles can be 
detected (this configuration is labeled as 'TPC'). 

The jet selection efficiency, defined as the ratio of the number of identified 7-tagged 
jets to the number of prompt photons found in PHOS was calculated using a Monte 
Carlo simulation. The efficiency for the configuration with EMCal is about 30%. For 
the configuration without EMCal we obtained an efficiency of 40-50%, because of: (i) the 
wider selection range, implying a lower identification quality; and {ii) the larger acceptance 
in azimuth of the central tracking system as compared to that of EMCal. 
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Jet fragmentation functions to be measured in a standard year of LHC running (for 
both pp and Pb-Pb collisions) were studied for identified 7-jet events in the Pt range 
from 20 to 100 GeV/c. The fragmentation functions obtained for jet-jet events misiden- 
tified as 7-jet events were also studied. For pp collisions, we obtained a signal (7-jet) to 
background (jet-jet) ratio of about 20 in the configuration without EMCal, and near to 
100% background rejection for the setup with EMCal. 

It is not advisable to use PHOS as a detector of jet neutral particles, because of its 
reduced acceptance. However, we may still consider another approach in which the prompt 
photon is detected in EMCal and jets are detected by the central tracking system. In such 
a setup, considering similar prompt photon identification features in PHOS and EMCal 
and the larger acceptance of the EMCal, the prompt photon detection would be enhanced 
by a factor 7 and consequently the statistical errors would be reduced by a factor 2.6. 
The EMCal granularity provides discrimination via shower shape in the range Pt ~ 
10-30 GeV/c. Due to the low j/n^ ratio, however, a robust 7-jet measurement requires 
additional hadron rejection from isolation cuts. The p^-reach up to pr ~ 30 GeV/c 
matches well the statistical reach in a standard year for 7-jet analysis [6] . 

Therefore an experimental study of the fragmentation function of photon-tagged jets 
(i.e. the distribution of charged hadrons within jets as a function of the variable z, defined 
as z = Pt/E^) will be feasible with EMCal in one year up to px ~ 30 GeV/c. 

Furthermore, in addition to single photon production and photon-jet (or photon- 
hadron) correlations, photon-photon correlations can be studied as well. The LO contri- 
butions to di-photon production are quark-antiquark annihilation (q -|- q — > 7 -|- 7) and 
gluon-gluon scattering (g -|- g ^ 7 -|- 7) . 

Processes where both photons originate from parton fragmentation or where one pho- 
ton is prompt and the other photon comes from the fragmentation of a recoiling parton 
also contribute in LO. In this way, many of the inputs entering the theoretical calculations, 
in particular the fragmentation functions |127j . can be tested. 

However, di-photon final states are not only interesting to perform tests of pQCD but 
they are also signatures for many new physics processes, such as Higgs production at the 
LHC or large extra dimensions. 

A further incentive to study prompt photons in pp collisions in ALICE comes from the 
need to provide a baseline against which medium effects observed in the measurements of 
prompt photons in heavy-ion collisions can be disentangled. 

Medium effects in nucleus-nucleus collisions modify the vacuum production cross sec- 
tions of prompt photons as measured in pp collisions: nuclear shadowing and in-medium 
parton energy loss lead to a suppression of the yield |128] , whereas the intrinsic transverse 
momentum distribution of the partons |129| 1130] and medium-induced photon radiation 
from quark jets |131j enhance the yield. 

On the other hand an expected signature of Quark Gluon Plasma formation in central 
heavy-ion collisions is an increased production of thermal photons, emitted in radiation 
processes, roughly in the Pt range 1 < Pt < 10 GeV/c; therefore it is important to 
understand non-thermal production mechanisms in pp collisions in the same energy range. 

Furthermore, in the case of Pb-Pb collisions, photons emerge almost unaltered from 
the dense strongly-interacting medium and provide a measurement of the original energy 
and direction of the parton emitted in the opposite direction. Therefore medium effects 
will be also identified through modifications of the jet fragmentation function, i.e. by the 
redistribution of the jet energy rather than by reduction of jet rate. A broadening of the 
distribution of the jet-particle momenta perpendicular to jet axis (jt), directly related to 
the colour density of the medium, is also expected |132] . 
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1.12 Exotica: mini black holes from large extra dimensions 



The concept of large extra dimensions provides a way of solving the hierarchy problem 
which concerns the weakness of gravity compared with strong and electro- weak forces. 
The extra space-dimensions, beyond the usual three dimensions, are assumed to be com- 
pactified. i.e. finite, so they are too small to be normally detected. A consequence of large 
extra dimensions is that mini black holes (BH) could exist at the greatly reduced Planck 
mass of around 1 TeV, and thus might be produced at the LHC in pp collisions. 

Quantitative calculations for BH production and detection in the ALICE experiment 
at the LHC have been presented in Ref . [13E[ I134j . In this study the BH event generator 
code CHARYBDIS [IS] has been used, that is coupled to PYTHIA code for parton 
evolution and hadronization. Taking advantage of the large-acceptance and high-precision 
tracking detectors available in ALICE, namely the Inner Tracking System (ITS), the Time 
Projection Chamber (TPC) and the Transition Radiation Detector (TRD), two event-by- 
event hadronic observables were used for BH studies: charged particle multiplicity and 
summed-pi^. 

The conclusions drawn from this study are that under the standard running conditions, 
with a minimum-bias trigger running for four months at the LHC initial luminosity, and 
with a maximum data acquisition rate in ALICE of 100 Hz, only a few BH events could be 
visible above the QCD background and only for a Planck mass Mp = 1 TeV, occurring 
for multiplicity above 200 and summed-pj. above 0.5 TeV/c. 

However it is possible to improve this situation, when applying a simple charged par- 
ticle multiplicity trigger to ALICE events, which is expected to greatly reduce the QCD 
background allowing for significant BH signals to be detected. 

For charged multiplicity, the sensitivity to Mp is raised to 2 TeV and hundreds of 
BH events above background corresponding to this case are expected for multiplicity 
greater than 250. An even better situation occurs for summed Pt distribution, since ten 
of thousands of BH events above background are expected for the Mp = 1 TeV case, and 
tens of BH events even for the Mp = 5 TeV case. 

The signature for BH creation from these simple distributions is seen to be an abrupt 
flattening of their slope, as the transition from QCD to BH-dominated charged particle 
production takes place. 

1.13 Concluding remarks 

We have presented in this document the potentialities of the ALICE detector in the field of 
pp physics. A special emphasis has been given to the minimum-bias pp physics programme, 
that is expected to dominate the start-up of the LHC operation. The importance of such 
a programme has been pointed out both for its intrinsic interest and also as a reference 
system for comparison with the nucleus-nucleus and proton-nucleus studies in ALICE. 
However, it has been shown that significant contributions can be given by ALICE also in 
other pp physics topics. 

The complete ALICE detector has significant advantages compared to other LHC de- 
tectors in pp physics attainable at the low luminosity stage of the LHC, mainly because of 
its low momentum thereshold, good momentum resolution and unique capacity to measure 
and identify a large spectrum of particles, including baryons and strange particles. On the 
other hand the ALICE momentum and angular resolution is at least comparable to the 
one of the other LHC experiments up to 10 GeV/c. Moreover ALICE has the capability 
to measure the transverse momentum of charged particles in the range |?7| < 1.5 and, by 
exploiting both the Silicon Pixel Detector and the Forward Multiplicity Detector, charged 
track multiplicity in the range —3.4 < rj < 5.1. Therefore in many essential ways the 
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ALICE pp programme described here is complementary to those possible with other LHC 
experiments. And especially at the early stage of the LHC operation, ALICE will be able 
to provide a significant contribution to this field. 
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2.15 Introduction 

The ALICE experiment |lj will study nucleus-nucleus collisions at the LHC, with a centre- 
of-mass energy per nucleon-nucleon collision ^/s^ = 5.5 TeV for the Pb-Pb system, in 
order to investigate the properties of QCD matter at energy densities of up to several 
hundred times the density of atomic nuclei. Under these conditions a deconfined state of 
quarks and gluons is expected to be formed. 

The measurement of open charm and open beauty production allows to investigate 
the mechanisms of heavy-quark production, propagation and hadronization in the hot and 
dense medium formed in high-energy nucleus-nucleus collisions. Of particular interest is 
the study of the effects of parton energy loss on c and b quarks. Believed to be at the 
origin of the jet quenching phenomena observed in Au-Au collisions at the Relativistic 
Heavy Ion Collider (RHIC), energy loss is expected to depend on the properties of the 
medium (gluon density and volume) and on the properties of the 'probe' (colour charge and 
mass). The open charm and open beauty cross sections are also needed as a reference to 
measure the effect of the transition to a deconfined phase on the production of quarkonia. 
Heavy-quark production measurements in proton-proton and proton-nucleus collisions 
at the LHC, besides providing the necessary baseline for the study of medium effects in 
nucleus-nucleus collisions, are interesting per se, as a test of QCD in a new energy domain. 



2.16 Heavy-flavour production from pp to Pb— Pb 

Heavy-quark pairs (QQ) are expected to be produced in primary partonic scatterings with 
large virtuality > (2mQ)^ and, thus, on small temporal and spatial scales. At ~ Ar ~ 
1/Q^O.l fm for nic = 1-2 GeV. In nucleus-nucleus reactions, this implies that the initial 
production process is not affected by the presence of the dense medium formed in the 
collision. Given the large virtualities, the baseline production cross sections in nucleon- 
nucleon collisions can be calculated in the framework of perturbative QCD (pQCD). For 
the estimate of baseline production yields in nuclear collisions (to be used for performance 
studies and preparation of the analysis strategies), scaling of the yields with the average 
number (A'coii) of inelastic nucleon-nucleon collisions (binary scaling) is usually assumed: 

d'^?A{pA)/dPtdy = (iVcou) X d^ATpQ /dptdy . (2.1) 
The expected cc and bb production yields for different collision systems at the LHC 



are reported in the first line of Table 2.1 [2j. These numbers, assumed as the baseline for 
ALICE simulation studies, are obtained from pQCD calculations at NLO |3j, including 
the nuclear modification of the parton distribution functions (PDFs) [3] in the Pb nucleus 
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Table 2.1: Expected QQ yields per event at the LHC, from NLO pQCD calculations 
For p-Pb and Pb-Pb, Modification of the PDFs in nuclei is taken into account and A'coU 
scaling is assumed. 



colliding system 


PP 


p-Pb 


Pb-Pb 




14 TeV 


8.8 TeV 


5.5 TeV 


centrality 




min. bias 


0-5% (j^'i^^ 


cc pairs 


0.16 


0.78 


115 


bb pairs 


0.0072 


0.029 


4.6 



(details on the choice of pQCD parameter values and PDF sets can be found in [2]). 
Note that the predicted yields have large uncertainties, of about a factor 2, estimated 
by varying the values of the calculation parameters. An illustration of the theoretical 



uncertainty bands for the D and B meson cross sections will be shown in section 2.17 
along with the expected sensitivity of the ALICE experiment. 

2.17 Heavy-flavour detection in ALICE 

The ALICE experimental setup, described in detail in [U |5], was designed in order to 
allow the detection of D and B mesons in the high-multiplicity environment of central 
Pb-Pb collisions at LHC energy, where a few thousand charged particles might be pro- 
duced per unit of rapidity. The heavy-flavour capability of the ALICE detector is provided 
by: 

• Tracking system; the Inner Tracking System (ITS), the Time Projection Chamber 
(TPC) and the Transition Radiation Detector (TRD), embedded in a magnetic field 
of 0.5 T, allow track reconstruction in the pseudorapidity range —0.9 < r/ < 0.9 with 
a momentum resolution better than 2% for pt < 20 GeV/c and a transverse impact 
parameteij^ resolution better than 60 //m for pt > 1 GeV/c (the two innermost layers 
of the ITS are equipped with silicon pixel detectors)]^ 

• Particle identification system; charged hadrons are separated via di?/dj; in the TPC 
and in the ITS and via time-of-flight measurement in the Time Of Flight (TOF) 
detector; electrons are separated from charged hadrons in the dedicated Transition 
Radiation Detector (TRD), and in the TPC; muons are identified in the muon spec- 
trometer covering the pseudo-rapidity range —4 < r] < —2.5 [6]. 

Detailed analyses [7|, based on full simulation of the detector and of the background 
sources, have shown that ALICE has a good potential to carry out a rich heavy-flavour 
Physics programme. In section 2.18| we describe the expected performance for the exclusive 



reconstruction of D'^ K~7r+ decays in pp, p-Pb and Pb-Pb collisions, and the estimated 
sensitivity for the comparison with pQCD predictions, for the pp case. In section 2.19 we 
present the perspectives for the measurement of beauty production in the semi-electronic 
channel in pp collisions. The expected performance for beauty production measurement 



using muons is described in section 2.20 



^The transverse impact parameter, do, is defined as the distance of closest approacli of the track to the 
interaction vertex, in the plane transverse to the beam direction. 

^Note that, for pp collisions, the impact parameter resolution may be slightly worse, due to the larger 
transverse size of the beam at the ALICE interaction point. This is taken into account in the studies 
presented in the following. 
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For all studies a multiplicity of dA'ch/dy = 6000 was assumed for central Pb-Pb col- 
lision^ We report the results corresponding to the expected statistics collected by AL- 
ICE per LHC year: 10'' central (0-5% o"™'^') Pb-Pb events at luminosity i2pb_pb = 
5 X 10^^ cm^^s^^ and 10^ pp events at Cp^^'^^ = 5 x 10^° cm^^s^\ in the barrel de- 
tectors; the muon spectrometer will collect about 40 times larger samples (i.e. 4 x 10^ 
central Pb-Pb events). 

2.18 Measurement of charm production in the D*^ K^tt^ 
channel 

One of the most promising channels for open charm detection is the D*^ K~7r+ decay 
(and its charge conjugate) which has a branching ratio (BR) of about 3.8%. The expected 
production yields (BR x AN/dy at y = 0) for D*^ (and D") mesons decaying in a K^vr^ 
pair in central Pb-Pb (0-5% cr™*^^) at ^Snn = 5.5 TeV, in minimum-bias p-Pb collisions 
at -^/snn = 8.8 TeV and in pp collisions at ^/s = 14 TeV are, in the order, 5.3 x 10~^, 
3.7 X 10~^ and 7.5 x 10~^ per event. 



Figure 2.6 (left) shows a sketch of the decay: the main feature of this topology is the 
presence of two tracks with impact parameters ~ 100 /im. The detection strategy to 
cope with the large combinatorial background from the underlying event is based on: 

1. selection of displaced- vertex topologies, i.e. two tracks with large impact parameters 
and small pointing angle 0p between the D*^ momentum and flight-line (see sketch 



in Fig.|2.6D; 

2. identification of the K track in the TOF detector; 

3. invariant-mass analysis (see pt-integrated distribution in Pb-Pb after selections in 



Fig. 2.6). 



*This value of the multiplicity can be taken as a conservative assumption, since extrapolations based 
on RHIC data predict diVch/dy ~ 2000-3000. 




Invariant Mass [GeV] 



Q 

1.78 1.8 1.82 1.84 1.86 1.88 1.9 1.92 1.94 1.96 

Invariant Mass [GeV] 

Figure 2.6: Schematic representation of the D*^ K~7r"'" decay (left). Kvr invariant-mass 
distribution corresponding to 10^ central Pb-Pb events (right); the background-subtracted 
distribution is shown in the insert. 
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Figure 2.7: Expected relative statistical errors for the measurement in ALICE of the 
production cross section of mesons in 0~5% central Pb-Pb and in pp collisions. 



This strategy was optimized separately for pp, p-Pb and Pb-Pb collisions, as a function 
of the D*^ transverse momentum [7 . 

Figure 2.7 shows the expected relative statistical errors on the measured D*' pt distri- 
bution for pp collisions at 14 TeV (10^ events, i.e. 7 months at nominal pp luminosity 
for ALICE) and central Pb-Pb collisions at 5.5 TeV (10'' events, i.e. 1 month at nominal 
Pb-Pb luminosity). The accessible pt range is 1-20 GeV/c for Pb-Pb and 0.5-20 GeV/c 
for pp (and p-Pb, not shown), with a point-by-point statistical error better than 15-20%. 
The statistical error on the cross section for pt > p™" is estimated to be of about 3% in pp 
and p-Pb {p^^ = 0.5 GeV/c) and of about 7% in central Pb-Pb (^j^" = 1 GeV/c). The 
systematic error (acceptance and efficiency corrections, centrality selection for Pb-Pb) is 
expected to be smaller than 20%. More details are given in [[7j. 

For the case of pp collisions, the experimental errors on the pt-differential cross sec- 
tion are expected to be significantly smaller than the current theoretical uncertainty from 
perturbative QCD calculations. In Fig. 2.8 we superimpose the simulated ALICE mea- 
surement points to the prediction bands from the MNR fixed-order massive calculation [3] 
and from the FONLL fixed-order next-to-leading log calculation [8], |9] . The perturbative 
uncertainty bands were estimated by varying the values of the charm quark mass and of 
the factorization and renormalization scales. The comparison shows that ALICE will be 
able to perform a sensitive test of the pQCD predictions for charm production at LHC 
energy. 



2.19 Measurement of beauty production in the semi-electronic 
decay channel 

The production of open beauty can be studied by detecting the semi-electronic decays 
of beauty hadrons, mostly B mesons. Such decays have a branching ratio of ~ 10% 
(plus 10% from cascade decays b ^ c e, that only populate the low-pt region in the 
electron spectrum). The expected yields (BR x dN/dy at y = 0) for b — > e + X plus 
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Figure 2.8: Sensitivity on d'^a^^' /dptdy, in pp at 14 TeV, compared to NLO pQCD pre- 
dictions from the MNR [3] and FONLL [8j calculations. The inner error bars represent 
the statistical errors, the outer error bars represent the quadratic sum of statistical and 
Pt-dependent systematic errors. A normalization error of 5% is not shown. 



b ^ c(^ e + X) + X' in central Pb-Pb(0-5% o-'"<=i) at ^/s^ = 5.5 TeV and in in pp 
collisions at ^/s = 14 TeV are 1.8 x 10^^ and 2.8 x 10^^ per event, respectively. 

The main sources of background for the signal of beauty-decay electrons are: decays of 
primary D mesons, which have a branching ratio of ~ 10% in the semi-electronic channels. 



and an expected production yield about 20 times larger than B mesons (see Table 2.1); 
decays of light mesons (mainly p, uj, K) and neutral pion Dalitz decays (vr*^ — > 'ye~^e~); 
conversions of photons in the beam pipe or in the inner layers of the ITS; charged pions 
misidentified as electrons. Given that electrons from beauty have an average impact 
parameter do ~ 500 //m and a hard momentum spectrum, it is possible to obtain a high- 
purity sample with a strategy that relies on: 

1. Electron identification with a combined dE/dx (TPC) and transition radiation se- 
lection, which is expected to reduce the pion contamination by a factor of about 10^ 
at low pt- 

2. Impact parameter (do) cut to reject misidentified vr^ and from Dalitz decays 
and 7 conversions (the latter have small impact parameter for pt ~ 1 GeV /c) and 
to reduce the contribution of electrons from charm decays. We have optimized the 
value of the impact parameter cut as a function of the transverse momentum in order 
to minimize the total errors (statistical -|- systematic). The typical value of the cut 
is do > 200 pm. 

The residual contamination of about 10%, mainly accumulated in the low-pt region, of 
electrons from prompt charm decays, from misidentified charged pions and 7-conversion 
electrons can be evaluated and subtracted using a Monte Carlo simulation tuned to repro- 



duce the measured cross sections for pions and D*^ mesons. Figure 2.9 shows the expected 
relative statistical errors for pp collisions at 14 TeV (10^ events, i.e. 7 months at nominal 
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Figure 2.9: Expected relative statistical errors for the measurement in ALICE of the 
production cross section of B-decay electrons in 0-5% central Pb-Pb and in pp collisions. 
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Figure 2.10: Sensitivity on da^{pt > p™")/dy, in pp at 14 TeV, compared to NLO pQCD 
predictions from the MNR [3] and FONLL [8] calculations. Error bars are defined as in 
Fig.lTSi 



LHC luminosity) and central Pb-Pb collisions at 5.5 TeV (10'' events, i.e. 1 month at 
nominal LHC luminosity). 

We tested the possibility to infer the ^^"-differential cross section for beauty mesons, 
da^{pt > pj"™)/dy, from the electron-level cross section using a procedure similar to that 
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developed by the UAl Collaboration [TU]. The method, described in detail in Ref. [7], 
is based on Monte Carlo simulation and it relies on the fact that the B meson decay 
kinematics, measured and studied in several experiments, is well understood. It has been 
shown [7] for the Pb-Pb case that , if electrons with pt> 2 GeV / c are used (below this limit , 
the correlation between the electron and B meson momenta is very poor), the additional 
systematic error is negligible with respect to the systematic uncertainties already present 
at the electron level. 



Figure 2.10 presents the expected ALICE performance for the measurement of the 
p™"-differential cross section of B mesons, Aa^{pt > p["™)/dy vs. pj^*" averaged in the 
range |y| < 1. For illustration of the sensitivity in the comparison to pQCD calculations, 
we report in the same figure the predictions and the theoretical uncertainty bands from 
the perturbative calculations in the MNR [H] and FONLL [SI IS] approaches. It can be seen 
that the expected ALICE performance for 10^ events will provide a meaningful comparison 
with pQCD predictions. 

2.20 Measurement of beauty production in the semi-muonic 
decay channel 

Beauty production can be measured also in the ALICE muon spectrometer, —4 < rj < 
—2.5, analyzing the single-muon pt distribution and the opposite-sign di-muons invariant 
mass distribution [7j. 

The main backgrounds to the 'beauty muon' signal are vr^, and charm decays. 
The cut Pt > 1.5 GeV/c is applied to all reconstructed muons in order to increase the 
signal-to-background ratio. For the opposite-sign di-muons, the residual combinatorial 
background is subtracted using the technique of event- mixing and the resulting distribution 
is subdivided into two samples: the low-mass region, M^+^- < 5 GeV, dominated by di- 
muons originating from a single b quark decay through b c{-^ /i^);U^ (BDsame)) and 
the high-mass region, 5 < M^+^^ < 20 GeV, dominated by bb /U^/i+, with each muon 
coming from a different quark in the pair (BBdiff). Both samples have a background 
from cc — > and a fit is performed to extract the charm- and beauty-component 

yields. The single-muon pt distribution has three components with different slopes: K 
and TT, charm, and beauty decays. The first component is subtracted on the basis of the 
identified hadron spectra measured in the central barrel. Then, a fit technique allows to 
extract a pt distribution of muons from beauty decays. A Monte Carlo procedure, similar 
to that used for semi-electronic decays, allows to extract B-level cross sections for the data 
sets (low-mass , high-mass and pt-binned single-muon distribution), each set 

covering a different B-meson pt > p™™ region. The results for central Pb-Pb collisions 



at ^Snn = 5.5 TeV using only the single muons are shown in Fig. 2.11, Since only 
minimal cuts are applied, the reported statistical errors (represented by the thickness of 
the horizontal bars) are very small and the high-pt reach is excellent. Similar performance, 
in terms of pt coverage, is expected for pp collisions at ^/s = 14 TeV. The main sources of 
systematic errors (vertical bars) are: corrections for acceptance and efficiency, subtraction 
of the background muons from charged pion and kaon decays, and fit procedure to separate 
the beauty and charm components. 

2.21 Conclusions 

We presented the performance of ALICE for the measurement of charm and beauty pro- 
duction in proton-proton collisions at = 14 TeV and central Pb-Pb collisions at 
-y/sNN = 5.5 TeV. For the pp case, these measurements will provide sensitive tests for per- 
turbative QCD in a new energy domain. They wiU also be essential for a comparison with 
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B pj"'" (GeV/c) 

Figure 2.11: Minimum-pt-differential production cross section per nucleon-nucleon colli- 
sion for B mesons with —4 < y < —2.5 in central Pb-Pb collisions, as expected to be 
measured from the single- muon data set. Statistical errors (represented by the thickness 
of the horizontal bars) corresponding to 4 x 10^ events and pt-dependent systematic errors 
(vertical bars) are shown. A normalization error of 10% is not shown. The line indicates 
the input cross section. 



the corresponding measurements in Pb-Pb collisions, for example for the investigation of 
c and b quark in-medium energy loss 

We conclude by comparing, for pp collisions at ^/s = 14 TeV, the ALICE envisaged 
acceptance for heavy-flavour production measurements to that of the other LHC exper- 
iments. Figure 2.12 (from the proceedings of the workshop "HERA and the LHC" |12J ) 
shows schematically the pt vs. rj acceptances for charm (c) and beauty (b) hadrons in 
the four experiments, as expected for one year of running at nominal luminosity (note 
that the value of the luminosity is different for each experiment: 10^^ cm~^s~^ for ATLAS 
and CMS, (2-5) x 10^2 cm-^s"! for LHCb, and 5 x 10^° cm'^s'^ for ALICE). ATLAS 
and CMS have similar acceptances for beauty measurements; the minimum accessible pt 
is larger than for ALICE because of the strong magnetic fields, which in turn, together 
with the high luminosity, allow to cover transverse momenta up to 200-300 GeV/c, and 
because of the larger material budget (amount of material) in the inner tracking detec- 
tors. In terms of acceptance for beauty measurements, ALICE overlaps with ATLAS and 
CMS at central rapidity and with LHCb at forward rapidity. The moderate magnetic field 
allows measurements down to transverse momenta of about 2 GeV / c for B mesons in the 
muon arm and in the barrel, and down to less than 1 GeV/c for D mesons in the barrel. 
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Figure 2.12: Schematic acceptances in transverse momentum and pseudorapidity for open 
heavy flavour hadrons (indicated as "Q-hadrons" ) in the four LHC experiments [12] • The 
high-pt coverages correspond to one year (i.e. 7 months) of running at nominal luminosity 
(see text). [Note that the acceptance of the ALICE muon spectrometer, indicated as 
—4 < r] < —2.5 in the ALICE coordinate system, is reported in the figure as 2.5 < 77 < 4 
to display the overlap with the acceptance of LHCb.] 
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Quarkonia detection with the 
ALICE Muon Spectrometer in pp 
collisions at 14 TeV and PDF 
sensitivity in the low x region 

D. Stocco for the ALICE Collaboration 

The ALICE Muon Spectrometer |lj is a forward detector, with acceptance in the polar 
angle interval 171° < < 178°. It consists of a composite absorber (~ 10Aj„(), made with 
layers of both high and low Z materials, starting 90 cm from the interaction vertex, a 
large dipole magnet with a 0.7 T magnetic field and 10 planes of high-granularity tracking 
stations. A second absorber (~ 7Xint of iron) at the end of the spectrometer and four 
more detector planes are used for muon identification and triggering. The spectrometer is 
shielded throughout its length by a dense absorber tube surrounding the beam pipe. The 
spectrometer was designed in order to detect quarkonia down to ~ in the rapidity 
region —4.0 < y < —2.5. 

The study of quarkonia production in pp collisions presents a twofold interest. On 
the one hand, pp measurements represent a baseline for quarkonia production in heavy- 
ion collisions. On the other, they have an intrinsic interest since they are expected to 
shed light on quarkonia production mechanisms by testing the existing theoretical models 
in an unexplored energy regime. In this respect, the relevant observables are quarkonia 
cross sections and pt distributions. In addition, the rapidity acceptance of the Muon 
Spectrometer for quarkonia will allow access to PDFs at very small x. 

The results of simulation studies of the ALICE Muon Spectrometer physics perfor- 
mance for quarkonia detection in proton-proton collisions at ^/s = 14 TeV are presented. 

The simulation input is provided by the Color Evaporation Model (CEM) [21 [3] pre- 
dictions. In this model the quarkonium production cross section is a measurable fraction 
{Fc) of all QQ pairs below the HH threshold (where H is the lowest mass heavy flavor 
hadron) without any constraints on the color or spin of the final state. The QQ pair then 
neutralizes its color by interaction with the collision-induced color field. At leading order, 
the production cross section of quarkonium state C in an AB collision is: 

fjCEM ^ / ds / dxi<lx2fi/A{xi,^I^)fj/B{x2,H^)cFij{s)5{s-XlX2s) (3.2) 

where A and B can be any hadron or nucleus, ij =qq or gg, aij{s) is the ij QQ 
subprocess cross section and fi/A{xi, n^) is the parton density in the hadron or nucleus. 
Finally, s and s are respectively the hadronic and partonic center of mass energies. The 
results presented here have been obtained with the set of parameters (from [3 ) listed in 
Table [3ll 
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The resulting cross sections in pp collisions, which will be referred to as prompt (and 
include direct production and feed-down from higher mass resonances within the same 
family), are shown in the same table. The values take into account branching ratios in 
the i-i^ n~ channel as well. 







^P' 
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a {ixh) 


3.18 


0.057 


0.028 


0.0069 
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Table 3.2: CEM parameters and resulting cross sections for quarkonia production in pp 
collisions at 14 TeV. Cross sections include feed-down from higher mass resonances and 
branching ratios in muon pairs. The adopted PDF comes from calculations at NLO 
precision by Martin- Roberts-Stirling-Thorne in 1998 (see |4j). 

In addition to prompt and tp', also those from B decay are taken into account in 
this study. These cross sections have been obtained from the open beauty cross section 
using the B J/ip + X and B ^ ip' + X branching ratios. 

The rapidity distributions for prompt production of the different quarkonia states 
are a parameterization of CEM predictions, while the pt distributions are obtained by 
extrapolating to LHC energies those measured by the CDF experiment at ^/s ~ 2 TeV O 
i. 

The invariant mass continuum from semileptonic decay of beauty and charm hadrons 
and from weak decay of pions and kaons was produced with PYTHIA. cc and bb pairs 
were produced with a cross section of 11.2 and 0.51 mb, respectively [3 [8]. 

The so obtained bb pairs were then used in order to get the J/ip and ^' from B decay. 
Since the Muon Spectrometer will not be able to distinguish in a direct way between the 
two main sources of charmonia (prompt and coming from B decay), both contributions 
were summed together to evaluate the expected yields. 

Since the full simulation of a sufficient number of events would require long computing 
times, a fast simulation was performed. Such a method is based on the parameterization 
of the whole spectrometer response at the single muon level. Given a muon of momentum 
p generated at the interaction point with polar and azimuthal angles 9 and ip, the fast 
simulation applies the smearing of the apparatus and gives the reconstructed p' , 9' and (p' 
together with the detection probability for that muon. 

At the trigger level, loose cuts on single muon pt are applied: a low cut for muons from 
charmonia resonances and a high cut for muons from bottomonia ones. Such trigger cuts 
are not sharp but they roughly correspond to pt 1 GeV/c and pt ^ 2 GeV/c, respectively. 



The simulations allow to calculate the global geometrical acceptances for quarkonia, 
integrated over the whole phase space, which are found to be of the order of 4% for both 
J/^ and T. The detection probabilities for the different onium states were computed 
by applying the trigger and tracking response to each muon of the pair. The rapidity 
(transverse momentum) dependencies of the detection probabilities for and T are 



shown in the right (left) panels of Fig. 3.13 These were computed as the ratio between 
the number of detected and generated quarkonia at given y (pt)0 The depletion in the pt 
detection probabilities at low pt is related to the trigger cuts applied. 



^The former ratio was computed by generating quarkonia on the whole pt range, while the latter was 
actually computed by generating quarkonia only in the rapidity interval —4.0 < y < —2.5 covered by the 
spectrometer. 
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Figure 3.13: pt (left) and y (right) detection probabilities for 3/tp (top) and T (bottom). 
The Pt detection probabilities are calculated in the Muon Spectrometer rapidity acceptance 
(-4.0 <y< -2.5). 



Quarkonia yields were computed for a data taking scenario of one year of pp data 
taking (assumed to be equivalent to 10'' s) at a luminosity of 3 x lO^^cm^^s^^ (cf. [9]). 

The resulting fi~^fJ.~ invariant mass distributions are shown in Fig. 3.14 both for the 
J/^ and T regions. As seen, all charmonium and bottomonium states are clearly resolved. 

Together with quarkonia, all sources contributing to the ^^/^~ invariant mass contin- 
uum were taken into account, including muons from correlated and uncorrelated decay of 
cc and bb pairs and from the decay of vr and K. The invariant mass dimuon continuum 
is dominated by correlated sources. 

The quarkonia signal is extracted from the total distribution by means of an interpo- 
lation. For each onium state, a Gaussian function was used for the central part of the 
peak and two more Gaussian with variable width were added to describe the tails. The 
correlated continuum was parametrized with two Gaussian functions with variable width, 
describing the low and high invariant mass regions, respectively. 

The total number of detected J/V' is of the order of 3 x 10^, while the statistics expected 
for T is about two orders of magnitude smaller (3 x 10^). The yields for all quarkonia 
states are summarized in Table 3.3 where the corresponding signal to background ratios 
and significances are also given. 

The obtained statistics will be high enough to allow extracting the dimuon yields per 
bin of transverse momentum and rapidity. To this aim the transverse momentum (rapidity) 
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M,, (GeV/c=) 

Figure 3.14: Opposite-sign dimuon mass spectra in pp collisions at y/s = 14 TeV for a 
running time of 10^ s at a luminosity of 3 x 10'^'^ cm^^s^^. The 3/ip and T mass regions 
are shown in the top and bottom panel respectively. 
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Table 3.3: Expected quarkonia signal and background yields. Numbers refer to an interval 
corresponding to it 1 FWHM around the resonance mass peak. Signal to background ratios 
and significances are also listed. All yields and significances are for a 10'' s running time 
with a luminosity of 3 x lO^'' cm~^s~^. 

of the detected opposite-sign muon pairs was computed and the complete sample of events 
was divided in several bins. For each pt {y) bin, the corresponding dimuon invariant mass 
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distribution was produced. From each of these distributions, the J/V" and T signals were 
extracted by fitting the invariant mass spectrum with the fitting procedure previously 
described. Then the raw number of detected resonances was corrected for the detection 
probability to obtain the differential cross section da/dpt (da/dy). The obtained da/dpt is 
normalized to the rapidity interval (—4.0 < y < —2.5) covered by the Muon Spectrometer. 
The results in Fig. 3.15 (Fig. 3.16) show that the statistical error bars on the measured 
differential cross sections are small, in particular for the due to the high expected 
statistics. 

For completeness the contribution of J/ip from B decay, though not directly measur- 
able, is also shown in the left panels of both figures. 
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Figure 3.15: BR^^da/dpt for (left) and T (right) expected to be measured in a 



data taking time of 10 s at a pp luminosity of 3 x 10 cm~ s~ . The cross section is 
referred to the rapidity window —4.0 < y < —2.5 covered by the Muon Spectrometer. The 
contribution oi J/ip from B decay is also shown in the left panel. 



_ x10'° 

D|>- : 
■°ro0.35- 

m 

0.3 r 
0.25 r 

0.2 ^ 
0.15 r 

0.1 ^ 
0.05 — 



— Total 

— jpsiFromB 



x10» 



D|>- 
EC 2.0 

m 



1.5 



1 



0.5 



H h- 



-3.8 -3.6 -3.4 -3.2 



-2.8 -2.6 



-3.8 -3.6 -3.4 



-3.2 



-2.8 -2.6 



Figure 3.16: BR^fj^da/dy for J/^|J (left) and T (right) measured in a data taking time of 
10'' s at a pp luminosity of 3 x 10^'^ cm^^s^^. The contribution of J/^ from B decay is 
also shown in the left panel. 

Leading order calculations show that in pp collisions at ^/s = 14 TeV, J/^^'s with 
rapidity higher than 3.0 are produced by gluon^ carrying a fraction x of the proton 
momentum lower than 10~^. 



"At high energy the gluon fusion becomes dominant in QQ production, while mechanisms involving 
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In such region, due to a lack of experimental data, the available gluon distribution 
functions rely on extrapolations, thus manifesting a significant disagreement. The feature 
is depicted in Fig. |3.17 showing a comparison between PDF sets calculated at Leading 
Order (LO) precision by different collaborations, namely Martin-Roberts-Stirling-Thorne 
(MRST98 m and MRSTOl [TO]) and the Coordinated Theoretical-Experimental Project 
on QCD (CTEQ5 fU] and CTEQ6 [l2]): the x- values explored by J/ip in the ALICE 
Muon Spectrometer acceptance (in yellow) partially sits in the region of extrapolation. It 
is worth noting that, due to its larger mass, T production is sensitive to x-values larger 
than 10-5. 




Figure 3.17: Comparison of the gluon distributions from MRST and CTEQ. x regions 
probed hy J/ip and T produced in pp collisions at y/s = 14TeV in the rapidity region 
—4.0 < y < —2.5 are shown. 

Starting from the Color Evaporation Model as a guideline, it is possible to show that 
the J/'i/' rapidity distribution, 

ArrCEM rp rim 



dy 



^ J2^" drsa.,{s)f,/A{^lley,^')f^/B{^^e-y,^') (3.3) 



is an observable sensitive to the PDF variation at low x. 

This is done by making the following two approximations: 

• in the elementary cross section aij , only the dominant contribution gg — > QQ is 
taken into account 

• both the elementary cross section agg and the gluon distributions are taken at leading 
order 

The calculation was carried out with different PDFs. The results are summarized by 



the curves in Fig. 3.18 We note that the rapidity distributions shown in this figure are 
normalized by setting equal to unit their integral from -4.0 to -2.5 rapidity units. This 
way of presenting the results emphasizes the fact that different behaviors of the gluon 



distribution functions at low x (see Fig. 3.17) lead to different shapes of the J/ip rapidity 



distribution (see Fig. 3.18) in the interval covered by the muon spectrometer. 



quarks can be neglected. Hence, in the following we will speak about gluon distribution functions instead 
of parton ones. 
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CTEQ5 LO 
MRST98 LO 
CTEQ6 LO 
MRST01 LO 
simulation (+J/\|/ from B) 
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Figure 3.18: Comparison between J/^/; rapidity distributions obtained with four different 
PDF sets (area in detector acceptance is normalized to 1). Simulation results are also 
shown. 
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Such behavior can be easily understood by analyzing Eq. 3.3 The rapidity distribution 
depends on: 

• the product of parton distribution functions 

• the elementary cross section (ygg{s) 

• the ratio F^j^/s (assumed to be constant in the model) 

The last term only affects normalization. Hence the shape of the distribution is related 
only to the first two quantities. However, since the elementary cross section doesn't depend 
explicitly on y and it is not expected to give rise to large variation of magnitude in the 
small range of integration (4mQ < s < Anijj), one can conclude that the shape of the 
distribution is mainly dependent on PDFs. 

The possibility of focusing the study on the shape of the rapidity distribution, disre- 
garding the absolute normalization, is extremely favorable from the experimental point of 
view, since lots of systematic errors (e.g. on luminosity, global acceptance effects, etc.) do 
not enter. 

The next (and last) step of the study presented here is aimed to show that the accuracy 
of data collected with the Muon Spectrometer will be good enough in order to resolve the 
different rapidity distributions. 



A direct comparison with the simulation results in Fig. 3.16 is not possible, since such 
data were obtained by adopting the MRST98 set computed at NLO, while the calculations 
presented in this section concern LO quantities. Hence the 3/^ rapidity distribution was 
re-obtained from simulation after adopting MRST98 LO as input. Such distribution (the 
simulation points in Fig. 



PDF sets (curves in Fig. 



3.18 



3.18 



is then compared with the calculations with different 
It is worth noting that, differently from calculations, 
ude contribution from B decay, which slightly change the 



the simulated data already inc 
shape of the distribution. However the change is limited as it can be seen by comparing 
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simulated data and the calculation of prompt J/V' distribution obtained with the same 
PDF set (red curve of Fig. 3.18). This figure shows that, due to the high statistics, the 
accuracy of the data that are expected to be collected by the ALICE Muon Spectrometer 
will be good enough to allow to discriminate among different gluon distribution functions 
in the region of x < 10~^ (at least in the frame of a leading order analysis). 
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Figure 4.19: 

The Q^-x kinematic plane for the LHC and previous experiments, showing the mass (M) 
and rapidity (y) dependence. 



Parton Densities at the LHC 

A. Tricoli 

4.22 Introduction 

The start up of the LHC machine is now imminent and theorists and experimentahsts are 
converging their efforts to enhance the LHC discovery potentiaL This imphes minimis- 
ing theoretical and experimental uncertainties. Among the theoretical uncertainties the 
knowledge of the proton structure plays a major role: the accurate evaluation of parton 
density functions (PDF's) is vital to provide reliable predictions of new physics signals 
(i.e. Higgs, Supersymmetry, Extra Dimensions etc.) and their background cross sections 
at the LHC. As shown in the contribution by C. Mariotti, E. Migliore and P. Nason , at 
hadron colliders the inclusive cross section for hard production processes is the convolution 
of the cross section at parton level, calculable at fixed order in perturbation theory, and 
the parton densities of the two interacting partons. 

Our knowledge of the proton structure is improving fast thanks to more experimental 
data being available and thanks to more precise and sophisticated theoretical calculations: 
PDF's are nowadays available up to the next-to-next-to leading order (NNLO) in pertur- 
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bative QCD and in recent years they have been also providing uncertainties which take 
into account experimental systematic errors and the correlations between data points that 
enter the global fits. Despite the great improvement on PDF's in recent years, their un- 
certainty dominates many cross section calculations for the LHC. As visible in fig. |4.19 



the LHC will probe kinematic regions in x (parton momentum fraction) and (hard 
scattering scale) never explored before, such as the very high-Q'^ and the very low-x re- 
gions. At low-x the current theoretical formalism (DGLAP) is at the edge of its supposed 
applicability. For the production of Z and W bosons the participating partons have small 
momentum fractions at central rapidity, x ~ 10~^, and in the whole measurable rapidity 
region, \y\ < 2.5, they are within the range 10~^ < x < 0.1. Thus, at the electro-weak 
scale the theoretical predictions for the LHC cross sections are dominated by low-x PDF 
uncertainty. At the TeV scale, where we expect new physics, the interacting partons have 
higher momentum fractions and very high (> 10^ GeV^). Thus, at the TeV scale the 
cross section predictions are dominated by high-x PDF uncertainty and rely on the ex- 
trapolation of the DGLAP equations. In both kinematic regimes the gluon density, which 
is in most regions the less well constrained density function, plays a major part: at low x 
the gluon density dominates the quark and anti-quark densities, at high the interacting 
partons get an important contribution from the sea, which is driven by the gluon density, 
via the g ^ qq splitting process. For a review on hard interactions of quarks and gluons 
at the LHC refer to [I] . 

Past and running experiments, such as HERA, have been providing vital information 
to improve our knowledge of the parton densities, however the broad kinematic region of 
the LHC forces (and offers a unique opportunity to) ATLAS and CMS experiments to use 
their own data to constrain the parton densities, in particular the gluon, in the kinematic 
regions where they are not sufficiently well determined. In section 4.25 it will be shown 
that significant improvement on PDF fits can be made with LHC data. 



4.23 Global fits and error analysis 

Perturbative QCD provides the evolution equations for the parton densities, DGLAP 
equations, but does not provide us with their analytic forms as function of x. The most 
common approach to extrapolate PDF's as function of x and consists in solving the 
DGLAP equations by parameterising the parton densities qi{x) at a fixed scale Qq = 
1 — 7 GeV^, applying assumptions and constraints derived from theory and measurements. 
Then, with the DGLAP equations, we numerically extrapolate the values of qi{x,Q'^) to 
different values of and a global fit of experimental data is performed. For valence 
quarks the parameterisations have usually this behaviour qy ~ x^(l — x)^, whereas for 
the gluon and sea quarks they are of this kind qs{g) ~ x~'^(l — x)"^. However there is 
no unanimous agreement on the parametric functions to use and on the number of free 
parameters. For a review refer to 0. 

Different regions in the x, plane and also different partonic components are probed 
by the available world experimental data. These include DIS data from fixed target 
experiments and HERA, Drell-Yan data, inclusive jet production and W charge asymmetry 
from Tevatron. 

There are various groups who are fitting the proton structure function data, among 
them CTEQ and MRST. Recent PDF sets include in their analyses up-to-date experi- 
mental data and attempt to provide coherent estimates of the uncertainties, including 
experimental correlated systematic errors. The differences between these PDF sets can be 
summarised in three categories: different choices of input data sets, different theoretical 
model assumptions and different error analyses. 

There are many sources of uncertainty which contribute to a global fit uncertainty. 
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Figure 4.20: 

Left: CTEQ6.5M set at ^ m^. Right: comparison between CTEQ6.5M (black) and 
MRST2OO4NLO (red) gluon PDF's and their uncertainties. 

These can have experimental and theoretical origins. The former are related to the data 
errors which enter the fit, the latter are due to the model uncertainties of the theoretical 
framework. The theoretical uncertainties concern both the non-perturbative (parameter- 
isations) and perturbative parts of the calculations: assumptions imposed to limit the 
number of free parameters, higher order truncations in the DGLAP formalism etc. 

The treatment of the experimental uncertainties, especially the correlated systematic 
uncertainties, is a complex subject which is partly still under debate. A modified version 
of the standard method is used to take into account non-Gaussian systematic errors and 
their correlations: ^ + AT^, where AT is the so-called "tolerance", a complicated 
mathematical expression that includes correlated systematic terms [2]. There are then 
two methods to compute the central values of the theoretical PDF parameters and their 
uncertainties: the offset and the Hessian method. In the offset method the correlated 
systematic errors affect only the determination of the PDF uncertainty, not the best fit. 
This method is used for ZEUS PDF's. Conversely in the Hessian method, used by CTEQ 
and MRST groups, the collective effect of the correlated systematic errors has also an 
impact on the best fit. 

For both, the offset and the Hessian methods, the PDF uncertainty is conventionally 
computed along the eigenvectors of the diagonalised covariance or Hessian matrices. The 
number of eigenvectors corresponds to the number of free parameters in the parton density 
parameterisations. Contemporary PDF sets provide a central value PDF set, correspond- 
ing to the best data fit, and two PDF sets for each uncertainty eigenvector, giving the 
upper and lower limit on the uncertainty. Given a PDF set, the upper limit of the PDF 
uncertainty is calculated for a physical observable by adding in quadrature the upward 
displacement eigenvectors, whereas the lower limit by adding in quadrature the downward 
displacement eigenvectors. MRST group has chosen 15 free parameters, leading to 30 



error sets; CTEQ6 has 20 free parameters and 40 error sets. Fig. 4.20 shows CTEQ6.5 fit 



for all parton densities at the scale ~ and its gluon uncertainty compared to the 



MRST2004NLO gluon best fit. 
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4.24 Impact of PDF uncertainty on LHC physics 



The experience from previous experiments teaches that the PDF uncertainties must be 
properly taken into account or features of the SM physics can be misinterpreted as evidence 
of new physics. For example an unexplained discrepancy between data and theory was 
originally found in the Tevatron Run-I jet data, which was subsequently reabsorbed within 
the theoretical uncertainty when a more accurate PDF error analysis was performed. 

G. Polesello's contribution on inclusive jet cross-section has shown that the PDF un- 
certainty is dominating for high Ej- jets over the renormalisation/factorisation scale and 
the experimental energy scale uncertainties: 10% at 1 TeV, 25% at 2 TeV, 60% at 5 TeV. 



Extra dimensions 

In extra dimensions models, if the compactification scale Mc is about few TeV[^ it is 
possible to observe the production of gravitons and Kaluza Klein (KK) excitations at the 
LHC. If gauge bosons can propagate in the extra dimensions, we also expect a violation 
of the SM logarithmic behaviour of the running couplings. In this scenario, if we consider 
the CTEQ6M PDF uncertainty on the di-jet cross-section, we see the extra dimensions 
prediction being absorbed within the SM prediction zone: the high-x gluon uncertainty 
can cause a decrease of the discovery reach from Mc = 5 (10) TeV to Mc < 2 (3) TeV, 
depending on the number of extra dimensions [3] . 



Higgs 

The accurate measurements of the Higgs production cross sections and decay branching 
ratios are crucial to explore all Higgs boson fundamental properties. At the same time, 
we need very precise estimates of the various theoretical uncertainties. 

It is found that the PDF uncertainty can be of the same order of magnitude or even 
higher than the other theoretical uncertainties. In fact the perturbative calculations of 
Higgs production cross section are becoming more stable as higher orders are included, 
leaving the PDF uncertainty as one of the largest contributions to the total theoretical 
uncertainty. For example for the dominant Higgs production channel, gg — > H, the PDF 
uncertainty on gg luminosity, can be larger than the factorisation and renormalisation scale 
uncertainty: in fact the differences in the gg luminosity prediction between MRST2002 and 
Alekhin2002 can be higher than 10% for low Higgs mass scenarios. Furthermore, studying 
the effect of three different PDF sets (i.e. CTEQ6M, MRST2001E and Alekhin2002) with 
their quoted uncertainties, on various Higgs productions channels, we see that the PDF 
uncertainty can be of the order of ^ 10 — 15% on the production cross-section 



High mass Drell-Yan 

Several new physics models predict events with two charged leptons originating from the 
decay of a massive object. A peak in the da/dM distribution is a clean signature of a 
new resonance: the identification and reconstruction of high-mass di-lepton final states 
can be done with high efficiency and the SM background can be small. However the 
shape and normalisation of the predicted observable distributions depend on PDF and its 
uncertainty. 



In fig. 4.21 we see the total CTEQ6.1 uncertainty on the distributions of the recon- 



structed rapidity y and invariant mass Mn of the lepton pair: 40 CTEQ error sets have 



been accounted for, applying the PDF reweighting technique (see sec. 4.26). The uncer- 



tainty is in the range 4 — 7% on both y and Mn up to 1 TeV. Excluding the bins at the 

^In this context the compactification scale is defined as Mc ~ i/Rc where Rc is the compactification 
radius of the extra dimensions on a hypersphere. 
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Figure 4.21: 

CTEQ6.1 uncertainty on distributions of the high-mass di-electron Mu (left and centre) 
and rapidity y (right). Herwig+Jimmy generation and ATLAS full simulation J^. N.B.: 
the drop in the low Mu spectrum is an artifact of the event selection in the Monte Carlo. 



edge of the rapidity distributions, where statistical fluctuations are present, we see that 
the largest PDF uncertainty is at y ~ 0. As explained in [5], a study shows that NLO 
QCD corrections, applied on Monte Carlo (MC) and on PDF, enhance the cross section 
with respect to the LO prediction by 24 — 36%, with the largest NLO corrections at y ~ 0. 
A discrepancy of about 6% is found between MRST-NLO and CTEQ-NLO PDF's. 

4.25 How to constrain PDF at LHC 

Several Standard Model processes are under study to constrain parton densities: 7, W 
and Z boson and inclusive jet production processes are equally important to constrain the 
parton densities and in particular the gluon density in complementary kinematic regions 
(see i). 

In G. Polesello's contribution we appreciate how the LHC jet data can be used to 
better constrain PDF fits: if the experimental systematic uncertainty is under control 
to < 10% level, LHC jet data can significantly contribute to constraining the high-x 
gluon density with 1 fh~^ luminosity. Other studies |i7^ have also shown that the prompt 
photon production process is extremely sensitive to PDF differences and can probe the 
perturbative theory of the gluon at high-x: the discrepancy between MRST2004-NLO, 
CTEQ6.1M and older PDF sets can be of the order of 16 — 18% on the photon 77 and Pt 
distributions. 

Furthermore, the hg — > Zh process is sensitive to the 6-quark content of the proton and the 
LHC predictions for the Z + b cross-section, using different PDF sets, are ±5 — 10% [8]. 

4.25.1 W rapidity distributions 

A few days of LHC running at the nominal low luminosity (10^^cm~^s~^) are sufficient 
to make the statistical uncertainty negligible with respect to the systematic uncertainties 
on W cross section. Among the systematic uncertainties there are experimental and 
theoretical contributions. 

The ATLAS strategy for selecting W bosons consist of identifying an isolated and 
highly energetic lepton, > 25 GeV, and requiring a large amount of missing energy in 
the event due to the neutrino escaping detection, > 25 GeV. The analysis oiW^ eu^ 
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Figure 4.22: 

HERWIG simulations of from decay, with CTEQ6.1M (red), MRST2001 (black) 
and ZEUS-S (green) PDF's and their quoted uncertainties (estimated with the PDF 
reweighting technique as in sec. 4-^6). The top plots are at generator level, the bottom 
plots at ATLFAST detector level. Left fig: e~ (left plots) and (right plots) rapidity 
spectra with NLO-QCD corrections. Right fig: electron- charge asymmetry E 



events fully simulated in the ATLAS detector, in the early data scenario, shows that the 
W boson is a very clean signature: the trigger and the electron off-line identification 
with the electron Et and fJx cuts leave a background contamination dominated by QCD 
events (less than 5%) and W tUt (about 0.5%). If a jet veto cut is added, the QCD 
background can be further reduced to a level of < 1% [9 . Therefore the W sector is an 
ideal environment to study and constrain theoretical and experimental systematics. 



Higher order corrections 

The differential cross section da/dy for W production has been calculated to the NNLO 
order in QCD with an energy scale uncertainty of < 1% |10j. With this level of preci- 
sion in perturbative QCD calculations, the electro- weak (EW) contributions are no more 
negligible. As presented in this workshop, leading order electro-weak contributions with 
multi-photon radiation introduce corrections of the order of few percent on W boson 
cross-sections. The EW corrections, computed by the program HORACE interfaced to 
HERWIG in the q:(0) scheme in the muon channel [Tl], are constant in rapidity and are 
about 3.5% for a cut on the muon transverse momentum of Pt > 25 GeV and can be up 
to 5.2% for loser cuts. The dependence on the muon charge is negligible (up to 0.4% 
for lose Pt cuts) p2]. Considering that these corrections in the muon channel are flat 
in rapidity and negligible on the muon-charge asymmetry, we can state that they do not 
have an impact on the PDF extraction, however they are relevant for luminosity measure- 
ments in order to achieve a precision of 6% or better. The electron channel needs further 
investigation. 



PDF uncertainty on W rapidity distribution. 



From fig. 4.22 we can see the full PDF uncertainties for three different PDF analyses, on the 
rapidity distribution of originating from decays. Their predictions are compatible 
within their uncertainties, which are in the range 4% — 12%, and are dominated by the 
gluon density. 
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In a previous paper [TJ] it is demonstrated that the LHC can improve the current 
constraint on the low-j; gluon parameter {xg{x) ~ a;"'^^) by more than 41% by fitting 
the e"*" and e~ rapidity distributions, if their experimental systematic uncertainties are 
kept under 5% level. 

In the lepton-charge asymmetry Ai = — ^ )/{% ~^ % ) ™ost of the gluon 

uncertainty cancel out leaving the valence up [uv) and down (dy) densities as main 
contributions to the total PDF uncertainty, which is reduced to ~ 5% at ~ 0. However 
a discrepancy of ~ 15% is present at ~ between the MRST2002 and other two PDF's, 
CTEQ6.1M and ZEUS-S [l5]. In fact the MRST PDF's prediction for uy - dy valence 
density is different from the other PDF's and is outside the quoted PDF uncertainty bands. 
This difference in current PDF fits comes from the lack of data on valence quantities at 
such low-x. The LHC can be the first experiment to perform such measurement in the 
kinematic region x ~ 10^^ and = M^. 

A posteriori inclusion of PDF's in NLO calculations. 

The MC computation of QCD final state observables to NLO is a lengthy process. In 
order to study the impact of PDF uncertainties on QCD cross section measurements in a 
faster way and allow for PDF fitting of these quantities, the technique of "a posteriori" 
inclusion of PDF's in NLO calculations has been developed for LHC processes [16] |17] . 
A MC run is used to generate a grid (in xi, X2 and Q) of cross section weights that can 
subsequently be combined with an arbitrary PDF set. This enables the decoupling of 
the lengthy calculation of perturbative MC weights from the convolution with the parton 
densities. Perturbative coefficients for jet (using NLOJET++), W and Z boson (using 
MCFM) production processes can be collected on a grid with an accuracy better than 
0.02%. 



4.26 PDF reweighting of Monte Carlo events 

The computation of the full PDF uncertainty on a physics process is a cumbersome pro- 
cedure. Given one PDF set, such as CTEQ or MRST, it requires the generation of twice 
as many MC samples as the number of free parameters in the global fit. Furthermore one 
error analysis might not be sufficient since, as seen above, there can be large discrepancies 
between the results of different error analyses. 

A PDF reweighting technique has been studied and tested, requiring only one Monte 
Carlo generation with one conventional PDF set |^ [13j |14j. 

This technique has been implemented using hard process parameters of the MC gen- 
eration: flavours {flavi and flav2) and momentum fractions of the interacting partons 
Xfiavii Xfiav2 and the energy scale Q. The PDF set used for the MC generation is named 
PDFi. 

The PDF reweighting technique consists of evaluating, on the event-by-event basis, the 
probability of picking up the same flavoured partons with the same momentum fractions 
Xfiavii Xfiav2i according to a second PDF set, PDF2, at the same energy scale Q, then 
evaluating the following ratio 

Event Weight = /^^^^(^/'"^i ' . fPDF^ixfiav^Q) _ ^^ ^^ 

fpDFi {Xflavi , Q) fPDFi {Xflav2 ) Q) 

After the Event Weight is applied on MC events generated with PDFi, they will 
effectively be distributed according to PDF2. 



*This techniques is not as reliable if the PDF set is as "unconventional" as MRST2003, i.e. the validity 
of its kinematic space is smaller than the one available to the LHC. 
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Figure 4.23: 

Left fig: and rapidity distributions at HERWIG generator level for events gen- 
erated with CTEQ6.1M (dashed lines) and for events generated with MRST2002 (solid 
lines) and their relative differences (at the bottom). The straight lines are the means 
of the points with uncertainty bands. Right fig: same as left hand side plots for events 
generated with CTEQ6.1M (dashed lines) and for events generated with MRST2002 and 
PDF-reweighted with CTEQ6.1 (solid lines). Similar results have been obtained reweight- 
ing between MRST2002 and ZEUS-S PDF's. 



This technique has been tested using HERWIG (for inclusive W production) and ALP- 
GEN interfaced to HERWIG (for W+jets production) as Monte Carlo generators and with 
various recent PDF sets. Similar results have been obtained with these two MC gener- 
ators and with different PDF sets, as discussed below. Fig. 4.23 shows the accuracy of 
this technique using HERWIG: the bias over the all y range is of the order of 0.5% or 
less and there is no evidence of y dependence. Comparing the bottom plots on the right 
and left hand sides of fig. 4.23 we see that the PDF reweighting technique corrects for the 
difference in normalisation between PDF\ and PDF2 and corrects for the y modulation. 

This technique can be used to estimate the full PDF uncertainty, starting from one 
sample of MC generated events, for distributions that are determined by the MC hard 
process. 
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NNLO Evolution of the Pdf 's and 
their Errors: Benchmarks and 
Predictions for Drell-Yan 

Alessandro Cafarella, Claudia Coriand, Marco Guzzi 

Abstract 

We quantify the impact of the next-to-next-to-leading order evolution on the Drell-Yan 
total cross section and on the corresponding rapidity distributions of the lepton pair and 
compute the corresponding errors. We base our analysis on Candia, a program that solves 
the DGLAP equations using the method of the x-space iterates. 

5.27 Introduction 

In the search for new physics at the LHC we need high precision in the determinations 
of the QCD background, possibly at next-to-next-to-leading order (NNLO) in the strong 
coupling constant a^. While it is expected that only a few processes will be computed in the 
near future at this order of accuracy in QCD, for Drell-Yan lepton pair production, some of 
these corrections - for instance those involving the invariant mass distributions - have been 
available for some time [1 j . The study of this process will be essential both in the search 
of extra neutral interactions and for partonometry, where the impact of the perturbative 
resummation |2] can be studied in detail given the large number of events expected at 
the LHC. More recently, following the computation of the invariant mass distribution 
da/dQ, where is the invariant mass of the lepton pair, also the rapidity distributions 
da/{dYdQ), implemented in VRAP [3] have been computed. A fully exclusive numerical 
computation has also been presented Next-to-leading order (NLO) analysis of the 
forward-backward asymmetries ApB on the Z resonance and NNLO charge asymmetries 
have also been determined [H] . The computation of the hard scatterings for some of these 
processes has been performed much before that the analytical computation of the NNLO 
evolution kernels needed for a consistent extraction of NNLO parton distributions (pdf 's) 
were available. Following the computation of the kernels [5], some benchmarks for the 
NNLO evolution have been presented [6], followed by a later update [7]. 

Testing the benchmarks by using independent approaches that solve the equations is 
not only a demanding numerical problem, but involves subtle issues concerning the types 
of solutions that are implemented in a given numerical algorithm. Specifically, we have 
shown in [HI |9] that the selection of a given ansatz - either in Mellin space or in x-space - in 
the solution of the DGLAP involves a specific arrangement of the logarithmic expansion 
that solves iteratively the equations. For instance, ansatze for the exact solutions perform 
automatically a resummation of the contributions identified by the simplest logarithmic 
ansatz (also called "truncated solutions"). These involve logarithms of the ratio of two 
couplings log(as((5^)/as((5g)) at two different scales Q and Qq. Exact solutions, instead. 
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replace these logarithmic expansions with more complicated functions of a^. Details can 
be found in m El HO] . 

These expansions apply generically both to forward and non- forward twist-2 operators, 
and converge to the solutions of the evolution equations with very high precision. Concern- 
ing the structure of the evolution codes, these are usually based either on the numerical 
Mellin inversion (using an ansatz in moment space) or on "brute force methods". In this 
second case the numerical solution is built by a discretization of the equations, reduced to 
a stable finite-difference scheme. The theoretical indetermination coming from the various 
approaches, as we have pointed out in [9], has to do with the the selected accuracy of the 
solution. 



5.28 The choice of the solution and the theoretical indeter- 
mination 

We have re-analized the issue of the initial state dependence of the predictions for the 
da/dQ and da / {dydQ) cross sections in |P, TU], focusing our attention near the Zq peak. 

Being the predictions coming from the inclusions of the NNLO corrections quite small, 
it is natural to worry whether the errors coming from the statistical fits in the determina- 
tion of the pdf's to the experimental data and the indeterminations due to the treatment of 
the evolution can be compared so that the size of the errors claimed in the hadronic cross 
sections are quantified consistently. We point out that the specific choice of the solution 
of the evolution equation brings in an intrinsic indetermination which is comparable with 
the error coming from the fits to the pdf's. This indetermination is of theoretical origin. 

In our analysis we will omit most of the theoretical details that can be found in |10j 
and come to a discussion of the numerical results, since these can be of interest for actual 
experimental searches. We also present some benchmarks for the NNLO evolution that can 
be useful for a consistent comparison with other codes. An extensive numerical analysis 
of the predictions on the Z resonance that illustrates the points summarized here can be 
found in [9 . 

As we are going to see, the variations induced by the choice of the solution of the 
DGLAP induce variations on the cross section of the order of 1% or so at NNLO, and 
clearly affect also the NNLO i^-factor for the total cross section. We have used for this 
study Candia [8], which implements the truncated and exact solutions of the DGLAP 
built without the numerical Mellin inversions. On the countrary of other programs that 
need initial conditions of a specific functional form, Candia allows any initial conditions 
to be studied, being based on x-space algorithms. An extended version of Candia, called 
CandiAdy allows to study the invariant mass distribution to NNLO. An interface with 
VRAP is also under development. 



5.29 Benchmarks and Numerical Results 

Now, we come to illustrate a comparison between the NNLO evolution performed by 
Candia and Pegasus [Hi, using the initial conditions of the Les Houches model. This 
model works as a benchmark and allows to compare various evolution codes. We show the 
behaviour of the gluon distribution and the valence up quark distribution (non-singlet). 
We observe that the differences in the singlet case are around 0.4 — 0.3% or less, while for 
the valence up quark they are around 3% at x = 10~^ and then decrease to 0.1 — 0.02% at 



X = 0.1 (see Table 5.4). We have denoted with x5f{x) the relative differences normalized 
to the Pegasus determination, i.e. x5f{x) = (a;/(x)pEGAsus - 2;/(a;)cANDiA)/a^/(a;)pEGAsus- 
We have set the factorization and renormalization scales to coincide and equal io Q = 100 
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Figure 5.24: Cross Sections in the region of the peak of the Z boson at LO, NLO, and 
NNLO obtained using the luminosities evolved respectively by Alekhin and MRST 
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Candia vs Pegasus PDFs at NNLO, Les Houches input, VFN scheme, Q = /xf = Mfl = 100 GeV 
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Table 5.4: NNLO pdf's determined with Candia and Pegasus using the Les Houches 
model. 



dffATjvLo/dQ [pb/GeV]. Candia vs Pegasus with Les Houches input. 
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"NNLO 
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Table 5.5: NNLO cross sections in the two evolution methods. 



GeV. In Candia we have used the "asymptotic solutions", which are similar to those 
obtained by PEGASUS in one of its 3 running modes. The two algorithms and their 
implementations compare very well in the x-region relevant for the LHC (around 0.1% at 
NNLO). Using the same benchmark we have calculated the NNLO cross sections using 
the two evolutions (see Tab. 5.5), and we observe that the differences between the two 
methods in the kinematical region that we are considering are around 0.2 — 0.3%. 

Coming t o a d escription of the NNLO cross section obtained from a realistic model, we 
show in Tab. 
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results for the invariant mass distributions at v "S" = 14 TeV according to 
Alekhin's model pL2]. Shown are also the errors which have been computed as discussed 
in [9]. A similar analysis has been performed using the MRST set [13^ (see Tab. 5.8), for 
which the errors can be obtained, at this time, only at NLO. 

The differences between our prediction and the MRST result for the total cross sections 
are around 1 per cent or below at LO, vary from 0.02% to 0.3% at NLO and are 2.6% and 
below at NNLO (see Tab. 5.7). In this case the maximum difference has been found for 
Q = 50 GeV. These differences, clearly, affect the values of the -fC-factors, as we are going 
to discuss below, which in our evolution are larger compared to those of MRST. 

We have summarized in Fig. 5.25 four plots of the behavior of the 3 IT- factors K = 
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Figure 5.25: Various if- factors obtained with the evolution performed by Candia and 
MRST. 



dcr/dQ in [pb/GeV] for Alekhin with Q'^ = fxj, = ^/S=14. TeV 


Q [GcV] 




<^NLO 


O'NNLO 


50 


6.22 ± 0.27 


7.48 ± 0.24 


7 A3 ± 0.21 


60.04 


3.72 ± 0.15 


4.50 ± 0.13 


4.49 ± 0.12 


70.1 


3.30 ± 0.12 


4.03 ± 0.11 


4.05 ± 0.10 


80.1 


6.65 ± 0.24 


8.20 ± 0.24 


8.19 ± 0.23 


90.19 


253 ± 8 


313 ± 9 


309 ± 8 


91.19 


415 ± 14 


514 ± 15 


506 ± 15 


120.07 


0.80 ± 0.02 


0.99 ± 0.02 


0.96 ± 0.03 


146.1 


0.225 ± 0.006 


0.277 ± 0.007 


0.269 ± 0.007 


172.1 


0.097 ± 0.002 


0.119 ± 0.003 


0.117 ± 0.003 


200 


0.047 ± 0.001 


0.058 ± 0.001 


0.058 ± 0.001 



Table 5.6: Cross sections derived from the best fits for the 3 orders with their errors for 
the set by Alekhin. 
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dcTATiVLo/dQ [pb/GeV]. Candia vs MRST evolution with MRST input, ixl = 1.25 GeV^ 
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Table 5.7: NNLO cross section for Drell-Yan obtained by Candia using the MRST input 
and the evolved MRST pdf 's 



da/dQ in [pb/GeV] for MRST with Q2 = = ^/S = U TeV 


Q [GeV] 




50 


6.77 ± 0.19 


60.04 


4.13 ± 0.10 


70.1 


3.79 ± 0.08 


80.1 


7.90 ± 0.14 


90.19 


305 ± 5 


91.19 


499 ± 8 


120.1 


0.952 ± 0.014 


146.1 


0.264 ± 0.003 


172.1 


0.113 ± 0.001 


200 


0.0556 ± 0.0007 



Table 5.8: Cross sections derived from the best fits at NLO with the errors for the MRST 
set. 
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(yNNLol^NLO, Ki = aNLo/o'LO and K2 = ctnnlo/c^lo obtained using Candia and the 
MRST evolution. These are shown as a function of Q, and evaluated at the center of mass 
energy of ^/S = 14 TeV. The dependence of the results on the evolution is significant. 



In fact, from Fig. 5.25 it is evident that while the shapes of the plots of the factors 
are similar, there are variations of the order 2%, in the results using the two different 
evolutions. Both in the evolution performed with Candia and in the MRST evolution 
we use the same MRST input, choosing the initial scale /Uq = 1.25 GeV^, and the same 
treatment of the heavy flavors. On the Z resonance we get 

K{Mz) = {gnnlo ® KmsT)/{<^NLO ® ^mwst) = 0-97 
K{Mz) = {aNNLO ® ^^a^ndPa)/(^^lo ® f ^i^DiA) = 0.95 

K{Mz) = {dNNLO ® ^Allthin)l{^NLO ® ^Atekhin) = 0-98 (5.5) 

which corresponds to a reduction by 2.7% of the NNLO cross section compared to the NLO 
result, (MRST evolution) and larger for the Candia evolution, 4.4%, while for Alekhin is 
1.5%. Prom the analysis of the errors on the pdf 's to NNLO, for instance for the Alekhin's 
set, the differences among these determinations are still compatible, being the variations 
on the iC- factors of the order of 4%. 



5.30 The rapidity distributions 

In this case the QCD cross section is given by 



dY 

ab 



E / / dxldx2f^^{x,,Q^/^^l,^iy^il)f^%x2,Qy^^l,^^y^il)x 



dcT^ 

^(xi,X2,gVM^,4/Ml). (5.6) 

Notice that the evolution implemented in Candia allows to analyze the renormaliza- 
tion/factorization scale dependence also in the evolution, which is not present in the 
MRST parameterizations. We have made explicit this dependences in ( |5.6| ). 

We have set the scales to be equal, fip = fJ-R and varied //i? in the interval 1/2Q < fip 1^ 
2Q, obtaining results which differ from those obtained in [3] by 2% due to the different 
implementation of the evolution. Using Candia and as initial condition the MRST grid 
input with /ig = 1.25 GeV^ the NNLO band and the NLO one are resolved separately. 
We have also found that with the inclusion of the A*^/a*|' effects in the pdf's evolution, 
the dependence on is quite sizeable at NLO, but is reduced at NNLO |9]. We show in 



Fig. (5.26) the plots of the variations of the rapidity distributions at the three orders and 
the corresponding errors on the pdf's for Alekhin's model and for MRST for Q = Mz- In 
both cases the reduction of the variation of the cross sections as we move toward higher 
orders is quite evident. We report also the errors on these distributions obtained in both 
models, which get systematically smaller as the accuracy of the calculation increases. 



5.31 Conclusions 

In the search for extra neutral currents precise theoretical determinations of lepton pair 
production via the Drell-Yan mechanism are going to play a very relevant role (see Ref. 

[laiEiiIIiiITiiigiig^). 

In our determination, the change in the value of the cross section from NLO to NNLO 
is around 4% on the Z peak, while the MRST and the Alekhin determinations are 2.6% 
and about 1.5% respectively. While these variations appear to be more modest compared 
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Alekhin pdf s error, LO ' — ' — 

LO 2 > ^ip/ Q > 1/2 

NLO pdf s error ^--x - 

NLO 2 > Hp/ Q > 1/2 

NNLO pdf s error ^ - - - 

NNLO 2 > |ip/ Q > 1/2 




J I i I i H i i H I n I H I i I J 
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Y 

(a) Alekhin's model 



MRST pdf s enor, NLO 
NLO 2 > |lp/ Q > 1/2 




Y 

(b) MRST model 



Figure 5.26: Plot of the rapidity distributions at LO, NLO and NNLO for Alekhin's model 
and MRST. Shown are also the bands due to the variation of the fip scale, and the errors 
on the cross sections at the corresponding orders. 
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to the analogous ones at a lower order (which are of the order of 20% or so), they are 
nevertheless important for the discovery of extra neutral currents at large invariant mass 
of the lepton pair in DY, given the fast falling cross section at those large values. The 
errors on the pdf 's induce percentile variations of the cross section as we move from NLO 
to NNLO of the order of 4% around the best-fit result, reducing the NNLO cross section 
compared to the NLO prediction and rendering these results compatible. 
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Combination of QCD and 
electroweak corrections to 
Drell-Yan processes 

G. Balossini, CM. Carloni Calame, G. Montagna, M. Moretti, O. Nicrosini, F. Piccinini, 
M. Treccani and A. Vicini 

Introduction 

Precision measurements of electroweak (EW) gauge boson production and properties will 
be a crucial goal of the physics program of proton-proton collisions at the LHC. W and Z 
bosons will be produced copiously and careful measurements of their observables will be 
important in testing the Standard Model (SM) and uncovering signs of new physics [1] |2] . 

Thanks to the high luminosity achievable at the LHC, the systematic errors will play a 
dominant role in determining the accuracy of the measurements, implying, in particular, 
that the theoretical predictions will have to be of the highest standard as possible. For 
Drell-Yan (D-Y) processes, this amounts to make available calculations of W and Z pro- 
duction processes including simultaneously higher-order corrections coming from the EW 
and QCD sector of the SM. Actually, in spite of a detailed knowledge of EW and QCD 
corrections separately, the combination of their effects have been addressed only recently 
[HI m IS] and need to be deeply scrutinized in view of the anticipated experimental accuracy. 

In this contribution, after a review of existing calculations and codes, we present 
the results of a study aiming at combining EW and QCD radiative corrections to D-Y 
processes consistently. We do not include in our analysis uncertainties due to factoriza- 
tion/renormalization scale variations, as well as uncertainties in the Parton Distribution 
Functions arising from diverse experimental and theoretical sources, which are left to a 
future publication. Some results already available in this direction can be found in [6]. 

Status of theoretical predictions and codes 

Concerning QCD calculations and tools, the present situation reveals quite a rich structure, 
that includes next-to-leading-order (NLO) and next-to-next-to-leading-order (NNLO) cor- 
rections to W/Z total production rate [HIH], NLO calculations for W, Z+1,2 jets signatures 
[SI [JUl (available in the codes DYRAD and MCFM), resummation of leading and next-to- 
leading logarithms due to soft gluon radiation [T2] (implemented in the Monte Carlo 
ResBos), NLO corrections merged with QCD Parton Shower (PS) evolution (in the event 
generators MC@NLO and POWHEG [H]), NNLO corrections to W/Z production in 
fully differential form [ 15| [T6] (available in the Monte Carlo program FEWZ), as well as 
leading-order multi-parton matrix elements generators matched with vetoed PS, such as, 
for instance, ALPGEN [H], MADEVENT [H], HELAC [H] and SHERPA ^20]. 
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As far as complete 0(a) EW corrections to D-Y processes are concerned, they have 
been computed independently by various authors in [211 EZl [23l EH l2S] for W production 
and in |26| 1271 EHl [29] for Z production. Electroweak tools implementing exact NLO 
corrections to W production are DK [21], WGRAD2 [22], SANG [S] and HORACE [25], 
while ZGRAD2 [26], HORACE [28] and SANG [29] include the full set of 0{a) EW 
corrections to Z production. The predictions of a subset of such calculations have been 
compared, at the level of same input parameters and cuts, in the proceedings of the 
Les Houches 2005 |30j and TEV4LHG [31] workshops for W production, finding a very 
satisfactory agreement between the various, independent calculations. A first set of tuned 
comparisons for the Z production process has been recently performed and is available 
in [M]. 

From the calculations above, it turns out that NLO EW corrections are dominated, in 
the resonant region, by final-state QED radiation containing large collinear logarithms of 
the form log(s/m^), where s is the squared partonic centre-of-mass (cm.) energy and mi 
is the lepton mass. Since these corrections amount to several per cents around the jacobian 
peak of the W transverse mass and lepton transverse momentum distributions and cause a 
significant shift (of the order of 100-200 MeV) in the extraction of the W mass Myy at the 
Tevatron, the contribution of higher-order corrections due to multiple photon radiation 
from the final-state leptons must be taken into account in the theoretical predictions, 
in view of the expected precision (at the level of 15-20 MeV) in the M\y measurement 
at the LHG. The contribution due to multiple photon radiation has been computed, by 
means of a QED PS approach, in [33| for W production and in [33] for Z production, and 
implemented in the event generator HORACE. Higher-order QED contributions to W 
production have been calculated independently in [35] using the YFS exponentiation, and 
are available in the generator WINHAC. They have been also computed in the collinear 
approximation, within the structure functions approach, in [36 . 

A further important phenomenological feature of EW corrections is that, in the region 
important for new physics searches (i.e. where the W transverse mass is much larger 
than the W mass or the invariant mass of the final state leptons is much larger than 
the Z mass), the NLO EW effects become large (of the order of 20-30%) and negative, 
due to the appearance of EW Sudakov logarithms oc -(a/vr) log^(s/M^), V = W, Z 
[2T| l22l l25l [26l [271 [28] . Furthermore, in this region, weak boson emission processes (e.g. 
pp — > e^Vey + that contribute at the same order in perturbation theory, can partially 
cancel the large Sudakov corrections, when the weak boson V decays into unobserved vi' 
or jet pairs, as recently shown in |37j . 

Theoretical approach 

A first strategy for the combination of EW and QCD corrections consists in the following 
formula 



da 
dO 



QCD&EW 



da_]^ 
dOj 



+ 



MCONLO 





da 




da 




{ 


dO 


EW 


dO 


Born } 



(6.7) 



where dcr/dOMCQNLO stands for the prediction of the observable da/dO as obtained by 
means of MC@NLO, da/dO^^^ is the HORACE prediction for the EW corrections to the 
da/dO observable, and da / dO^^^^ is the lowest-order result for the observable of interest. 
The label HERWIG PS in the second term in r.h.s. of eq. (6.7) means that EW corrections 



are convoluted with QCD PS evolution through the HERWIG event generator, in order 
to (approximately) include mixed 0{aas) corrections and to obtain a more realistic de- 
scription of the observables under study. However, it is worth noting that the convolution 
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of NLO EW corrections with QCD PS implies that the contributions of the order of aag 
are not rehable when hard non-cohinear QCD radiation turns out to be relevant, e.g. for 
the lepton and vector boson transverse momentum distributions in the absence of severe 
cuts able to exclude resonant W/Z production. In this full 0{aas) calculation 

would be needed for a sound evaluation of mixed EW and QCD corrections. Full 0{a) 
EW corrections to the exclusive process pp ^ W + j (where j stands for jet) have been 
recently computed, in the approximation of real W bosons, in |3Sjl39j, while one-loop weak 
corrections to Z hadro-production have been computed, for on-shell Z bosons, in [ID]. It 



is also worth stressing that in eq. (6.7) the infrared part of QCD corrections is factorized. 



whereas the infrared-safe matrix element residue is included in an additive form. It is 
otherwise possible to implement a fully factorized combination (valid for infra-red safe 
observables) as follows: 



da 
dO 



QCD0EW 



[da/dO] 



1 + 



^ { dOEW } 



MC@NLO 



[^"^/^^Iherwig PS 



[da/dO] 



Born 



(6.e 



HERWIG PS 



where the ingredients are the same as in eq. (6.7) but also the QCD matrix element 



residue in now factorized. Eqs. (6.7) and (6.8) have the very same 0{a) and 0{a 



content, differing by terms of the order of aos- Their relative difference has been checked 
to be of the order of a few per cent in the resonance region around the W/Z mass, and 
can be taken as an estimate of the uncertainty of QCD and EW combination. 



Numerical results: W and Z production 

In order to assess the phenomenological relevance of the combination of QCD and EW 
corrections, we study, for definiteness, the charged-current process pp — > fi^ + X 
at the LHC, imposing the following selection criteria 

a. > 25 GeV, > 25 GeV, \r]^\ < 2.5, 

b. the cuts as above > 1 TeV, (6.9) 

where and t]^ are the transverse momentum and the pseudorapidity of the muon, 
is the missing transverse energy, which we identify with the transverse momentum of the 
neutrino, as typically done in several phenomenological studies. For set up b., a severe 
cut on the W transverse mass is superimposed to the cuts of set up a., in order to 
isolate the region of the high tail of MjT, which is interesting for new physics searches. 
We also consider the neutral-current reaction pp ^ j,Z ^ e'^e' +X, selecting the events 
according to the cuts 

pf > 25 GeV, |7/^*| < 2.5, M^+^- > 200 GeV. (6.10) 

The granularity of the detectors and the size of the electromagnetic showers in the calorime- 
ter make it difficult to discriminate between electrons and photons with a small opening 
angle. We adopt the following procedure to select the event: we recombine the four- 
momentum vectors of the electron and photon into an effective electron four-momentum 
vector if, defining 

AR{e, 7) = V^vie, 7)' + A0(e, 7)2, (6.11) 

AR(e, 7) < 0.1 (with Aij, the distances of electrons and photons along the longitudinal 
and azimuthal directions). We do not recombine electrons and photons if rj^ > 2.5 (with 
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rj.y the photon pseudo-rapidity). We apply the event selection cuts as in Eq. (6.10) only 
after the recombination procedure. 

The parton distribution function (PDF) set MRST2004QED tilj has been used to 
describe the proton partonic content. The QCD factorization/renormalization scale and 
the analogous QED scale (present in the PDF set MRST2004QED) are chosen to be 
equal, as usually done in the literature [21] |22l |25l EH |28], and fixed at = fip = 

\J iPi^)^ ~^ "^M^^M charged-current case), where M^y^ is the //f^ invariant mass, 

and at = fip = \J {p±Y + ^e+e- neutral-current case), where Mg+g- is the 

invariant mass of the lepton pair. 

In order to avoid systematics theoretical effects, all the generators used in our study 
have been properly tuned at the level of input parameters, PDF set and scale to give the 
same LO/NLO results. The tuning procedure validates the interpretation of the various 
relative effects as due to the radiative corrections and not to a mismatch in the setups of 
the codes under consideration. 



LO HORACE + HERWIG PS ■ 
HCSNLO - 
HCaNLO+HORACE - 
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Figure 6.27: Upper panel: predictions of MC@NLO, MC@NLO-hHORACE and leading- 
order HORACE+HERWIG PS for the (left) and (right) distributions at the LHC, 
according to the cuts of set up a. of Eq. (6.9). Lower panel: relative effect of QCD and 



EW corrections, and their sum, for the corresponding observables in the upper panel. 



A sample of our numerical results is shown in Fig. 6.27 for the W transverse mass 
Mj^ and muon transverse momentum distributions according to set up a. of Eq. (6.9), 
and in Fig. 6.28| for the same distributions according to set up b. In Fig. 6.27 and Fig. 



6.28, the upper panels show the predictions of the generators MC@NLO and MC@NLO 



-|- HORACE interfaced to HERWIG PS (according to eq. (6.7)), in comparison with the 
leading-order result by HORACE convoluted with HERWIG shower evolution. The lower 
panels illustrate the relative effects of the matrix element residue of NLO QCD and of full 
EW corrections, as well as their sum, that can be obtained by appropriate combinations 
of the results shown in the upper panels. More precisely, the percentage corrections shown 
have been defined as 6 = ((Tnlo - o-Born+HERWiGPs) /o-Bom+HERWiG PS , where ctnlo stands 
for the predictions of the generators including exact NLO corrections matched with QCD 
PS. 

From Fig. 6.27 it can be seen that the QCD corrections are positive around the W 

and tend to compensate the negative effect due to EW 



jacobian peak, of about 10-zu7o 
corrections. Therefore, their interplay is crucial for a precise Mw extraction at the LHC 
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and their combined contribution can not be accounted for in terms of a pure QCD PS 
approach, as it can be inferred from the comparison of the predictions of MC@NLO versus 
the leading-order result by HORACE convoluted with HERWIG PS. It is also worth noting 
that the convolution of NLO corrections with the QCD PS broadens the sharply peaked 
shape of the fixed-order NLO QCD and EW effects. 

The interplay between QCD and EW corrections to W production in the region in- 
teresting for new physics searches, i.e. in the high tail of and distributions, is 
shown in Fig. 6.28 For both and pj_, the QCD corrections are positive and largely 
cancel the negative EW Sudakov logarithms. Therefore, a precise normalization of the 
SM background to new physics searches necessarily requires the simultaneous control of 
QCD and EW corrections. 
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Figure 6.28: The same as Fig. 6.27 according to the cuts of set up b. of Eq. (6.9). 



Figure 6.29 shows the combination of QCD and EW corrections for the di-lepton 
invariant mass in the neutral-current D-Y process — > 7, Z ^ e~^e^ + X, according to 
the cuts of Eq. (6.10 1 |l2]. The QCD corrections are quite flat and positive with a value of 



about 15% over the mass range 200-1500 GeV. The EW corrections are negative and vary 
from about —5% to —10% and thus partially cancel the QCD contribution. Therefore, as 
for the charged-current channel, the search for new physics in di-lepton final states needs 
a careful combination of EW and QCD effects. 



Conclusions 



During the last few years, there has been a big effort towards high-precision predictions for 
D-Y-like processes, addressing the calculation of higher-order QCD and EW corrections. 
Correspondingly, precision computational tools have been developed to keep under control 
theoretical systematics in view of the future measurements at the LHC. 

We presented some original results about the combination of EW and QCD corrections 
to a sample of observables of W and Z production processes at the LHC. Our investigation 
shows that a high-precision knowledge of QCD and a careful combination of EW and strong 
contributions is mandatory in view of the anticipated experimental accuracy. We plan, 
however, to perform a more complete and detailed phenomenological study, including the 
predictions of other QCD generators and considering further observables of interest for 
the many facets of the W/Z physics program at the LHC. 
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Figure 6.29: Left panel: the di-electron invariant mass distribution according to 
the leading-order and NLO EW predictions of HORACE, of MC@NLO and of 
MC@NLO+HORACE at the LHC, using the cuts of Eq. ( |6.10[ ). Right panel: relative 
effect of QCD and EW corrections, and their combination. 
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Subtraction at NNLO and Higgs 
boson production at hadron 
colliders 

M. Grazzini 

7.32 Introduction 

The dynamics of hard scattering processes involving hadrons is nowadays remarkably well 
described by perturbative QCD predictions. Thanks to asymptotic freedom, the cross 
section for sufficiently inclusive reactions can be computed as a series expansion in the 
QCD coupling as. Until few years ago, the standard for such calculations was next-to- 
leading order (NLO) accuracy. Next-to-next-to-leading order (NNLO) results were known 
only for few highly-inclusive reactions (see e.g. Refs. [U El E]). 

The extension from NLO to NNLO accuracy is important to improve QCD predictions 
and to better assess their uncertainties. In particular, this extension is essential in two 
cases: in those processes whose NLO corrections are comparable to the leading order (LO) 
contribution; in those 'benchmark' processes that are measured with high experimental 
precision. Such a task, however, implies finding methods and techniques to cancel the 
infrared (IR) divergences that appear at intermediate steps of the calculations. 

Recently, a new general method [3], based on sector decomposition fSj, has been pro- 
posed and applied to the NNLO QCD calculations of e+e~ — > 2 jets |i6j, Higgs [7^ and 
vector [H] boson production in hadron collisions, and to the NNLO QED calculation of 
the electron energy spectrum in muon decay The calculations of Refs. [71 [8] allow us 
to compute the corresponding cross sections with arbitrary cuts on the momenta of the 
partons produced in the final state. 

The traditional approach to perform NLO computations is based on the introduction 
of auxiliary cross sections that are obtained by approximating the QCD scattering am- 
plitudes in the relevant IR (soft and coUinear) limits. This strategy led to the proposal 
of the subtraction [10\ and slicing methods. Exploiting the universality properties 
of soft and collinear emission, these methods were later developed in the form of general 
algorithms |12| [T3l ITS] . that make possible to perform NLO calculations in a (relatively) 
straightforward manner, once the corresponding QCD amplitudes are available. In recent 
years, several research groups have been working on general NNLO extensions of the sub- 
traction method [151 [161 [T71 [181 [19] . Results have been obtained in some specific processes: 
e+e" 2 jets ^ [21] and, more recently, e+e" 3 jets [221 [23]. 

In Ref. [2^ we proposed an extension of the subtraction method to NNLO for a specific, 
though important, class of processes: the production of colourless high-mass systems in 
hadron collisions. We presented a formulation of the subtraction method for this class 
of processes, and we applied it to the NNLO calculation of Higgs boson production via 
the gluon fusion subprocess gg H. The calculation has now been implemented in 
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the numerical program HNNLO, which includes all the relevant decay modes of the Higgs 
boson for this production subprocess, namely, H ^ j'y [23], H — > WW Ivlv and 
H^ZZ^Al[2b\. 



This contribution is organized as follows. In Sect. 7.33 we discuss the version of the 



subtraction formalism we use. In Sect. 7.34 we present a selection of numerical results 



that can be obtained by our program. In Sect. 7.35 we summarize our results. 

7.33 The method 

We consider the inclusive hard-scattering reaction 

hi + h2^ F{Q)+X, (7.12) 

where the collision of the two hadrons hi and /12 produces the triggered final state F. 
The final state F consists of one or more colourless particles (leptons, photons, vector 
bosons, Higgs bosons, . . . ) with momenta qi and total invariant mass Q . Note that, since 
F is colourless, the LO partonic subprocess is either qq annihilation, as in the case of the 
Drell-Yan process, or gg fusion, as in the case of Higgs boson production. 

At NLO, two kinds of corrections contribute: i) real corrections, where one parton 
recoils against F; ii) one-loop virtual coriections to the LO subprocess. Both contributions 
are separately IR divergent, but the divergences cancel in the sum. At NNLO, three kinds 
of corrections must be considered: i) double real contributions, where two partons recoil 
against F; ii) real-virtual corrections, where one parton recoils against F at one-loop 
order; iii) two-loop virtual corrections to the LO subprocess. The three contributions are 
still separately divergent, and the calculation has to be organized so as to explicitly achieve 
the cancellation of the IR divergences. 

Our method is based on a generalization of the procedure used in the specific NNLO 
calculation of Ref. |26]. We first note that, at LO, the transverse momentum of the 
triggered final state F is exactly zero. As a consequence, as long as / 0, the (N)NLO 
contributions are actually given by the (N)LO contributions to the triggered final state 
F -|- jet(s). Thus, we can write the cross section as 

d'^(N)NLo\<lT7^0 = ^^{JV)LO • ("^-13) 

This means that, when qT ^ 0, the IR divergences in our NNLO calculation are those in 
dcT^^Q^^: they can be handled and cancelled by using available NLO formulations of the 
subtraction method. The only remaining singularities of NNLO type are associated to the 
limit qT — > 0, and we treat them by an additional subtraction. Our key point is that the 
singular behaviour of da^j^-^^^ when — > is well known: it appears in the resummation 
program |27| |28| [29] of logarithmically-enhanced contributions to transverse-momentum 
distributions. Then, to perform the additional subtraction, we follow the formalism used 
in Ref. |30| |3T] to combine resummed and fixed-order calculations. 

We use a shorthand notation that mimics the notation of Ref. [3D]. We define the 
subtraction counterternj^ 

da^^ = dalo » ^^{Qt/Q) d^T- (7.14) 

The function Ti^{qT/Q) embodies the singular behaviour of da^'^^^'^^ when qx — > 0. In 
this limit it can be expressed as follows in terms of QT-independent coefficients Yi^i"^'^^)-, 

®The symbol ® understands convolutions over momentum fractions and sum over flavour indeces of the 
partons. 
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The extension of Eq. (7.13) to include the contribution at qx = is finally: 



da 



{N)NLO 



n 



{N)NLO 



da 



LO 



+ 



da 



F+jets 
(N)LO 



"'^{N)LO 



(7.16) 



Comparing with the right-hand side of Eq. (7.13), we have subtracted the truncation of 
Eq. (7.14) at (N)LO and added a contribution at qx = needed to obtain the correct 
total cross section. The coefficient "H^^^^^q does not depend on qx and is obtained by 
the (N)NLO truncation of the perturbative function 



n 



TT 



asy 

TT / 



(7.17) 



The counterterm of Eq. (7.14) regularizes the singularity of da^~^^^^^ when qx 0: the 
term in the square bracket on the right-hand side of Eq. (7.16) is thus IR finite (or, better, 
integrable over qx). Note that, at NNLO, da^^^j^Q acts as a counterterm for the sum of 

the two contributions to da^~^^'^^^: the double real plus real-virtual contributions. We also 
note that the counterterm function Tj^{qx/Q) can be defined in different ways: the only 
property we require is that in the small qx limit it must take the form given in Eq. (7.15), 
so as to match the singular behaviour of da^^'^'^^^. Note that the perturbative coefficients 
Y^F{n\k) g^j,g universal: more precisely, the NNLO coefficients S^^^'^^ and E^^^'^^ have a 
non-universal contribution that, nonetheless, is proportional to the NLO coefficient Tl^^^\ 
The above coefficients only depend on the type of partons (quarks or gluon) involved in the 
LO partonic subprocess {qq annihilation or gg fusion). We finally note that the simplicity 
of the LO subprocess is such that final-state partons actually appear only in the term 
da^^^^^^ on the right-hand side of Eq. (7.16). Therefore, arbitrary IR-safe cuts on the jets 
at (N)NLO can effectively be accounted for through a (N)LO computation. Owing to this 
feature, our NNLO extension of the subtraction formalism is observable-independent. 

At NLO (NNLO), the physical information of the one-loop (two-loop) mrtita/ correction 
to t he LO subprocess is contained in the coefficients 7i^^^ {'H^'^^). Once an explicit form of 
Eq. (7.14) is chosen, the hard coefficients TC^^""^ are unique ly ide ntified (a different choice 
would correspond to different Ti^^"^^). According to Eq. (7.16), the NLO calculation of 
da^ requires the knowledge of TL^^^"^ and the LO calculation of da^'^^'^^^. The general 
(process-independent) form of the coefficient TL^^^^ is basically known: the precise relation 
between 7i^^^^ and the IR finite part of the one-loop correction to a generic LO subprocess 
is explicitly derived in Ref. [32;. 

At NNLO, the coefficient Ti^^'^^ is also needed, together with the NLO calculation of 
da^^^"^^^. The coefficients Ti^^"^^ for Higgs boson production in the large- M^p limit have 
been computed [33]. Since the NLO corrections to gg — > i?-|-jet(s) are available Pl] in the 



same limit, we are able to apply Eq. (7.16) at NNLO. We have encoded our computation 
in a parton level Monte Carlo program, in which we can implement arbitrary IR-safe cuts 
on the final state. 



7.34 Results 

In the following we present numerical results for Higgs boson production at the LHC. 
We use the MRST2004 parton distributions [35] , with densities and as evaluated at each 
corresponding order (i.e., we use {n -\- l)-loop as at N^LO, with n = 0, 1,2). The renor- 
malization and factorization scales are fixed to the value /i/j = fip = Mh, where Mh is 
the mass of the Higgs boson. 

7.34.1 H ^-f-f 

We consider the Higgs boson decay in the H ^ ^'j channel and follow Ref. [36j to apply 
cuts on the photons. For each event, we classify the photon transverse momenta according 
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to their minimum and maximum value, PTmin and PTmax- The photons are required to 
be in the central rapidity region, \rj\ < 2.5, with pTmm > 35 GeV and pTmax > 40 GeV. 
We also require the photons to be isolated: the hadronic (partonic) transverse energy in 
a cone of radius R = 0.3 along the photon direction has to be smaller than 6 GeV. When 
Mh = 125 GeV, by applying these cuts the impact of the NNLO corrections on the NLO 
total cross section is reduced from 19% to 11%. 
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Figure 7.30: Distributions in pTmin o-nd PTmax for the diphoton signal at the LHC. The 
cross section is divided by the branching ratio in two photons. 



In Fig. 7.30 we plot the distributions in pTmm and pTmax for the gg ^ H ^ 
signal. We note that the shape of these distributions sizeably differs when going from 
LO to NLO and to NNLO. The origin of these perturbative instabilities is well known 
[57] . Since the LO spectra are kinematically bounded by px < Mh/'^^-, each higher-order 
perturbative contribution produces (integrable) logarithmic singularities in the vicinity of 
that boundary. More detailed studies are necessary to assess the theoretical uncertainties 
of these fixed-order results and the relevance of all-order resummed calculations. 



In Fig. 7.31 we consider the (normalized) distribution in the variable cos 9* , where 9* is 
the polar angle of one of the photons in the rest frame of the Higgs boson. At small values 
of cos^* the distribution is quite stable with respect to higher order QCD corrections. We 
also note that the LO distribution vanishes beyond the value cos^J^j^^ < 1- The upper 
bound cos ^j^ax is due to the fact that the photons are required to have a minimum px of 



35 GeV. As in the case of Fig. |7.30 in the vicinity of this LO kinematical boundary there 



is an instability of the perturbative results beyond LO. 
7.34.2 H WW lulu 

We now consider the production of a Higgs boson with mass = 165 GeV in the decay 
mode H — > WW — > lulu [25] • We apply a set of selection cuts taken from the study of 
Ref. |38j . The charged leptons are classified according to their minimum and maximum 
Pt- The PTmin should be larger than 25 GeV, and PTmax should be between 35 and 50 
GeV. The charged lepton rapidity should fuUfil \rj\ < 2. The missing px of the event is 
required to be larger than 20 GeV and the invariant mass of the charged leptons is smaller 
than 35 GeV. The azimuthal separation of the charged leptons in the transverse plane 
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Figure 7.31: Normalized distribution in the variable cos6*. 



(A0) is smaller than 45°. Finally, there should be no jet with pj^* larger than p^^° 

In Table 7.9 we report the corresponding cross sections in the case of p^^° = 30 GeV. 



a(fb) 


LO 


NLO 


NNLO 


tJ-F = f^R = Mh/2 


17.36 ±0.02 


18.11 ±0.08 


15.70 ±0.32 


IJ-F = l^R = Mh 


14.39 ± 0.02 


17.07 ±0.06 


15.99 ±0.23 


IJ'F = IJ'R = "^Mh 


12.00 ±0.02 


15.94 ±0.05 


15.68 ±0.20 



Table 7.9: Cross sections for pp ^ H + X ^ WW + X ^ Ivlu + X at the LHC when 
selection cuts are applied and = 30 Ge V. 

The cuts are quite hard, the efficiency being 8% at NLO and 6% at NNLO. The scale 
dependence of the result is strongly reduced at NNLO, being of the order of the error 
from the numerical integration. The impact of higher order corrections is also drastically 
changed. The i('-factor is now 1.19 at NLO and 1.11 at NNLO. As expected, the jet veto 
tends to stabilize the perturbative expansion, and the NNLO cross section turns out to be 
smaller than the NLO one. 

7.34.3 H ^ ZZ ^ e+e-e+e- 

We now consider the production of a Higgs boson with mass Mh = 200 GeV [25]. In 
this mass region the dominant decay mode is — > ZZ — > 4 leptons, providing a clean 
four lepton signature. In the following we consider the decay of the Higgs boson in two 
identical lepton pairs. When no cuts are applied, the NLO i^-factor \s K = 1.87 whereas 
at NNLO we have K = 2.26. We find that the interference contribution is smaller than 
1% in this region of Higgs boson masses. 
We consider the following cuts [36] : 

1. For each event, we order the transverse momenta of the leptons from the largest 
{pti) to the smallest {pta)- They are required to fulfil the following thresholds: 
PTi > 30 GeV PT2 > 25 GeV Pts > 15 GeV pr4 > 7 GeV ; 

2. Leptons should be central: \y\ < 2.5; 

^°Jets are reconstructed with the kr algorithm [39] with jet size D = 0.4. 
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3. Leptons should be isolated: the total transverse energy Et in a cone of radius 0.2 
around each lepton should fulfil Et < 0.05 px', 



4. For each possible e'^e pair, the closest (mi) and next-to-closest (771-2) to Mz are 
found. Then mi and m2 are required to be 81 GeV < mi < 101 GeV and 40 GeV 
< m2 < 110 GeV. 

These cuts are designed to maximize the statistical significance for an early discovery, but 
to keep the possibility for a more detailed analysis of the properties of the Higgs boson. 



The corresponding cross sections are reported in Table 7.10 



a(fb) 


LO 


NLO 


NNLO 


tJ-F = fJ'R = Mh/2 


1.541 ±0.002 


2.764 ±0.005 


3.013 ±0.023 


IJ'F = IJ'R = Mh 


1.264 ±0.001 


2.360 ±0.003 


2.805 ±0.015 


I^F = IJ'R = "^Mh 


1.047 ±0.001 


2.044 ± 0.003 


2.585 ±0.010 



Table 7.10: Cross sections for pp ^ H + X ^ ZZ + X ^ e+e'e+e- + X at the LHC 
when cuts are applied. 

Contrary to what happens in the H — > WW — > luli^ decay mode, the cuts are quite 
mild, the efficiency being 63% at NLO and 62% at NNLO. The NLO and NNLO K-factors 
are 1.87 and 2.22, respectively. Comparing with the inclusive case, we conclude that these 
cuts do not change significantly the impact of QCD radiative corrections. We also find 
that the effect of lepton isolation is mild: at NNLO it reduces the accepted cross section 
by about 4%. 



In Fig. |7.32| we report the pT spectra of the charged leptons. We note that at LO, 
without cuts, the pxi and pT2 are kinematically bounded by Mh/2, whereas pT3 < Mh/^ 
and PT4: < Mh/4:. Contrary to what happens in the ^ 77 decay mode (see Sect. 7.34.1[ ) 



the distributions smoothly reach the kinematical boundary, and no perturbative instability 
is observed beyond LO. 



7.35 Summary 

We have illustrated an extension of the subtraction formalism to compute NNLO QCD 
corrections to the production of high-mass systems in hadron collisions. We have con- 
sidered an explicit application of our method to the NNLO computation oi gg ^ H at 
the LHC, including the decay of the Higgs boson in all the relevant decay modes, namely, 
H — > 77, H — > WW lulv and H — > ZZ — > 4 leptons. We have presented few selected 
results, including kinematical cuts on the final state. In the case of the H ^ ■yy and 
H — > WW Ivlv decay modes, our computation parallels the one of Refs. [71 HQ], but 
it is performed with a completely independent method. In the quantitative studies that 
we have carried out, the two computations give results in numerical agreement. In our 
approach the calculation is directly implemented in a parton level event generator. This 
feature makes it particularly suitable for practical applications to the computation of dis- 
tributions in the form of bin histograms. The calculation is implemented in the numerical 
program HNNLO, which can be downloaded from |41| . 
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Figure 7.32: Tranverse momentum spectra of the final state leptons for pp ^ H + X ^ 
ZZ + X — > e~^e~e~^e~ + X, ordered according to decreasing pT, at LO (dotted), NLO 
(dashed), NNLO (solid). 
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QCD final states: resummation 
and Monte Carlo simulations 



A. Banfi 



8.36 Introduction 

Any short distance dominated cross section da in QCD can be written as a formal series 
in the QCD coupUng Os- 

da = dao + as ■ dai + ■ da2 + ■ ■ ■ (8.18) 

where dao is the leading order (LO) or Born contribution, dai the next-to-leading order 
(NLO) contribution, da2 the next-to-next-to-leading order (NNLO) and so on. In spite 
of the smallness of as, the coefficients of the expansion may be large. This happens 
typically when a process is characterised by two widely separated scales Q and Qq, where 
Q represents the hard scale of the process and Qo ^-n energy resolution. In this case large 



logarithms L = In Q/Qq arise at any order in the perturbative (PT) series eq. (8.18), and 
only after an all-order resummation can one give meaning to the PT expansion. In many 
cases resummation makes it possible to rewrite da as an exponent 

da = C{as) exp {Lgi{asL) + g2{asL) + asQsiasL)) + suppressed terms , (8.19) 

where gi resums the leading logarithms (LL, a^L""*"^), §2 the next-to-leading logarithms 
(NLL, a'^L^), gs the next-to-next-to-leading logarithms (NNLL, a"L"~^), and so on. 

The physical origin of large logarithms is the incomplete cancellation of soft and 
coUinear (SC) singularities between real and virtual contributions. In particular we dis- 
tinguish between double logarithms agL"^ arising from soft and coUinear emissions, and 
single logarithms a^L from hard coUinear or soft large-angle emissions. SC singularities 



factorise pi from hard matrix elements and build up the exponent in eq. (8.19). Finite 
virtual corrections and the exact treatment of the phase space in the SC limit give the 
multiplicative constant C{as), while hard emission contributions are suppressed by powers 
of Qo/Q. 

The above discussion can be visualised with the help of Lund diagrams. Consider for 
instance the well-known example of vector or Higgs boson production in hadron-hadron 
collisions. A generic contribution to the total cross section for this process in the SC limit 
is illustrated in fig. |8.33 Each dot represents an emitted parton, identified via its rapidity 



{rj, on the x axis) and transverse momentum {In kt/Q, on the y axis) with respect to the 
beam. The hard vertex is the origin of the axes. The yellow bands represent the coUinear 
limit 7] < ln(2Ei/kt), where Ei is the energy of emitting hard parton (leg) i = 1,2. Hard 
coUinear emissions (blue) are kinks on the two bands, since they reduce the emitting 
parton energy by a significant fraction. Soft large-angle emissions (red) are along the line 
77 = 0, while all remaining black dots are soft and collinear emissions. Since for fixed as 
emissions are distributed uniformly in In kt/Q and rj, an area in the picture corresponds 
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to a double logarithmic contribution, a line to a single logarithm, while points represent 
0{as) corrections. 

Virtual corrections are universal and can be shown to exponentiate The con- 
tribution of real emissions is instead observable dependent and can be represented as a 
vetoed region, where real emissions are forbidden, and only virtual contributions survive. 
According to the way the veto condition is imposed one can distinguish between 

1. inclusive observables: no hadrons are directly observed, QCD radiation is restricted 
via energy-momentum conservation. In this case, after an integral transform, real 
contributions exponentiate to all (logarithmic) orders. An example of an inclusive 
observable is the cross section for the production of a non-QCD particle, for instance 
a Higgs boson [3] ; 

2. final-state observables: one measures final-state hadron momenta, a typical exam- 
ple being event shape distributions and jet rates. In this case there is no general 
statement concerning the level of accuracy at which exponentiation holds. 



A lnk,/0 

hard scattering 



hard collinear 



soft collinear 



soft large angle 



hard collinear 



soft collinear 



In (1 „/Q 



Figure 8.33: The Lund diagram for a generic contribution to the total cross section for the 
production of a vector or Higgs boson at hadronic colliders. Here fiQ is an infrared cut-off 
and the hard scale Q is the boson mass. 

Inclusive observables have been discussed during the workshop by Massimiliano Grazzini. 
Here I will concentrate on final-state observables, introducing specific classes of observables 
and discussing to what extent a resummation for those observables is feasible. 



8.37 Global observables 

Let us consider a global variable that measures hadrons everywhere, for instance the thrust 
in e~^e~ annihilation, defined as 

r = max^ij%^. (8.20) 

" Ei \Pi\ 

The thrust, as all other event shapes, is a measure of the geometrical properties of hadron 
energy-momentum flow. For instance, for pencil-like events we have T ~ 1, for planar 
events T~2/3, while for spherical events T~ 1/2. This variable has the property of in- 
frared and collinear (IRC) safety, i.e. its value does not change after emission of extremely 
soft particles and/or quasi-collinear splittings. IRC safety is precisely the condition that 
ensures that observables related to the thrust, i.e. distributions and mean values, can be 
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computed using the quark-gluon language, in spite of the fact that the variable defini- 
tion involves hadrons, and the difference between parton and hadron level predictions is 
suppressed by powers of the hard scale Q. 

The basic quantity we are interested in is S(t), the probability that 1 — T < r, the 
differential distribution a^^da/dr being the derivative of S(r). Fixed order QCD pre- 
dictions are reliable as long as r is large, but fail as soon as events approach the Born 
limit T = 0. In the small r region one needs to resum logarithmic enhanced contributions 
to all orders, and the resulting resummed distribution has the same shape as the data. 
However, to get on top of the data, one needs to add a further correction that can be inter- 
preted as the difference between parton and hadron level. This hadronisation correction 
can be estimated using Monte Carlo (MC) event generators like herwig |i4j, pythia [5l|6] 
or ARIADNE [7], taking the ratio of the distributions obtained with MC's before and after 
hadronisation, and estimating hadronisation uncertainties using different event generators. 
This has lead to a successful description of IRC safe event shape distributions and jet rates 
in e~^e~ annihilation, giving one of the most accurate measurements of Og (see |8J for a 
recent review). The validity of this procedure relies strongly on the fact that MC event 
generators contain the physics that is needed to describe the main features of final-state 
observables. This statement is in general true for variables whose LL exponentiate, as we 
shall see in the following. 

Consider then a generic final-state variable V{{ki}), a function of final-state momenta 
{ki}, and its rate the probability that V{{ki}) <v. A generic contribution to 



can be represented by the Lund diagram on the left hand side of fig. 8.34 The grey area 
corresponds to the vetoed region where no real emissions are allowed, and only virtual 
corrections survive, and one can write in general: 

^(v) = e-^^^'^J^iv) , (8.21) 

where R{v) is the exponent representing virtual corrections up to the scale vQ, while 
real emission outside the vetoed region and the remaining virtual corrections build up the 
function J-{v). The variable V is said to exponentiate if all leading (double) logarithms 
are contained in the exponent R{v) and J-{v) is a pure NLL function, usually denoted by 
J-'{R'), with R'{v) = —vdR/dv. There are two basic conditions for this to happen, which 
go under the name of recursive infrared and coUinear (rlRC) safety conditions [9 : 

1. the variable must scale in the same fashion with multiple emissions as with a single 
emission. Formally, parametrising the momentum ki of each emission in terms of 
y(fcj)J^the value the variable V would have if only emission ki were present, and 
defining V{ki) = vQ, the first rlRC safety condition states that the following limit: 

lim nhKO.....M^C„)) 

has to be finite and non-zero. This ensures that the boundary of the vetoed region 



in fig. 8.34 does not change substantially whatever is the number of emissions; 



2. the variable's scaling property (8.22 ) must not be altered by the addition of extra-soft 



particles or by quasi-collinear splittings, formally: 

[ hm , hm] nkl{vCl),...,K{vCn),kn^,{vCn^,)) ^ ^ 

where the only non-trivial part of the commutator is the one where one takes the 
limit Cn+i after the limit v ^ 0, the other part being equal to the limit in 

For simplicity, we will always write V{ki, . . . , k„) instead of the more correct form V^({j5}, fci, . . . , 
where {p} denotes final-state hard parton momenta. 
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eq. (8.22) due to IRC safety. An analogous condition should hold also for coUinear 



splittings. This implies that in fig. 8.34 one can eliminate all emissions in grey, far 
from the boundary of the vetoed region, without altering the value of V, and for all 
emissions with V{ki) ~ V{ki, . . . ,kn) one can replace clusters of emissions close in 
rapidity with a single emission having the total momentum of the cluster. 




Figure 8.34: The Lund diagram for emissions contributing to a generic final-state observ- 
able (left) and its simplified version in case of a rlRC safe observable (right). 

In the end, for a rlRC safe variable, a generic contribution to T,{v) can be represented by 
a Lund diagram like the one on the right hand side of fig. 8.34[ where one has a vetoed 



area, giving rise to the LL exponent R{v), and real emissions contributing to J-{R') at 
NLL accuracy are both soft and coUinear, well separated in rapidity and confined in a 
narrow region of width Inl/e <C Inl/f close to the boundary of the vetoed area. The fact 
that emissions are well separated in rapidity makes it possible to exploit QCD coherence, 
and consider soft gluons as radiated independently (like in QED) from the hard legs. This 
simplification of multi-gluon soft matrix elements makes it possible to compute J'{R') with 
a MC procedure, where emissions are ordered in V{ki) = vQ, with Q < d-i, V{ki) is 
fixed at v, and the probability of emission of gluon ki coUinear to leg I with rapidity r]i 
and azimuth (with respect to leg £) is 

dPiUC^,V^,C^^,e))=R',^^^f(^-^y , E^^ = ^'- (8-24) 

The function T{R') can then be computed as the following average: 



where the limit v — > ensures that the result contains no NNLL contributions. Eq. (8.25) 
is an example of application of MC techniques used in parton shower event generators to 
obtain exact QCD results, and is one of the building blocks of the automated resummation 
program CAESAR [9j. 

We can now discuss what level of accuracy can be achieved by parton shower event 
generators. MC parton showers produce emissions in the whole of the phase space (for 



instance all emissions in the diagram on the left hand side of fig. 8.34), with approximated 
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matrix elements that are exact in the colhnear hmit but mistreat the soft large-angle 
region, both because they do not have full interference terms, and because they are correct 
only in the large- A'^c limit. However, for rlRC safe observables, LL and NLL contributions 
are determined only by emissions that are well separated in rapidity and close to the 
boundary of the vetoed region. MC event generators correctly describe such emissions, so 
that one expects that they reproduce not only LL, but most NLL contributions to rlRC 
safe observables. 



8.38 Non-global observables 



Non-global variables measure hadrons in a restricted part of the phase space. The most 
relevant example is the hadron transverse energy flow in a region J7 away from the hard 
jets, defined as [10] 



ti ) 



! 

Jo 



^^^^ 1 da 

dEt — 
a dEt 



].26) 



The distribution E(Qq) is sensitive only to soft emissions at large angles, so that LL are 
single logarithms with L = In Q/Qq. To achieve a LL resummation for T,{Qq) one 

cannot rely on an independent emission picture of QCD radiation, because configurations 
like the one in fig. |8.35 where one has two emissions close in rapidity and the kt of the 
harder emission can take any value from Q to Q^, give a single logarithm |10| . These 



In k,/Q 



In Qj/Q 



Figure 8.35: A configuration giving rise to non-global logarithms. The region $7 corre- 
sponds to a rapidity gap of width Arj between the jets. 

non-global logarithms are ruled by the correlated part of soft gluon emission matrix ele- 
ments, and can be resummed only in the large- A'c limit. Furthermore, due to the fact that 
one needs to take into account the exact form of 0, for an arbitrary number of emissions, 
resummation can be performed only numerically. In the MC procedure of ref. [10], emis- 
sions are ordered in the variable t = ^ "2!^'' ~ aslnQ/fcj, with ti > ti-i, and, given 
an ensemble of n— 1 colour connected soft gluons forming n dipoles, a softer gluon ki is 
emitted from dipole j with probability: 

dP{ki{ti,v^A^,j)) = 2CAdvi e - ^ dU e-2C^A^^-(*-*-i) , (8.27) 

where r]i and (pi are the rapidity and azimuth of emission ki in the frame in which partons 
forming dipole j are back-to-back, Arjj represents a collinear cutoff and Ar]tot = Ylj ^Vj- 
The MC starts from a quark-antiquark dipole, and continues emitting gluons as long as 



99 



one emission falls into Q. If one bins the value of t corresponding to this last emission, 
when the MC stops one has reconstructed dT,(t)/dt at LL accuracy, with tr-^aslnQ/Q^. 

Given the fact that the away-from-jet Et flow is used to tune MC models of the 
underlying event in hadron-hadron collisions, it is crucial to understand whether MC 
parton shower generators are able to describe If this is not the case, missing 

PT contributions, which are non-universal, will be included by tuning hadronisation and 
underlying event parameters, which are supposed to be universal. In ref. [llj it was first 
investigated the difference between a full treatment of soft radiation as in eq. (8.27), 
and an angular ordering (AO) approximation resulting from a free azimuthal average of 
dP{ki) (see for details). The difference between the full and the AO distributions was 
found to be less than 10% in the whole range of values of t relevant for phenomenology 
{t < 0.15, corresponding to Qn > IGeV at the LEPl energy), thus supporting the idea 
that MC generators implementing angular ordering should be able to give a reasonable 
description of T,{Qq). Then a direct comparison of a LL resummation for a~^da/dEt in 
e~^e~ annihilation and the corresponding predictions obtained with herwig and pythia at 
parton level was performed, the away-from-jet region O being a rapidity gap of width Ar] 
between the two jets. The comparison showed that herwig, whose evolution variable is 
the angle of each branching, is in good agreement with the LL resummation, discrepancies 
being less than 10%, both for large and small gaps. On the contrary, the old pythia 
shower [5], which uses the virtuality as an evolution variable and rejects a posteriori 
configurations not respecting AO, lies below the resummation, about 50% less at t = 0.15. 
The new pythia shower [6] is in agreement with the resummation for small gaps, Ar/ = 1, 
while for large gaps, Ar/ = 3, the two distributions have different shapes, and the point 
where they start to deviate seems to be exponentially related to the gap size. While it is 
known that the AO requirement in the old pythia shower places too strong a veto on soft 
emissions, this problem should be overcome with the new pythia shower. Therefore the 
origin of the disagreement is unclear and needs further investigation. 

Unlike rlRC safe global observables, resummation of non-global observables cannot 
rely on general approximations for soft radiation. For instance, the distribution for the 
away-from-jet hadron energy flow in eq. (8.26) has the form S(t) = exp[—R(t)]S{t), where 
S{t) contains non-global logarithms. Real emissions from the hard legs cancel with virtual 
corrections at scales larger than t, leaving only the exponent R{t). Therefore within LL 



accuracy there is no analogous of the function T{R') of eq. (8.25). This might be argued 
from the fact that for global observables J'{R') is sensitive only to soft and coUinear 
emissions, while here all relevant emissions are soft and at large angles. This however is 
not the case if one considers the away-from-jet energy flow of mini-jets, defined by replacing 
the sum over hadron with the sum over jets. This observable was introduced to reduce the 
impact of non-global logarithms, because emissions close in rapidity tend to be clustered 
together by jet algorithms ^L2i- Unfortunately, recombination spoils the complete real- 
virtual cancellation of primary emissions, so that one needs to introduce a new correction 
factor C{t) [13 . Also the statement that inter-jet energy flows are sensitive only to soft 
large-angle emissions might not be completely true. In ref. |13] it has been shown that, 
in hadron-hadron collisions, if one assumes that kt is the ordering variable for virtual 
corrections, a gluon outside the gap collinear to one of the incoming legs and emitting a 
softer gluon inside the gap gives a double-logarithmic contribution if it is accompanied 
by two-loop non-cancelling Coulomb phases. This higher order contribution to S(t) is 
a^vr^Ar^L^, representing thus a super-leading logarithm. The presence of super-leading 
logarithms depends on the ordering variable used for virtual corrections, so that at the 
moment one cannot clearly state whether these logarithms are actually present, and in 
case they are how they can be resummed to all orders. 
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8.39 Conclusions and outlook 



Final-state observables are very sensitive to QCD radiation. They are therefore extremely 
valuable tools to test our understanding of QCD dynamics. In particular, comparison of 
data, resummed analytical calculations and MC simulations can be used to improve our 
description of multi-soft gluon radiation. We hope that such a study, which was very 
successful at LEP, will continue at hadron colliders. 

One research direction could be then trying to use global rlRC safe observables for 
the tuning of models of the underlying event. The advantage of this class of observables 
is that they are correctly described by MC event generators at LL accuracy. This is 
not the case for the variable that is traditionally used for this purpose, the away-from-jet 
energy flow, which is instead non-global. However, globalness may be a problem at hadron 
colliders, since one cannot measure hadrons too close to the beam pipe. This problem has 
been addressed in ref. [15] , where one can find a list of global rlRC safe event-shapes and 
jet resolution parameters in hadronic dijet production. We look forward to experimental 
investigations in this direction. 

Concerning non-global logarithms, in hadronic collisions they appear in a variety of 
contexts. First of all, they contribute significantly to the distribution of the Et flow away 
from the hard jets, for which even a LL calculation is missing. Furthermore, non-global 
logarithms give NLL effects in distributions of non-global event shapes such as the ones 
defined in [16], which measure emissions only in a central rapidity region. These vari- 
ables are preferred from an experimental point of view, since hadron momenta can be 
measured by combining central tracker and calorimeter information, thus considerably 
reducing systematic uncertainties. Non-global logarithms will appear also in more inclu- 
sive distributions, like jet transverse momentum spectra or dijet azimuthal correlations. 
Theoretical studies in these directions are in progress. 
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Bottom-quark fragmentation: 
resummations and Monte Carlo 
simulations 



G. Corcella 

Heavy-quark phenomenology is currently one of the main fields of investigation in 
theoretical and experimental particle physics. In the following, we shall study i3-hadron 
production in e^e~ — > bb annihilation, top (t — > bW) and Higgs {H — > bb) decays. We 
will describe 6-quark production using resummed calculations as well as Monte Carlo 
generators, and get non-perturbative information by tuning hadronization models to ex- 
perimental data from SLD [1] and LEP [21 |3l S] . We shall also use a recently proposed 
non-perturbative model O El [7] , which includes power corrections in an effective strong 
coupling constant and does not introduce any further tunable parameter. We first consider 
bb pair production at the pole in the next-to-leading order (NLO) approximation, 

e+e- ^ Z^{q) ^ b{pb)b{pi) {g{pg)) , (9.28) 

and define the 5-quark energy fraction 

X, = (9.29) 



The energy spectrum of a massive b quark is given by: 



1 da as 



Pgg(xb)ln^ +^(X6) 



+ 0[^ \ , (9.30) 



where ao is the cross section of the Born process e^e'^ qq, p > 1, A{xb) is a function 
independent of the b mass, Pqq(xf,) is the Altarelli-Parisi (AP) splitting function. The large 
logarithm ~ asln(m^/?7i^) can be resummed by the use of the perturbative fragmentation 
formalism. 

Following [8j , the b spectrum is expressed as the convolution of a coefficient function, 
corresponding to the emission from a massless parton, and a perturbative fragmentation 
function D(mb, fJ-p), associated with the transition of a massless parton into a heavy b: 

1 dcTfe dz ^ "' 



-[Xh,mz,rnb, , ^ . 



f dz \ 1 dai 1 f^b ^ 

Jxt z [ao dz J \ z J 

+ OiimblmzY) . (9.31) 



In Eq. (9.31), dui/dz is the coefficient function for the production of a massless parton 



«, after subtracting the coUinear singularity in the MS factorization scheme. Neglecting 

-)M~ 



g —>■ bb splitting, i = b and expresses the fragmentation of a massless b into a massive 



b. 
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The perturbative fragmentation function follows the DGLAP evolution equations [21 
[To] and its value at a scale fip can be obtained once an initial condition at /Uqf is given. In 
[8] the NLO initial condition hqf, mj,) was calculated and its process-independence 

was established in [11]. Solving the DGLAP equations for an evolution from fi^p to fip, 
with a NLO kernel, allows one to resum leading (LL) agln"'(^|,///Qp,) and next-to-leading 
(NLL) ln"'~^(/i|.//iQ^) logarithms (collinear resummation) . Setting fiop ~ mb and 
Hp ~ mz, one succeeds in resumming the logarithms ln(m^/?n,^) appearing in the massive 



spectrum (9.30). 



Furthermore, both initial condition and coefficient function [8] present terms, ~ 1/(1 — 
Xb)+ and ~ [ln(l — Xb)/{1 — Xb)] + , which become large for Xb — > 1. The large-x;, limit 
corresponds to soft- or collinear-gluon radiation. Such contributions are usually resummed 
in Mellin moment space, where they correspond, at 0{as), to single (~ as' In A'^) and 
double (~ Qgln^ A^) logarithms of the Mellin variable N (soft or threshold resummation). 
In [11] threshold resummation was implemented in the NLL approximation; in [6] even 
large- A" NNLL contributions were resummed. To NNLL accuracy, terms ~ a'^ln'^^^ N 
(LL), ~ a^ln" A (NLL) and ~ agln"~^ A are kept in the resummed exponent. 

As for Monte Carlo parton shower algorithms, implemented in event generators such as 
HERWIG fl2l and PYTHIA [13] , they rely on the universality of the elementary branching 
probability for soft or collinear radiation. Referring, e.g., to parton cascades in e~^e~ — > qq 
processes, the probability of soft or collinear emission reads: 

.P=|5^PW..M||^. (3.32) 



In (9.32) P{z) is still the AP splitting function, z is the energy fraction of the emitted 
parton, is the shower ordering variable. 

In HERWIG, is an energy- weighted angle, equivalent to angular ordering in soft 
limit [13] . In PYTHIA [13] , is the momentum squared of the radiating parton, with an 
option to veto branchings that do not fulfil angular ordering. Moreover, the latest PYTHIA 
version offers, as an alternative, the possibility to order final-state showers according to 
the transverse momentum of the emitted parton with respect to the emitter [18] . It was 
found out [19] that the PYTHIA transverse-momentum-ordered showers yield a better 
treatment of angular ordering, although its implementation is still not as accurate as 
it is in HERWIG. Hereafter, we shall use PYTHIA 6.220, whose cascades are ordered 
in virtuality with the option to reject non-angular-ordered showers turned on, and the 
version 6.506 of HERWIG. In ( |9.32[ ) As{Ql,Ql) is the Sudakov form factor, expressing 



the probability of evolution from Qf to Q2 with no resolvable emission. 



For multiple emissions, iterating the branching probability (9.32) allows one to resum 
soft- and collinear-enhanced radiation: as discussed in [20], parton shower algorithms 
resum leading logarithms in the Sudakov exponent, and include a class of subleading 
NLLs as well. 

Calculations based on the perturbative fragmentation formalism are supplemented by 
non-perturbative fragmentation functions to yield hadron spectra. Up to power correc- 
tions, the i?-hadron spectrum reads: 

i^(x^,Q,m,) = i r '^'p^iz,Q,mb)D--r-^), (9.33) 
a axB cr J^^ z az \ z / 

where is the B energy fraction and D^p the non-perturbative fragmentation function. 
In the following, we shall use the the NLO/NLL perturbative calculation in Ref. |11] along 
with Kartvelishvili model |15| : 

Z)"P(x; 7) = (1 + 7)(2 + 7)(1 - x)x\ (9.34) 
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and fit 7 to experimental data. 

We shall also use a non-perturbative model, based on an extension of [5], consisting in 
including power corrections in an effective coupling constant. Such a model was presented 
in detail in Refs. [6l [7]; here we just point out that it employs the following effective 
coupling constant: 

as{k^) = — f ds Disc, ^Mz^^ (9.35) 
where in the integrand function one sets: 







1 f°° ds 

^^''^^^ ^ 2TTi J s + k"^ ^^^^"^ OLs{-s), (9.36) 

with as being the standard coupling constant. qs(A;^) and asikP') are usually called [B] 
time-like and space-like effective coupling constants, respectively. As discussed in [6l[7], 
one can prove that the effective coupling as{k'^) is free from the Landau pole and that 
at small momenta it includes power-suppressed effects. Also, it is remarkable that the 



model based on Eq. (9.35) does not introduce any extra tunable parameter. As in |^, the 
effective coupling will be implemented in the NNLO approximation, and used along with 
a calculation based on the perturbative fragmentation formalism, with NLO coefficient 
function and initial condition, NLL DGLAP evolution and NNLL threshold resummation. 

As far as HERWIG and PYTHIA are concerned, their parton showers terminate when 
a scale Qq, of the order of 1 GeV, is reached. The hadronization is simulated according to 
the cluster |16J and string |l7j models, respectively. 

We shall now consider data on i?-hadron spectra at the pole, collected by the SLD 
[I], OPAL [3] and ALEPH [2] collaborations. The ALEPH data contain only 6-flavoured 
mesons, the OPAL and SLD ones a small fraction of baryons as well. As in [6l |2T], when 
using resummed calculations, we limit the comparison to < 0.92, in order to avoid 
very large-xs data, where our computation is still not completely reliable and the spectra 
become negative or oscillate. We convolute the NLO/NLL calculation of [H] with the 



Kartvelishvili model (9.34) and find that, in the considered range, the best fit is obtained 
for 7 = 17.178 ± 0.303, with xV^of = 46.2/53. 

The effective-coupling model does not have any free parameter, but nonetheless in [6] 
all quantities entering in the perturbative computation were varied within conventional 
ranges, in order to gauge the theoretical uncertainty on the prediction. For the sake of 
brevity, we do not present here all the plots shown in we just point out that the best 
comparison with the data is obtained for /ioF = "i-fe/2, where /Uqf is the factorization scale 
in the initial condition of the perturbative fragmentation function. We obtain x^/dof = 
103.0/54, which is a quite reasonable value, since we are not tuning any parameter. It was 
also shown in Ref. [6j that setting ^qf = fni, and = 5.3 GeV, a mass value characteristic 
of a meson, leads to an excellent description of the ALEPH data, with x^/dof = 11.9/16. 

For the purpose of HERWIG and PYTHIA, as in plj, we fit only the non-perturbative 
parameters of the cluster and string models; in fact, the default parametrizations yield 
rather poor fits of the 6-fragmentation data, as we obtain x^/dof = 739.4/61 for HERWIG 
and xVdof = 467.9/61 for PYTHIA. In HERWIG, we change CLSMR(l) and CLSMR(2), 
ruling the Gaussian smearing of the hadron direction with respect to the original con- 
stituent quarks; PLSPLT(2), which determines the mass distribution of 6-flavoured clus- 
ter decays; DECWT, affecting the relative weight of decuplet and octet baryons; and 
CLPOW, to which the heavy-cluster yield and the baryon/meson ratio are sensitive. The 
fitted values are: CLSMR(l) = 0.4 (default 0), CLSMR(2) = 0.3 (0), PSPLT(2) = 0.33 
(1), DECWT = 0.7 (1), CLPOW = 2.1 (2). After the tuning, the agreement with the 
data is still not very good, but it is much better than with the default parametrization: 
X^/dof = 222.4/61. In PYTHIA, we modify the values of the fragmentation parameters 
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Figure 9.36: B-hadron spectrum measured by ALEPH, OPAL and SLD experiments, along witli 
the HERWIG and PYTHIA predictions, after fitting the hadronization parameters to the data. 
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Figure 9.37: As in Fig. 9.36 but comparing the data witli the NLO/NLL calculation using the 
Kartvelishvili hadronization model and the NLO/NNLL computation with the effective strong 
coupling constant. 



PAR J (41) and PARJ(42), which control the a and h parameters of the Lund symmetric 
fragmentation function, and PARJ(46), which modifies the endpoint of the Lund function 
according to the Bowler hadronization model [22*. Our tuning gives: PARJ(41) = 0.85 (de- 
fault value 0.3), PARJ(42) = 1.03 (0.58), PARJ(46) = 0.85 (1). After our tuning, PYTHIA 
matches the e^e~ data very well, and we obtain x^/dof = 45.7/61 from the fit. We have 
checked that our tuning works well also for the new model implemented in PYTHIA 6.3, 
which orders parton showers in transverse momentum; we found x^/dof = 46.0/61 from 
the comparison with the xb data. In Fig. 9.36 we compare LEP and SLD data with 
HERWIG and PYTHIA; in Fig. 9.37 we present the experimental spectra along with the 
NLO/NLL calculation using the Kartvelishvili model and the NLO/NNLL one with the 
analytic coupling constant. All approaches use the best-fit parametrizations. 

We note in Fig. 9.36 that PYTHIA, after the tuning, gives a good description of the 
experimental spectra, while HERWIG's distribution is broader, below the data around 
the peak and above them at small xb- From Fig. 9.37 we learn that the NLO/NLL 



calculation using the Kartvelishvili model reproduces the data quite well, but it becomes 
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Figure 9.38: _B-hadron spectrum in top quark decay according to HERWIG, PYTHIA and the 
NLO/NLL calculation which includes non-perturbative corrections via the Kartvelishvili hadroniza- 
tion model. 



negative at large xb- The plot relying on the effective-coupling model lies above the data 
around the peak and approaches zero more rapidly at large xb- In any case, even this 
result is acceptable, considering that, when modelling power corrections by means of the 
effective coupling constant, we are not tuning any free parameter to the data. 

Using the fits to LEP and SLD data, we can predict the i?-hadron spectrum in other 
processes, such as top-quark decay (t — > bW) and the decay of the Standard Model Higgs 
boson H 



bb. In Figs. 



9.38 



and 



9.39 



we show the predictions yielded by HERWIG, 
PYTHIA and the resummed calculation based on the perturbative fragmentation approach 
for i?-hadron production in top {t bW) and Higgs {H bb) decays. We parametrize 
cluster, string and Kartvelishvili models using the best fits to LEP and SLD data. As for 
the resummation, we use the NLO/NLL calculations in Refs. [231 121] for top decay and 
in Ref. [25] for H ^ bb processes. In our plots we have set the top and Higgs masses to 
the values mt = 175 GeV and rriH = 120 GeV. The results in Figs. 9.38 and 9.39 exhibit 



similar features to the comparison presented in Figs. 9.36 and 9.37: PYTHIA and the 



NLO/NLL calculation using the Kartvelishvili model are in good agreement, while the 
spectra yielded by HERWIG show some discrepancy, as they are broader than the other 
two predictions and lie below them at small xb and above at large xb- 

Finally, we present the same comparison in moment space, where the moments of the 
differential cross section are defined as follows: 



UN 



dz z 



).37) 



In Ref. [1], the DELPHI collaboration presented the first five moments for B produc- 
tion in e^e~ annihilation. From the point of view of resummed calculations, working in 
moment space [2B] presents several advantages. In A^-space, convolutions become ordi- 
nary products, and the relation between parton- and hadron-level cross sections becomes: 

~ ^N^^^ where cj^ and are the moments of the b and B cross sections, and 
is the A^-space counterpart of the non-perturbative fragmentation function. Therefore, 
there is no need to assume any functional form for the non-perturbative fragmentation 



function in x^-space. The results of our A^-space analysis are summarized in Table 9.11 



We note that, after the fits to LEP and SLD data, HERWIG and PYTHIA agree with the 
DELPHI moments, within the experimental uncertainties. As for the calculation based on 
the effective coupling constant, it is able to reproduce the experimental moments within 
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Table 9.11: Moments from DELPHI [4 , and A^-space results in e+e annihilation, Higgs (H) 
and top (t) decay, according to NLO/NLL calculations, HERWIG (HW) and PYTHIA (PY). Also 
presented are the TV-space e~^e~ results obtained using the effective coupling constant as along 
with the theoretical errors 





{x) 








e+e data 


0.715±0.005 


0.540±0.006 


0.424±0.007 


0.341±0.006 


e+e- NLL a% 


0.780 


0.644 


0.548 


0.476 




0.917 


0.839 


0.773 


0.716 


e+e- HW 


0.711 


0.535 


0.418 


0.335 


e+e- PY 


0.716 


0.541 


0.424 


0.340 


e+e- NNLL+Q!5 


0.687 ±0.040 


0.5019 ±0.047 


0.381 ±0.046 


0.298 ±0.046 


t-dec. NLL 


0.723 


0.556 


0.443 


0.363 


i-dec. HW 


0.733 


0.570 


0.461 


0.381 


t-dec. PY 


0.722 


0.559 


0.449 


0.369 


H-dec. NLL 


0.695 


0.517 


0.402 


0.321 


H-dcc. HW 


0.684 


0.504 


0.388 


0.308 


iJ-dec. PY 


0.688 


0.508 


0.391 


0.310 
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the theoretical errors which were calculated in [H] and quoted in Table 9.11 As done for 
the x^-space analysis, we also present in Table [9TT] the moments of the differential width 
for the production of B hadron in top or Higgs decays, using the moments taken from 
the fits to the DELPHI data. 

In summary, we studied bottom-quark fragmentation in e'^e" annihilation, top and 
Higgs decays using resummed calculations and the HERWIG and PYTHIA parton shower 
models. We fitted the Kartvelishvili, cluster and string models to i?-hadron spectra mea- 
sured at LEP and SLD, and then predicted the i?-energy distribution in other processes. 
We also presented the results yielded by a model which incorporates non-perturbative 
power corrections via an effective strong coupling constant. The analysis was finally ex- 
tended to Mellin moment space. We believe that the study here presented can be a useful 
starting point to address 6-quark fragmentation at present and future colliders, as it sets 
some benchmarks for the hadronization models which are typically used along with Monte 
Carlo generators and resummed computations. 



109 



Bibliography 



[1] SLD Collaboration, K. Abe et al., Phys. Rev. Lett. 84 (2000) 4300. 

[2] ALEPH Collaboration, A. Heister et al., Phys. Lett. B512 (2001) 30. 

[3] OPAL Collaboration, G. Abbiendi et al., Eur. Phys. J. C29 (2003) 463. 

[4] DELPHI Collaboration, G. Barker et al., DELPHI 2002-069, CONF 603. 

[5] D. Shirkov, Nucl. Phys. Proc. Suppl. 152 (2006) 51. 

[6] U. Aglietti, G. Corcclla and G. Fcrrcra, Nucl. Phys. B775 (2007) 162. 

[7] G. Corcella and G. Ferrera, JHEP 0712 (2007) 029. 

[8] B. Mele and P. Nason, Nucl. Phys. B 361 (1991) 626. 

[9] G. AltarelH and G. Parisi, Nucl. Phys. B126 (1977) 298. 

[10] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438; 
L.N. Lipatov, Sov. J. Nucl. Phys. 20 (1975) 95; 
Yu.L. Dokshitzer, Sov. Phys. 46 (1977) 641. 

[11] M. Cacciari and S. Catani, Nucl. Phys. B617 (2001) 253. 

[12] G. Corcella, I.G. Knowles, G. Marchesini, S. Moretti, K. Odagiri, P. Richardson, 
M.H. Seymour, B.R. Webber, JHEP 0101 (2001) 010. 

[13] T. Sjostrand, S. Mrcnna and P. Skands, JHEP 0605 (2006) 026. 

[14] G. Marchesini and B.R. Webber, Nucl. Phys. B238 (1984) 1; ibid. B310 (1988) 461. 

[15] V.G. Kartvelishvih, A.K. Likehoded and V.A. Petrov, Phys. Lett. B78 (1978) 615. 

[16] B.R. Webber, Nucl. Phys. B238 (1984) 492. 

[17] B. Andcrsson, G. Gustafson, G. Ingclman, T. Sjostrand, Phys. Rept. 97 (1983) 31. 

[18] T. Sjostrand and P.Z. Skands, Eur. Phys. J. C39 (2005) 129. 

[19] A. Banfi, G. Corcella and M. Dasgupta, JHEP 0703 (2007) 050. 

[20] S. Catani, G. Marchesini and B.R. Webber, Nucl. Phys. B349 (1991) 635. 

[21] G. Corcclla and V. Droningcr, Nucl. Phys. B730 (2005) 82. 

[22] M.G. Bowler, Z. Phys. Cll (1981) 169. 

[23] G. Corcella and A.D. Mitov, Nucl. Phys. B623 (2002) 247. 

[24] M. Cacciari, G. Corcella and A.D. Mitov, JHEP 0212 (2002) 015. 

110 



[25] G. Corcella, Nucl. Phys. B 705 (2005) 363 [Erratum-ibid. B 715 (2005) 609]. 
[26] M. Cacciari and P. Nason, Phys. Rev. Lett. 89 (2002) 122003. 



Ill 



Monte Carlo simulations of 
top-quark pair production in 
hadronic collisions 

M. Treccani 

10.40 Introduction 

In view of the starting of the LHC and the accumulated statistics at Tevatron, there ap- 
pears the need for further improvement in the accuracy of theoretical predictions. One 
of the most interesting fields refer to the class of events with multiple final states, giving 
rise to multiple jets with complicated topologies. There exists different strategies to tackle 
this problem, with distinct features and points of strength. The main problem is how to 
consistently compose the contributions due to Matrix Element (ME) calculations with the 
contributions of the Monte Carlo showering codes (MC), in order to exploit their comple- 
mentariety and avoid at the same time the so-called double counting phenomenon [1]. 
One of these strategies, known as MC@NLO, put the emphasis on achieving the next-to- 
leading-order (NLO) accuracy in the description of the inclusive rates for a given final 
state F, accompanied by the exact leading-order (LO) description of the emission of one 
extra jet (-F-|-jet). For a detailed explanation of this approach and its implementation in 
several cases, see [21 [3l Hj. 

One alternative approach relies on a consistent leading-logarithmic (LL) accuracy in the 
prediction of a final state F accompanied by a varying number of extra jets. The removal 
of double counting of jet is achieved by the so-called matching algorithm for matrix ele- 
ments and parton shower. It is understood that the matching algorithm approach cannot 
improve the intrinsic LL accuracy of the predictions; however it will give a better accuracy 
in the prediction of the observables more sensible to the production of two or more jets in 
addition to F. 

In this note, we study in detail the so-called MLM matching [5.,, 6] embedded in the the ME 
generator ALPGEN [7i in the ti pair production at hadron colliders. First we will address 
its stability w.r.t its internal parameters, and after we will perform detailed numerical 
comparison between MLM matching and the MC@NLO code. 

In particular, in Section 2 we will perform some robusteness test on the ALPGEN cal- 
culations, comparing predictions obtained with different parameters and discussing the 
related uncertainties. Section 3 covers the detailed comparison between ALPGEN and 
MC@NLOpredictions, and in Section 4 we will present our conclusions. 

10.41 Consistency studies of the matching algorithm 

In this section we study the overall consistency of the matching algorithm applied to the 
case of tt final states. We shall consider tt production at the Tevatron (pp collisions at 
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y/S = 1.96 TeV) and at the LHC {pp collisions at = 14 TeV). 

The generation parameters for the light partons are defined by the following kinemati- 
cal cuts: the default values for the event samples at the Tevatron (LHC) are given by: 
^min_2Q ^^Q-j QcVand i?„i„=0.7 (0.7), while they are considered only in the geometrical 
region defined by r] < 4(5). 

The top particle is assumed to be stable, and therefore all jets coming from the de- 
cay of top quarks are neglected. For the shower evolution we use HERWIG, version 
6.510 [H]. We stopped the evolution after the perturbative phase, in order to drop down 
all the common systematics that could smooth out any possible discrepancy between 
the various simulations. For all generations we chose the parton distribution function 
set MRST2001J^, with renormalization and factorization scales squared set equal to: 
A'r — f^F ~ Z]i=t fjets ["''i + (Pt)^]- observables are built out of the partons emerging 
form the shower in the rapidity range < 6 and adopting the cone algorithm GETJET|10]. 
The jet cone size is set to i?cone = 0.7 and the minimum transverse momentum to define 
a jet at the Tevatron(LHC) is 15(20) GeV . 

Having defined the environmental parameters of these studies, we then explore the sys- 
tematic uncertainties due to the variation of the internal parameters. These uncertainties 
reflect the underlying fact that this approach relies on the LO evaluation of the hard ME 
and on the LL accuracy in the removal of double counting and in the description of the 
shower evolution. In this section we shall show that the size of the resulting uncertainties 
is consistent with what can be expected in such a LL approach in the case of ti production. 



To our analysis, the important feature of the whole procedure is the presence of two 
set of parameters: the generation cuts and the matching cuts (see 13 |6]). The first set 
is necessary to avoid the Infra-Red (IR) and collinear singularities:/?™*" , the minimum 
transverse momentum of the extra-parton(s) to be generated, and Rmin, the minimum 
separation between extra-partons in the (r/, (p) plane. Along with these parameters, there 
exist an analogous set, but with slightly different meaninings : the matching cuts E^"° and 
Rmatch- It's worth to stress that the latter parameters are necessary to effectively separate 
the phase space, but the prediction should be stable against (slight) modifications of them, 
together with the choice of the particular cone jet algorithm adopted in the matching 
procedure. 

In our examples here we consider two independent variations of the generation and of 
two of the matching cuts, as in table 10.12, keeping fixed our definition of the physical 
objects (the jets) and of the observables. 

Then we proceed to study some distributions for the Tevatron, showing the observables 
dominated by contributions with up to 1 hard parton in fig. 10.40 and those relative to 



multijet final states in fig. 10.41 We find that these distribution are stable against rea- 
sonable variations of the internal parameters, with relative differencies confined well below 
few percents, both in matching and generation parameter variations. Angular observables, 
such as Ai?, are more sensible, since they are directly related to the matching variables, 
and their agreement is within 10%. 

The analysis at the LHC, which will not be shown here, leads to qualitatively and 
quantitatively similar results. 



10.42 Comparisons with MC@NLO 

We shall now compare in detail the description of ti events as provided by ALPGEN and 
MC@NLO. For consistency with the MC@NLO approach, where only the 0{al) ME effects 
are included, we use ALPGEN samples obtained by stopping the ME contributions only 
to 1 extra-parton besides the ti pair. This strategy allow to highlight the different features 
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Figure 10.40: Comparison between the tliree alternative sets of generation (left) 



and matching (right) parameters given in table 10.12 at the Tevatron 
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Figure 10.41: Comparison between the three alternative sets of generation and 
matching parameters given in table [T0.12| for multijet distributions at the Tevatron. 
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Table 10.12: Variations of the generation and matcliing parameters used for the 
study of the systematics. 



of the two alternative approaches applied to same set of contributions. It is understood 
that a homogeneous comparison can only be done through the introduction of a proper 
K factor, determined by the ratio of the total rates of the two predictions. We adopt the 
same simulation setup, modifying only the same factorization and renormalization scale 
in order to match MC@NLO's default: 

The upper two rows of plots in figs. 10.42 refer to inclusive properties of the tt sys- 
tem, namely the transverse momentum and rapidity of the top and anti-top quark, the 
transverse momentum of the tt pair, and the azimuthal angle At/)" between the top and 
anti-top quark. The overall agreement is good, once ALPGEN is corrected with the proper 
K-factor (1.36 for the Tevatron, and 1.51 for the LHC), and no large discrepancy is seen 
between the two descriptions of the chosen distributions. The most significant differencies 
(10 to 20%) are seen in the p^'' distribution, ALPGEN's one being slightly softer. 

The study of jet quantities reveals instead one important difference: the rapidity of 
the leading jet, yi, is different in the two descriptions, where MC@NLO exhibits a dip at 
Ui = 0. This difference is particularly marked at the Tevatron, but is very visible also at 
the LHC. This is shown in the right figure of the third row in fig. 10. 42 [ Visible differences 
are also present in the distribution of the 1st and 2nd jet separation in (r/, (p) space, Ai2i^2- 

To understand the difference in the rapidity distribution, we look in more detail in 
fig. 10.43 at some features in the MC@NLO description of the leading jet. For the px of the 
leading jet, Pt,i, we plot separately the contribution from the various components of the 
MC@NLO generation: events in which the shower is initiated by the LO tt hard process, 
and events in which the shower is initiated by a tt + q{g) hard process. In this last case, we 
separate the contribution of positive- and negative- weight events, where the distribution 
of negative events is shown in absolute value. The plots show that for MC@NLO the 
contribution of the ti + q{g) hard process is almost negligible over most of the relevant 
range and becomes appreciable only for very large values of Pt,i- This hierarchy is stronger 
at the LHC than at the Tevatron. 

Upper set of fig. 10.44 shows the various contributions to the rapidity distribution 
7/i for different jet Pt thresholds. It appears that the Ui distribution resulting from the 
shower evolution of the tt events in MC@NLO has a strong dip at ?/i=0, a dip that cannot 
be compensated by the more central distributions of the jet from the tt+q{g) hard process, 
given its marginal role in the overall jet rate. 

That the dip at yi=0 is a feature typical of jet emission from the tt state in HERWIG 
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Figure 10.42: Comparison of ALPGEN (histogram) and MC@NLO (plot) distri- 
butions, at tlic Tcvatron. The ALPGEN results are rescaled to MC@NLO, using 
the K factor of 1.36. The relative difference (MC@NLO-ALPGEN)/ ALPGEN) is 
shown at the bottom of each plot. 
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Figure 10.43: Contributions to the transverse momentum of the leading jet in 
MCQNLO. Tevatron (left) and LHC (right). 




Figure 10.44: Rapidity of the leading jet j/i at Tevatron for for various jet 
thresholds. Upper set: MC@NLO, with partial contributions. Central set: HER- 
WIG. Lower set: comparison between ALPGEN, MC@NLO, and the parton level 
predictions 
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is shown in central set of fig. 10.44 obtained from the standard HERWIG code rather than 
from MC@NLO. We speculate that this feature is a consequence of the dead-cone descrip- 
tion of hard emission from heavy quarks implemented in the HERWIG shower algorithm. 
To complete our analysis, we show in lower set of fig. |10.44 the comparison between the 
ALPGEN, MC@NLO and the parton-level yi spectra, for different jet Pt thresholds. We 
notice that at large where the Sudakov effects that induce potential differences between 
the shower and the PL results have vanished, the ALPGEN result reproduces well the PL 
result, while still differing significantly from the MC@NLO distributions. 



10.43 Conclusions 

The study presented in this paper examines the predictions of ALPGEN and its matching 
algorithm for the description of tt-|-jets events. Several checks of the algorithm have shown 
its internal consistency, and indicate a rather mild dependence of the results on the param- 
eters that define it. The consistency of the approach is confirmed by the comparison with 
MC@NLO. In particular, inclusive variables show excellent agreement, once the NLO/LO 
K factor is included. 

Instead we found a rather surprising difference between the predictions of two codes for the 
rapidity distribution of the leading jet accompanying the ti pair. In view of the relevance 
of this variable for the study at the LHC of new physics signals, it is important to further 
pursue the origin of this discrepancy, with independent calculations, and with a direct 
comparison with data. Preliminary results obtained with the new positive-weight NLO 
shower MC introduced in [TTl [12] and presented in this Meeting [13] , appear to support 
the distributions predicted by ALPGEN. 



119 



Bibliography 



[1] S. Hoche, F. Krauss, N. Lavesson, L. Lonnblad, M. Mangano, A. Schalicke and 
S. Schumann, arXiv:hep-ph/0602031. 

[2] S. Prixione and B. R. Webber, arXiv:hep-ph/0601192. 

[3] S. Frixione and B. R. Webber, JHEP 0206 (2002) 029 [arXiv:hep-ph/0204244]. 

[4] S. Frixione, R Nason and B. R. Webber, JHEP 0308 (2003) 007 [arXiv:hep- 
ph/0305252]. 

[5] M.L. Mangano, presentation at the FNAL Matrix Element /Monte 

Carlo Tuning Working Group, 15 Nov 2002, http://www- 
cpd.fnal.gov / personal / mrenna / tuning/nov2002 / mlm.pdf . 

[6] M. L. Mangano, M. Moretti, F. Piccinini and M. Treccani, JHEP 0701, 013 (2007) 
[arXiv:hep-ph/0611129]. 

[7] M. L. Mangano, M. Moretti, F. Piccinini, R. Pittau and A. D. Polosa, JHEP 0307 
(2003) 001 [arXiv:hep-ph/0206293]; 

F. Caravaglios, M. L. Mangano, M. Moretti and R. Pittau, Nucl. Phys. B 539 (1999) 
215 [hep-ph/9807570]. 

[8] G. Marchesini and B. R. Webber, Nucl. Phys. B 310 (1988) 461. G. Marchesini, 
B. R. Webber, G. Abbiendi, I. G. Knowles, M. H. Seymour and L. Stanco, Comput. 
Phys. Commun. 67 (1992) 465. G. Corcella et ai, JHEP 0101 (2001) 010 [hep- 
ph/0011363]. G. Corcella et al, "HERWIG 6.5 release note," arXiv:hep-ph/0210213. 

[9] A. D. Martin, R. G. Roberts, W. J. Stirling and R. S. Thorne, Eur. Phys. J. C 23 
(2002) 73 [arXiv:hep-ph/0110215]. 

[10] F.E. Paige and S.D. Protopopescu, in Physics of the SSC, Snowmass, 1986, Colorado, 
edited by R. Donaldson and J. Marx. 

[11] P. Nason, JHEP 0411 (2004) 040 [arXiv:hep-ph/0409146]. 

[12] P. Nason and G. Ridolfi, JHEP 0608 (2006) 077 [arXiv:hep-ph/0606275]. 

[13] P. Nason, presented in this Meeting 



120 



Phenomenology of the Standard 
Model Higgs boson at the LHC 



G. Corcella and D. Rehuzzi 

The Higgs boson plays a crucial role in the Standard Model (SM) of electroweak 
interactions, as it is responsible of the mechanism of mass generation. However, this 
particle has not yet been experimentally discovered. Searches for the Higgs boson will be 
one of the main goals of the Large Hadron Collider, which will be capable of exploring the 
Higgs mass spectrum from 100 GeV to about 1 TeV. In order to accurately perform such 
searches, the use of precise QCD calculations and reliable Monte Carlo event generators 
will be mandatory. 

We study several observables related to the phenomenology of the SM Higgs boson 
and compare the predictions yielded by HERWIG and PYTHIA the two most 
popular event generators, as well as QCD computations resumming the large logarithms 
appearing in the Higgs transverse momentum spectrum [3]. In fact, some differences 
between HERWIG and PYTHIA are to be expected, since they implement parton showers 
[H [5], matrix-element matching [6 , 7J and hadronization [HI |9] in a different fashion (see 
|10] for some discussions and comparison between the two Monte Carlo codes). Hereafter 
we shall use the HERWIG 6.510 and PYTHIA 6.403 versions. As parton distribution 
functions (PDFs), the leading-order (LO) CTEQ6L1 set [TT will be employed. 

We shall consider Higgs production via gluon-gluon fusion (GGF), which is the domi- 
nant channel at the LHC, and through vector-boson fusion (VBF). In HERWIG, the same 
user-defined process simulates both gg ^ H and qq H, while in PYTHIA the two sub- 
processes can be run separately. In the following, for the sake of comparison, we shall use 
a modified version of HERWIG, with the qq ^ H subprocess turned off. We checked that, 
with the PDF set [11] , HERWIG simulates about 6-7% of events according to qq — > H. 

Indeed, if we compare the HERWIG and PYTHIA total cross sections for gg — > H, 
using the default parametrizations, we find meaningful discrepancies, about 15-20% for a 
Higgs mass 110 GeV < niH < 190 GeV. We investigated the possible causes determining 
such differences and understood that they are mostly due to the value of the strong 
coupling constant asimn) implemented in the two programs. Both codes tuned the QCD 
parameter A to LEP data, along with other quantities, such as the shower cutoff, quark 
and gluon effective masses and the hadronization non-perturbative parameters. However, 
such fits led to pretty different results for the strong coupling constant at the Z mass: 
as{mz) ~ 0.116 in HERWIG and as{mz) ^ 0.127 in PYTHIA. While the value of 
HERWIG is consistent with the world average, i.e. as{mz) = 0.118±0.002 [12], PYTHIA 
uses a somewhat higher value. Different values of asirnz) clearly do not affect the total 
LO e^e" qq cross section, but do have an impact on the LO gg ^ H one, which 
is C'(a|(m/f)). Another difference between the two codes is the implementation of the 
QCD beta function, which is at two loops in HERWIG and at one loop in PYTHIA. At 
hadron colliders, while HERWIG still uses its own value of asimz), PYTHIA employs, 
for the hard process and the initial-state parton cascade, the same value of as{mz) as 
the one in the chosen PDF. The LO CTEQ6L1 set uses as^mz) = 0.130, even larger 
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than the PYTHIA default, and much above the world average. In any case, we point 
out that, since HERWIG and PYTHIA yield only LO rates, they should not be used 
to calculate the Higgs production total cross section, which is currently available up to 
next-to-next-to- leading order (NNLO) |13j . Parton showers predict instead more reliably 
differential distributions, which exhibit milder dependence on as{rnz), since they are 
equivalent to resummations (see for some comparison between Monte Carlo algorithms 
and resummed computations). Throughout our analysis we shall nonetheless use the 
same value of asimn) in both HERWIG and PYTHIA. For example, we can employ, 
as a reference point, the value as^mz) = 0.130, as in the parton density [11], leading 
to as{mH) — 0.123 using one-loop evolution from mz to mn- In HERWIG, in order 
to increase asijnz), we shall have to increase the parameter QCDLAM (default value 
0.18 GeV), which roughly corresponds to QCD quantity A in the MS renormalization 
scheme at high momentum fractions [T3]. However, QCDLAM cannot be made arbitrarily 
large without modifying other HERWIG non-perturbative parameters, such as the shower 
cutoffs for gluons and quarks, namely VGCUT (default 0.10 GeV) and VQCUT (default 
0.48 GeV). We found that a possible combination yielding 05(771^) = 0.123 in HERWIG 
is the following: 



QCDLAM = 0.378 GeV ; VGCUT = 1.00 GeV ; VGCUT = 1.50 GeV. 



(11.38) 



It can be seen that using the parametrization (11.38) in HERWIG the two programs give 
approximately the same rates, as shown in Table 11.13[ where we quote the total cross 
section given by defaul t PYTHIA (PY) and HERWIG (HW), and by HERWIG tuned 
according to Eq. ( 11.38| ) (HW*). Of course, one can also modify the PYTHIA parameters 
in order to have the same as{m}{) as in HERWIG. Before moving on, we point out that 
we are aware that changing the default parameters of a Monte Carlo program can be 
dangerous, as this may likely spoil the agreement with the data taken into account in the 
fits. In fact, we are not recommending that one should use a different parametrization. 



such as Eq. (11.38), but just trying to understand the reason of the discrepancy and 



whether it is possible to improve the agreement between HERWIG and PYTHIA. 



Table 11.13: Cross sections for Higgs production in the gluon-fusion channel, according to 
PYTHIA (PY) and HERWIG (HW), using their default parametrization, and according 



to HERWIG, tuned as in Eq. (11.38) (HW 





a(PY) 


a (HW) 


a (HW*) 


110 GeV 


20.7 pb 


16.6 pb 


20.2 pb 


130 GeV 


15.5 pb 


13.2 pb 


15.5 pb 


150 GeV 


12.2 pb 


10.2 pb 


12.2 pb 


170 GeV 


10.3 pb 


7.9 pb 


10.7 pb 


190 GeV 


7.9 pb 


6.6 pb 


8.1 pb 



As far as the decays of the Higgs boson are concerned, we investigated both the total 
width as well as the partial rates into a few given channels.. The total rates yielded by 



HERWIG and PYTHIA, even after tuning HERWIG as in Eq. (11.38), are a bit differ- 



ent, as quoted in Table 11.14 Several issues contribute to the discrepancy exhibited by 
Table [11.14 For example, the treatment of the decay H — > ZZ, on which our study 
will be later on mostly concentrated, is quite different in the two codes. While PYTHIA 
implements the general case, where both Z bosons are allowed to be off-shell, in HERWIG 
at least one Z is forced to be on-shell. We checked that if we selected events where both 
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Table 11.14: Total Higgs decay width according to PYTHIA and HERWIG, using its 
default parameters and the tuning (11.38|). 



niH 


r (PY) 


r (HW) 


r (HW*) 


110 GeV 


2.5 MeV 


3.5 MeV 


3.0 MeV 


130 GeV 


4.4 MeV 


5.4 MeV 


4.8 MeV 


150 GeV 


15.8 MeV 


16.4 GeV 


15.6 GeV 


170 GeV 


355.2 MeV 


328.6 MeV 


337.3 MeV 


190 GeV 


981.6 MeV 


919.3 MeV 


919.6 MeV 



Z's are on the mass shell, e.g. within five widths, HERWIG and PYTHIA rates would 
agree up to a good accuracy level. 

We present in Table |11.15| the H ^ ZZ rates according to PYTHIA, HERWIG and 
HDECAY, a code, based on Ref. [15], computing the total and partial Higgs decay rates. 



possibly including higher-order radiative corrections. We note reasonable agreement be- 
tween PYTHIA and HDECAY, which also permits that both Z's are off-sheh, while HER- 
WIG yields slightly lower widths. 

Among the other Higgs decay modes, major differences between HERWIG and PYTHIA 
are present especially in the channels into heavy quarks, such as H ^ cc or H ^ bb, where 
the discrepancies between the two default codes can be up to ~ 50%. Considering, e.g., the 
decay H — > bb, both HERWIG and PYTHIA implement the 6-quark MS mass rhh{mH), 
which is an appropriate mass definition for bb production at the Higgs mass scale. However 
HERWIG, unlike PYTHIA, also includes the leading-logarithmic (LL) resummation of the 
large contributions ~ a§(m/i-) ln"(m///m;,) and NLO corrections to Removing 
such higher-order corrections, and still using the tuning (11.38), we expect that HER- 



WIG and PYTHIA should agree. The partial rates of the two Monte Carlo generators 
can again be compared with the results of the HDECAY code. In the H ^ bb mode, 
HDECAY includes the NNLO corrections to the MS 6-quark mass, the NLO ones to the 
massive rate T(H bb) near threshold, and even the NNNLO ones, in the massless ap- 



proximation rrih <C niH, far above threshold. We present in Table 11.16 the widths for 
H —>bb according to HERWIG, PYTHIA and HDECAY for few values of niH and using 
ctsi'iTT'H) = 0.123 everywhere. HERWIG yields the largest rate, even above the HDECAY 
result; hence, we may conclude that some of the higher-order corrections that HDECAY 
implements, whereas HERWIG does not, have negative sign. We finally remark that, 
unlike HERWIG, PYTHIA includes the modes H ^ gg and H ^ jZ*. 



Table 11.15: Width T{H for different values of nin according to PYTHIA, 



HERWIG with the parametrization (11.38) and HDECAY, for several values of the Higgs 
boson mass 



niH 


r (PY) 


r (HW*) 


r (HDECAY) 


110 GeV 


0.012 MeV 


0.011 MeV 


0.012 MeV 


120 GeV 


0.053 MeV 


0.051 MeV 


0.054 MeV 


130 GeV 


0.186 MeV 


0.173 MeV 


0.189 MeV 


140 GeV 


0.530 MeV 


0.503 MeV 


0.541 MeV 


150 GeV 


1.357 MeV 


1.311 MeV 


1.374 MeV 



As anticipated, in our phenomenological analysis we shall consider SM Higgs produc- 
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Table 11.16: As in Table 



11.15| but for the decay H — > hh. 



rriH 


r (PY) 


r (HW*) 


r (HDECAY) 


110 GeV 


1.88 MeV 


2.46 MeV 


2.23 MeV 


130 GeV 


2.17 MeV 


2.82 MeV 


2.55 MeV 


150 GeV 


2.46 MeV 


3.18 MeV 


2.85 MeV 


170 GeV 


2.74 MeV 


3.51 MeV 


3.15 MeV 


190 GeV 


3.02 MeV 


3.84 MeV 


3.44 MeV 



tion in GGF and VBF and concentrate ourselves mostly on the decay channel H ZZ 
AL We shall study the following distributions: the Z- and ff-boson mass spectrum, the 
Higgs and Z transverse momentum [qT,H and qT,z) and pseudorapidity [rjH and r]z). Such 
spectra are presented in Figs. 11.45f|Tl.49 for a Higgs mass tuh = 130 GeV. 

The ruH distributions of HERWIG and PYTHIA (Fig. |11.45[ ) look compatible. On the 
contrary, the fact that PYTHIA allows both Z's to be off-shell has an evident impact on 



Fig. 11.46: in the intermediate Z-mass range, say 40-80 GeV, where both Z's are off-shell, 
and PYTHIA yields more events. Once again, if we set a filter allowing only Z's near the 



mass shell, the discrepancies in the intermediate mz range will disappear. In Fig. 11.47 
we instead compare the Z transverse momentum distributions: in GGF PYTHIA predicts 
more events than HERWIG at large qT, while in VBF the two codes roughly agree. As 
we shall discuss later in more detail, for the time being, the default version of PYTHIA 
includes matrix-element corrections to Higgs production in GGF, while HERWIG does not. 
Such corrections are responsible of the simulation of a few events with a Higgs of large 



transverse momentum, whose decays still yield Z's at large qT- In Fig. 11.48 we instead 
present the Higgs pseudorapidity distributions: we clearly note an asymmetry, about 5%, 
for VBF in HERWIG, with more events simulated at positive rather than negative rjH- 
PYTHIA yields instead a symmetric spectrum. This asymmetry exhibited by HERWIG is 
currently under investigation [16] and should be clarified in a forthcoming publication |17] . 
As for the r]H distribution in GGF, both HERWIG and PYTHIA spectra are symmetric, 
although some discrepancy is still present, with PYTHIA leading to more events around 
77// = 0. The HERWIG and PYTHIA Z-pseudorapidity spectra, presented in Fig. |11.49 
are instead rather similar. 
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Figure 11.45: Higgs invariant-mass distribution for H ZZ decays and Higgs production 
in gluon-gluon (left) and vector-boson (right) fusion, according to HERWIG (solid line) 
and PYTHIA (dashes). 
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Figure 11.46: As in Fig. 11.45 but showing the Z-mass distribution. 




Figure 11.47: Z transverse momentum spectrum in GGF (left) and VBF (right). 




Figure 11.48: Higgs pseudorapidity distribution in GGF (left) and VBF (right). 



We finally wish to present results on the Higgs transverse momentum distribution, 
which has been thoroughly investigated using Monte Carlo generators as well as re- 
summed calculations. In fact, the NLO Higgs qt spectrum exhibits contributions, ~ 
a^ln^{m'jj/q'^), with k < n + 1, which are large for small values of the Higgs transverse 
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Figure 11.49: As in Fig. 11.48 but presenting the Z pseudorapidity. 



momentum. Small values of qt correspond to initial-state soft or collinear parton radi- 
ation. Such logarithms have been resummed in [3^ up to next-to-next-to-next-to-leading 
logarithmic accuracy (NNLL) in a Sudakov-like exponential factor. The authors of Ref. [3] 
also released a computing code, named HqT, implementing numerically their resummation. 

In detail, the LLs are ~ agln"+^(m2^/g2 ), the NLLs ~ ln"(m2^/4), the NNLLs 
~ a^lii'^~^ {m^ / q^) and so forth. While at small qt the resummation works fine, at large 
qt one should rely on fixed-order calculations. The HqT code allows one to match the 
resummation with the fixed-order spectrum up to NNLO accuracy As for HERWIG and 
PYTHIA, their standard algorithms [4, 5j are reliable for soft or collinear parton radiation, 
i.e. for Higgs production at small qx- In fact, it can be shown [14, that parton showers 
include all soft /collinear leading logarithms, plus some NLLs as well. Higgs production 
at large qT corresponds instead to hard or large-angle initial-state radiation: therefore, it 
cannot be simulated by standard parton shower algorithms, but must be described by the 
use of the exact tree-level NLO matrix-element |18| [T9] . 

In order to allow hard and large-angle radiation, HERWIG and PYTHIA have been 
provided with matrix-element corrections [6l [7], although their actual implementation is 
indeed somewhat different. HERWIG splits the phase space into two regions: a region 
corresponding to soft /collinear emission, where one uses the parton shower approximation, 
and a region associated with hard or wide-angle radiation, the so-called 'dead zone' of the 
standard algorithm, where the exact tree-level amplitude is used. Moreover, the exact 
matrix element is also employed to correct the radiation in the HERWIG parton-shower 
region any time an emission is capable of being the hardest so far, i.e. it has the largest qx 
with respect to the radiating parton. Unlike HERWIG, PYTHIA uses instead the parton- 
shower approximation in all physical phase space and corrects with the exact amplitude 
only the first branching [7]. As discussed in [6], correcting only the first emission will 
however lead to an unphysical dependence of the hard-emission probability on the infrared 
cutoff, appearing in the Sudakov form factor. 

The latest version of PYTHIA [2] does include matrix-element corrections to Higgs 
production in gluon-gluon fusion, though treating the top quark in the loop in the infinite- 
mass limit. The official version of HERWIG [II does not include yet the corrections to 
Higgs production, although a preliminary code is available, based on the paper [20] , which 
extends the earlier work for Drell-Yan processes [21] . With respect to PYTHIA, Ref. [20] 
fully includes top-quark mass effects and possibly corrects even to qq H, whenever 



^^Throughout this paper, by NLO we always mean corrections of ©(aal), which are NLO for the total 
//-production cross section. For the qr spectrum, such contributions are instead LO. 
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this subprocess is turned on. We also compare the Monte Carlo spectra with the HqT 
code, which we run in the NLL approximation, matched to the NLO total cross section 
in the point qt = mn- The results are presented in Fig 11.50 where we also show the 



result yielded by the HERWIG 'Higgs+jet' process, where the hard-scattering process is 
always generated according to one of the tree-level corrections to gg H, i.e. gg Hg, 
qg Hq, etc. We set in all codes asimu) to the same value, however, at large qx still 
different scales // for the strong coupling are used. In fact, by default, HERWIG uses 

the Higgs transverse mass /x = ^J~(^^^\-rr}^ , which is a reasonable scale since the dead 
zone corresponds to gr > mn', PYTHIA uses instead fj, = q^', HqT sets fi = mn- For 
the sake of comparison, we set // = in ah codes for qx > mn- this way, all spectra 



should roughly agree at large transverse momentum. From Fig. 11.50 we learn that the 
Higgs-l-jet spectrum roughly agrees with HqT at large qx- The small discrepancy can be 
due to the fact that HqT also includes bottom quarks in the loop and to different choices 
of the PDF factorization scale, which is set to fip = in HqT and fip = qx in HERWIG 
and PYTHIA. The HERWIG and PYTHIA qx spectra, though in agreement with each 
other, are a bit below the other two curves, even though the same tree-level NLO matrix 
element has been used. Such a discrepancy obviously needs further investigation [T7j . 
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Figure 11.50: Transverse momentum distribution of the Higgs boson, produced in gluon fu- 
sion, according to HqT (solid line), HERWIG (dashed), PYTHIA (dot-dashed) and the HERWIG 
'Higgs-l-jet' process (dotted). 

In summary, we performed a study on SM Higgs boson production and decay at the 
LHC, mainly using PYTHIA and HERWIG, the two most popular Monte Carlo generators. 
We found several differences between the two codes: we understood the causes of most 
discrepancies, whereas a few are still under investigation. In a forthcoming work [T7] we 
shall present an even more detailed analysis of Higgs boson phenomenology, where we shall 
also study the impact of the NLO corrections to Higgs production in GGF, implemented 
in the MC@NLO code |22]. 
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12.44 Introduction 

According to the Minimal Supersymmetric Standard Model (MSSM), two isospin Higgs 
doublets have to be introduced. After electroweak symmetry breaking, five Higgs scalar 
mass eigenstates remain: one CP-odd neutral scalar boson A, two charged scalars H^, 
and two CP-even neutral scalars h and H. At the tree level, the Higgs sector is completely 
defined by only two parameters: they are usually chosen as the ratio of the vacuum 
expectation values (tan/3 = V2/V1) and the Higgs boson A mass {Ma}- The tree level 
hierarchies {Mh < Mz, Ma < Mh and Mw < Mh) are modified by large radiative 
corrections: the leading one-loop correction is proportional to mf and the upper bound of 
Mh is shifted to Mh < 135 GeV/c^ . 

Varying Ma in the range 91 GeV/c^ < Ma < 1 TeV/c^ we can distinguish three 
different regimes: 

• Decoupling regime. If M4 » M"^"-^ then the Higgs bosons H, A and are very 
heavy and almost degenerate in mass, while h has a mass very close to MJ^""^ and 
becomes SM-like. H and A, besides their masses, are degenerate also in width and 
cross section. 

• Low Ma regime. If Ma < MJ^'^^ the behavior of the two CP-even neutral Higgs 
bosons h and H is swapped with respect to the decoupling regime: h is almost 
degenerate in mass, width and cross section with A, and H is the SM-like Higgs, 
with a mass close to Mj^'^^. 

• Intense coupling regime. This occurs for Ma ~ MJ^^""^ and high tan /3 jll [2] and 

it leads to similar, but not degenerate, masses for the three neutral Higgs bosons. 
This property, in principle, allows to detect the three neutral Higgs separately. 

The LEP experiments have excluded the Higgs masses Ma < 91.9 GeV/c^ , Mh^n < 
91 GeV/c2 [3] and M^± < 78.6 GeV/c^ g]. 

In this report the discovery potential of the MSSM neutral Higgs boson with the CMS 
detector at LHC is presented. These analysis are also described in [5]. 



12.45 Neutral Higgs bosons searches 

The production cross-section for the MSSM neutral Higgs bosons is strongly dependent 
on the value of the tan /? parameter. All neutral MSSM Higgs production cross sections 



including NLO QCD corrections are shown in Fig, 12.51 
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12.45.1 Large tan (3 

In the region with large tan /3 values (> 15) the neutral Higgs bosons are mainly produced 
in association with b-quarks: pp — > qq/gg h/A/H + bb. The presence of a bb pair is 
important to suppress the very large background from Drell-Yan processes. The Higgs 
bosons mainly decay in a bb pair (90%) and in a rr pair (10%). In CMS six channels have 



been studied: 






• A/H^ 






• A/H^ 


TT - 


e + jet + X 


• A/H^ 


TT - 


fi + jet + X 


• A/H^ 


TT - 


->■ jet + jet + X 


• A/H^ 


TT - 


■^e + n + X 


• A/H^ 


bb 





The muon final state, with respect to the other channels, has a much lower branching 
ratio (~ 3 X 10~^), but the event is very clean and Higgs masses and widths can be 
reconstructed precisely. Moreover it is possible to exploit the theoretical relation between 
the Higgs decay width and tan (3 {Tfj oc tan^ /5) to perform a direct measurement of this 
latter quantity. 

The rejection strategy is mainly based on identification of isolated muons and on b- 
tagging. This latter selection is particular important: b jets from signal events are mainly 
produced in the forward region with lower pT with respect to the b jets coming from tt 



background. Figure 12.52 (left) shows the reconstructed dimuon invariant mass for signal 
and background. 

The tau channels, on the other hands, have a better signal to background ratio and 



can reach larger discovery region in the plane (Myi,tan/3), as can be seen in Fig. 12.54 



(left). Fig. 12.52 (right) shows the reconstructed invariant mass for signal and background 
that can be obtained with the A/H tt ^ e+jet channel. Indeed, despite the escaping 
neutrino, the Higgs boson mass can be reconstructed also for these channels, exploiting the 
coUinearity approximation: the neutrino is assumed to be emitted along the r direction. 

Finally, the bb channel must take into account the huge QCD background and, to 
perform a discovery, one needs to know in advance masses and widths of the Higgs bosons. 
Thus this channel can be considered mainly as a cross-check for the discovery. 

12.45.2 Small tan /3 

Concerning low tan (3 values, the dominant neutral MSSM Higgs production mechanism 
is the gluon fusion gg — > h/A/H, which can be mediated by top and botton loops (as in 



the SM case), but also by stop and sbottom (Fig ]12.5T ). 
In CMS two channels have been investigated: 

• A^ Zh^ £+£-bb 

. A/H^x^X^^^i + E^''' 

The first channel provides an interesting way to detect A and h simultaneously. The 
cross section increases with decreasing of tan /?, while the mass range is niz + rnh < ttia < 
2mtop- However results are strongly dependent on the MSSM parameters /i and M2, 
because the Higgs boson decay ^ — > XiX? may become dominant (the best results being 
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M^^ (GeV/c') m,, (GeV/c^) 



Figure 12.52: Reconstructed dimuon mass for the main background and for the signal 
sample with Ma = 150 GeV/c^ and tan/3 = 40 for the A/H ^ fifi channel (left) and 
with Ma = 300 GeV/c^ and tan (3 = 25 for the A/ H ^ tt ^ e + jet channel (right). 



obtained for large values of fx and M2). Fig. 12.53| (left) shows the discovery region for 
this channel for 30 and 60 fb~^ of integrated luminosity. 

To increase the discovery region in the low and intermediate region of tan (3, the sec- 
ond channel has been studied which takes into account the decay modes of the neutral 
Higgs bosons to supersymmetric particles. The final state studied by this channel is par- 
ticular clean (four leptons plus missing transverse energy). The analysis is performed in 
three benchmark points of the minimal Super Gravity constrained version of the MSSM 
(mSUGRA): these points are obtained varying the parameters niQ and mi for tan/3 = 5, 
10, signifi) = + and Aq = 0. Fig, 12.53 right) shows the discovery region in the (mo, mi) 

I V 2 

plane, for an integrated luminosity of 30 fb~^ . 
Fig, 



12.54 



summarize the 5a discovery region that can be obtained with 30 fb of 



integrated luminosity. 



12.46 Charged Higgs bosons searches 

Three channels, depending on the final state and on the Higgs boson mass, have been 
investigated for charged Higgs bosons: 

• TVr, with Mh < Mtop 

• TUr, with Mh > Mtop 

• ^ tb, with Mh > Mtop 

For Higgs boson masses below Mtop, the main production mechanism is through top 



decay, t H^b, and the branching ratio of the tu channel is about 98% (Fig, 12.55 'left)). 
The study is performed considering the leptonic decay of the W: ti — > H^W^bb — > 



TUT-iu^bb, T — > hadrons. The discovery region, as can be seen in figure 12.55 'right), covers 
almost the entire allowed region in the {Ma, tan (3) plane. 

If Mh > Mtop, the charged Higgs bosons are mainly produced in association with a 
top-bottom pair, gg — > tbH^. The final state for this channel is very clean {H^ — > r^z/, 
r — > hadrons + v and — > jj) and, after the selection cuts, almost background free. 
The characteristics for this channel is the presence of large missing transverse energy and 
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Figure 12.53: (left) The 5a discovery region contours for 30 and 60 fb^^ integrated 
luminosity for the A — > Zh channel. The effect of underestimation or overestimation of 
the background systematic uncertainty can be seen in the courve of 30 fb~^ . (right)The 
5a discovery region contours for 30 fb~^ integrated luminosity for the A/H — > X2X2 
channel. 




Figure 12.54: The 5a discovery regions for the neutral Higgs bosons (p {4'=h, H, A) 
produced in the association with b quarks pp — > bbi;^ with the (p /i/U and (p ^ tt decay 
modes (left plot) and for the light, neutral Higgs boson h from the inclusive pp h+X 
production with the h — > 77 decay and for the light and heavy scalar Higgs bosons, h and 
H, produced in the vector boson fusion qq qqh(H) with the h(H) — > rr — > £+jet decay 
(right plot). The m™*^^ scenario is used. 



133 




Figure 12.55: Branching ratios for charged Higgs boson decaying to different final states 
for tan P = 20 (left). The 5(T-discovery regions for the charged Higgs boson with the tv 
decay mode with 30 fb~^ of integrated luminosity (right). 



the T helicity correlations favouring the — > decay over the —>■ t^u decay. A 



Two 



large sector of the (M^,tan/3) plane can be covered (Fig. 12.55 'right)). 

Finally, for masses above Mtop + Mhottom, the channel — > tb opens up. 
production mechanism are considered: 

• gb^ tH^ ttb W+W~bbb qq'^iv^bbb 

• gg^ tH^b ttbb W+W^bbbb qq'fii^^bbbb 

Unfortunately no sensitivity is obtained in the MSSM parameter space with this anal- 
ysis, due to the large background and the resulting large effects of systematic uncertainties 
on its knowledge. 



Fig, 12.55 right) summarize the 5a discovery region that can be obtained with 30 fb 



of integrated luminosity. 



12.47 Conclusions 

Many channels have been studied to estimate the discovery potential of MSSM Higgs 
bosons at CMS. A large area in the {M^, tan f3) plane will be explored: the most promising 
channels are (j) ^ tt {(j) = h,A,H) and rv. 
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Physics studies at the LHC with 
PHANTOM 

G. Bevilacqua 

13.48 Introduction 

To shed light on the nature of Electroweak Symmetry Breaking (EWSB) is one among 
the main purposes of the LHC experiments. The Standard Model provides the most 
economical explanation in terms of the Higgs mechanism. The remarkable agreement 
between its predictions and precision electroweak data is consistent with a light Higgs 
boson with a mass lower than 200 GeV[T], while direct searches have established a lower 
bound of 114.4 GeV[2! at 95% CL. 

The Higgs mechanism is essential to ensure the renormaliz ability of the Standard 
Model. Moreover, without a Higgs the theory is not self-consistent and violates pertur- 
bative unitarity at high energies, which implies that for several processes the amplitude 
grows indefinitely with energy. This unphysical behaviour is particularly manifest in the 
scattering of longitudinally polarized vector bosons, whose cross section goes over the uni- 
tarity limit at about one TeV. If the Higgs hypothesis is not realized, effects of new physics 
are expected to restore unitarity at this scale. 

On the other side, if a massive Higgs boson exists, a resonance will be observed in 
the VV invariant mass spectrum in correspondence of the Higgs mass. Vector Boson 
Fusion represents the second most important contribution to the cross section for Higgs 
production at LHC, moreover H VV is the preferred decay channel above about 150 
GeV. Final states in which the vector bosons decay into four leptons or two leptons plus 
two jets provide one of the cleanest signatures for Higgs detection. 

It should be mentioned that even if a light Higgs will be discovered, the gauge boson 
scattering could nonetheless deviate from the Standard Model predictions in the high 
energy domain. Several models have been proposed, whereby the Higgs is a composite 
Goldstone boson associated to some new strong dynamics active at the TeV scalejH 1^. 
The l^y-scattering cross section is expected to be enhanced at high Myv as a consequence 
of the new strongly interacting sector. 

In brief, Boson-Boson Scattering has a great potential for probing the mechanism of 
EWSB at LHC, independently of its particular realization. Unfortunately, no beam of 
on-shell bosons will be available. The only consistent way to extract some information 
from data is to rely upon a complete calculation of final states with six fermions. 

A Monte Carlo event generator, PHANT0M[5 , has been recently developed to provide a 
complete description at 0{a%j^,f) + 0{a^^a1) of all six-parton final states in the Standard 
Model, including for the first time all EW and mixed EW-I-QCD contributions at this 
perturbative order. Exact matrix elements at tree level are evaluated efficiently. A good 
coverage of phase space is achieved thanks to a new approach [6^ which combines the 
best features of multichannel with the adaptivity of the VEGAS algorithm[7J. Apart from 
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Boson-Boson Scattering and Higgs search, the program gives access to the full tree-level 
description of tt production which is important both in itself and as a background in 
connection with the previous topics. 

In this note we present the first parton-level studies based on PHANTOM. The analysis 
is focused on Boson-Boson Scattering and Higgs search via Vector Boson Fusion in the 
semi-leptonic channels, that is final states characterized by four quarks and two charged 
leptons or one charged lepton plus one neutrino. A light Higgs scenario is compared with 
no-Higgs {Mh — > oo) results, to be intended at this stage as a minimal benchmark in the 
search for new physics beyond the Standard Model. 

13.49 Boson-Boson Scattering signature and its irreducible 
background 

The problem of how to define a signal for Boson-Boson Scattering has been debated since 
a long time. Ideally, one should isolate the contribution of diagrams like the one depicted 
m Fig{l3^a), in which the incoming partons emit two space-like vector bosons which 
scatter among themselves and finally decay into fermions. The resulting amplitude should 
then be deconvoluted from the PDF's and finally projected to on-shell bosons. In practice, 
this procedure turns out to be problematic due to gauge invariance and cancellation effects. 
In fact the subset of boson-boson fusion diagrams is not separately gauge invariant and 
there are relevant negative interferences |8j with the rest of the diagrams, some of which are 



illustrated in Fig, 13.56 ^b,c,d). Previous studies have evidentiated substantial differences 



among results based on this approximation and exact tree-level predictions on the VV 
mass spectrum[9l [ipj, showing that a complete evaluation of all contributions at 0{a%j^,j) 
is necessary. 







(a) (b) (c) (d) 

Figure 13.56: The Boson-Boson Scattering topology (a) and other typical contributions 
to the amplitude of a six-fermion final state (b,c,d). 

In the two-step approach we follow, we first study the characteristics of Boson-Boson 
Scattering from a six-fermion point of view. As no contribution to VV scattering can be 
found outside 0{a^j^j), it is convenient to start with a complete study of all processes at 
this perturbative order. In our intention, this represents a benchmark for physics analyses 
and does not claim to be fully realistic. With the help of appropriate selection cuts, we 
isolate a sample of events which best exemplifies the signature of the physical process 
under investigation. This means that candidates for top or three-boson production are 
subtracted to the advantage of events characterized by two tag jets widely separated in 
pseudorapidity and two couples of jets and leptons resonant atWovZ mass. Subsequently, 
further procedures are examined to enhance the contribution of longitudinally polarized 
vector bosons, which are most sensitive to effects of EWSB. The inclusion of the relevant 
QCD background represents the step forward which is essential to get meaningful results 
at parton level and completes our analysis. 

Some examples of contributions to the QCD irreducible background are illustrated in 



Fig, 13.57 Top production, either from single t or tt, is a huge source of background for this 
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kind of studies. In particular tt is dominated by ©(a^j^^al) contributions due to diagrams 
shown in Fig, 13.57 ^a,b). VV + 2 jets at 0(a|J^,^a|) is another important background with 
a signature characterized by two outgoing vector bosons, quite similar to the scattering 
one and therefore particularly challenging to subract. 




(a) (b) (c) (d) 

Figure 13.57: Examples of contributions to the QCD irreducible background: ti production 
(a,b) and VV + 2 jets (c,d) 

In the case of totally leptonic channels, i.e. final states with four leptons, 0{a'^j^ja^) 
encompasses the full QCD background at tree level for six-fermion analyses. On the 
contrary, semileptonic channels get additional contributions from ©(al^^a^) diagrams 
which are responsible for the V + 4 jets background. The latter must be included as well 
in a complete analysis and can be covered by other Monte Carlo event generators, such 
as AlpGen[ll or MadEvent |T2] . It should be noticed nevertheless that these contributions 
have quite different kinematical features with respect to the scattering signature, therefore 
they are expected to be easily suppressed by means of appropriate selection cuts. Of course, 
the final word is left to the results of a complete study. 



13.50 0{a^j^j^j) results in the semileptonic fi^fi channel 

This section is devoted to the analysis at 0{a%^) of all the processes of type qq — > 
qqqqH'^fJ'~- All our results have been obtained using the CTEQ5L[13^ PDF set with scale 

1 ^ 

Q^ = M^ + -^^p^T„ (13.39) 

i=l 

where pTi denotes the transverse momentum of the i-th final state particle. 

As already mentioned, the choice to start from pure EW results is motivated by the fact 
that no QCD diagram contributes to the scattering topology. At this stage we are mainly 
concerned with isolating as much as possible the VV scattering signature in the spirit of 
the discussion of the previous section. Following this approach, the selection procedure 
we adopt in this section makes use of flavour information and is not fully realistic. Results 
are found anyway to be not too sensitive to the details of the selection cuts, as shown 



in Fig, 13.58 We compare different methods of reconstructing the VZ invariant mass. 
Results based on flavour selection and predictions obtained using cuts on the invariant 
mass only are shown. The two distributions differ by about 20% at small invariant masses 
but agree quite nicely above 800 GeV, showing that results based on flavour information 
are not seriously degraded when selection procedures closer to the actual experimental 
practice are adopted. 

It is not possible to separate the contribution of the various subprocesses {ZW — > ZW, 
WW ZZ, ZZ ZZ, tt, single t) to a given six-fermion final state. Still, it is possible 
to isolate different subgroups of reactions including different subsets of diagrams. For 
example, only the subprocess ZZ ZZ can be identified inside the scattering set of the 
diagrams describing the reaction uu uussix^ fj,~ , hence the latter will be unambigously 
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Figure 13.58: Invariant mass distribution of the lepton pair and the two jets from boson 
decay for the no-Higgs case. The sohd(dashed) hne is obtained identifying the two most 
central jets as the vector boson decay products before(after) top vetoing. The dotted line 
is obtained requiring the correct flavour content for the jets identified as decay products of 
both the vector boson and the top. The dot-dashed lines is produced using solely invariant 
mass informations to identify the vector boson and the top decay products. 



assigned to the ZZ — > ZZ subgroup. If more than one subprocess contributes to a given 
reaction, the same will appear in different subgroups. This is the case, for instance, with 

ZZ\ The plot on the left in Fig, 13.59 shows 



ud udssfi+fi- {ZZ ZZ, WW 

the invariant mass distribution of the two most central quarks and of the two leptons for 
all reactions which contain the different subprocesses as well as the distribution for the 
complete set of processes. A light Higgs with Mh = 150 GeV has been assumed. In order 
to comply with typical acceptance and trigger requirements, the standard acceptance cuts 



in Tab, 13.17 have been applied. It should be clear that the total cross section in Fig, 13.59 
is smaller than the sum of the cross sections for the various groups as a consequence of 
double counting. 



Acceptance cuts 


Selection cuts 


Pt{(^) > 20 GeV 


\M{bW) - Mtop\ > 15 GeV 


|ry(£±)| <3 


\M{e+e~) - Mz\ < 10 GeV 


E{q) > 30 GeV 


\M{qiq2 from W) - Mvk| < 10 GeV 


PT{q) > 20 GeV 


\M{qiq2 from Z) - Mz| < 10 GeV 


\v{q)\ < 5 


\M{qm) - Mv^-I > 10 GeV 


M {£+£-) > 20 GeV 


|M(q3Q4) -Mz\>W GeV 


M{qq) > 60 GeV 




Ar/(tag-quarks) > 3.8 



Table 13.17: Standard acceptance and selection cuts applied in results on the fJ-^iJ-~ chan- 
nel. M{qiq2 from W(Z)) is the invariant mass of the two quarks with the correct flavour 
content to be produced in a W{Z) decay. If more than one combination of two quarks sat- 
isfies this requirement, the one closest to the corresponding central mass value is selected. 
M{q^q4) denotes the invariant mass of the two remaining quarks. 



For the selected type of final states, top background enters only through single t pro- 
duction. Top candidates are rejected requiring a b (b) quark in the final state together 
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Figure 13.59: Invariant mass distribution of the two charged leptons and the two most 
central quarks, for different sets of processes. The plot on the left includes the set of 
acceptance cuts described in the left part of Tab 13.171 In the plot on the right a further 
cut for vetoing top production is applied. 



with two other quarks of the right flavour combination to be associated to a (VV~) 
decay, with an invariant mass M{bW) between 160 and 190 GeV. The plot on the right 
in Fig jl3.59 shows the results after top subtraction. The subprocess ZW — > ZW provides 
the most relevant contribution. 

We intend to investigate the differences between the Standard Model and alternative 
scenarios in the high invariant mass region. In this analysis we have considered the Stan- 
dard Model without Higgs sector, hereafter denoted no-Higgs scenario, which acts as an 
upper limit for new physics effects in models like those mentioned in Ref.^4i. Despite 
longitudinally polarized gauge bosons become strongly interacting in absence of the Higgs 
particle, detecting any deviation still remains a challenge at LHC because the rise of the 
cross section is masked by the decrease of parton luminosities at high energies. An analysis 
of selection cuts capable to enhance the difference between the no-Higgs and light-Higgs 
case could provide some guidance in the search for signals of new physics. As shown in 



Fig jl3.60 simple requirements of centrality of the final-state bosons achieve this result at 
0(«Im)- 

The procedure we adopted for separating jets imposes a minimum invariant mass of 
60 GeV for each pair of quarks. An alternative approach consists in requiring a minimum 



AR separation among coloured particles and is discussed in Tab jl3.18| We find out that 
the no-Higgs case is more sensitive than light-Higgs to cuts based on AR. Indeed, in the 
model without Higgs the outgoing vector bosons are favoured to be more central and they 
have a larger pT- As a consequence of Lorentz boost they decay producing two quarks 
with small relative angle which are most likely to merge into one jet unless a sufficiently 
small AR is selected. The larger is the minimum separation required, the smaller the 
number of expected events. This is an indication that low AR cuts should be studied to 
evidentiate new physics. 

On the other hand, alternative procedures for jet reconstruction based on the K± 
algorithm have been recently proposed |14j. which may prove useful in connection with 
this kind of studies as they lead to encouraging results in identifying hadronic decays of 
heavy bosons via a cut on the sub-jet separation scale. 
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Figure 13.60: On the left: invariant mass distribution for M{V Z) > 800 GeV. The full 
and long-dashed line refer to the no-Higgs case, the short-dashed and dot-dashed ones to 
Mh = 200 GeV. All results satisfy the standard acceptance and selection cuts reported 
in Tab jl3.17 For the long-dashed and dot-dashed histograms we have further required 
\ri{i^i~)\ < 2 and \r](qq from y)| < 2. On the right: pseudorapidity distribution of the Z 
reconstructed from leptons. 



13.51 0{a^j^j^j)+0{a'^j^a1) results in the semileptonic /iz/^ chan- 
nel 

We will now investigate how much the sensitivity to the VV scattering signal is affected by 
the QCD irreducible background, which is expected to considerably dilute the differences 
between the light-Higgs and no-Higgs scenarios evidentiated in the previous section. To 
this purpose, we consider the full set of parton- level processes involved at ©(a^^^) + 

o(«iM«l), 



9Q 



qqqqUUf, 
gqqq^u^. 



99 
qq 



qqqqiiiy^, 
99qqf^^fi 



together with a selection procedure as close as possible to the actual experimental practice, 
without resorting to any flavour information other than the one which a typical 6-tagging 
algorithm is able to provide. The pz of the neutrino is approximately reconstructed by 
imposing the invariant mass of the two leptons equal to the W boson nominal mass: 



(p'' + p' 



u\2 



w 



This equation has two solutions. 



ap'z 



Pz 



■p. 



fM2 



where 



a 



+p'iPx+p'iPl 



(13.40) 



(13.41) 



(13.42) 



The event is rejected if the discriminant of Eq.(13.41) is negative, otherwise both the 
solutions are required to pass the selection cuts. 

All the results presented in this section have been obtained using the CTEQ5L PDF 
set and the QCD coupling constant running at the scale 



= Mlp + pT{topf 



(13.43) 
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Mcut 


no 


Higgs 


Mh 


= 200 GeV 


Ratio 


800 GeV 


31 


(14,17) 


12 


(7,5) 


2.59 


900 GeV 


25 


(12,13) 


8 


(5,3) 


3.12 


1.0 TeV 


19 


(9,10) 


6 


(4,2) 


3.16 


AR = 0.4 


Mcut 


no 


Higgs 


Mh 


= 200 GeV 


Ratio 


800 GeV 


18 


(8,10) 


10 


(6,4) 


1.80 


900 GeV 


12 


(5,7) 


6 


(4,2) 


2.00 


1.0 TeV 


8 


(4,4) 


4 


(2,2) 


2.00 


AR = 0.5 


Mcut 


no 


Higgs 


Mh 


= 200 GeV 


Ratio 


800 GeV 


12 


(5,7) 


8 


(5,3) 


1.50 


900 GeV 


8 


(4,4) 


5 


(4,2) 


1.60 


1.0 TeV 


5 


(2,3) 


3 


(2,1) 


1.60 



Table 13.18: Number of events as a function of the minimum invariant mass of the ZV — > 
IJi'^lJL~ jj pair for C = 100 fb~^, having apphed the cuts in Tab, 13.17 15^. All events satisfy 
\'q{l!.^i~)\ < 2 and {"qi^qq from V)\ < 2. In brackets we show the contribution of the {ZW, 
ZZ) final states. 



where pT{top) is the transverse momentum of the reconstructed top, for all processes in 
which a i or t can be produced. For all the other processes the scale has been evaluated 



as in Eq.( 13.39). 



Mh = 200 GeV 




^EW+QCD 


all events 


0.89 pb 


80.8 pb 


top events 


0.52 pb 


71.6 pb 


ratio top/all 


0.58 


0.89 



Table 13.19: Contribution of W/single t to the total cross section with standard acceptance 



cuts only (see the left part of Tab, 13.20). Comparison between results at 0{a%f^i) (EW) 
and 0{a%j^j) + 0[a%j^ja'l) (EW+QCD). Interferences between the two perturbative orders 
are neglected. 



It should be clear from the results shown in Tab 13.19 that suppressing the top back- 
ground is the primary objective to achieve. In this analysis we assume the possibility to 
tag 6-jets in the central region with 0.8 efficiency for |7y| < 1.5, which allows to discard part 
of the events containing h quarks in the final state. We impose additional cuts against top 
on the invariant mass of triplets of type {jjj} and {j'^i/}, where j denotes any final-state 
quark or gluon. In order to isolate two vector boson production, kinematical cuts are 
applied on the invariant mass of the two most central jets, which are associated in our 
analysis to aW oi Z decaying hadronically. The VV fusion signature is further isolated 
by requiring a minimum A-q separation between the two forward/backward jets. 

At this stage, however, any attempt to appreciate differences between Higgs and no- 
Higgs scenarios at high invariant masses would still be vain. This is essentially due to 



the fact that the contribution of the QCD diagrams depicted in Fig, 13.57) ^c,d) is not sub- 



stantially affected by the above-mentioned selection criteria. Investigating the differences 
between the kinematics of VV scattering and VV + 2 jets, we have identified additional 
cuts that serve our purpose. As the background dominates in the phase space regions 
characterized by one vector boson and the forward/backward jet produced with small 
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invariant mass, a viable method of taming VV + 2 jets consists in applying cuts on the 
Pt and rj of the W reconstructed from leptons as well as on the invariant mass of the W 
plus one of the two tag jets. All details about the selection cuts applied are reported in 
Tab JTMl 



Acceptance cuts 


Selection cuts 


PT{i^,j) > 10 GeV 


6-tagging for < 1.5 (80% efficiency) 


E{£^,j) > 20 GeV 


\M{jjj;je^iyrec)-Mtop\ > 15 GeV 


ivini < 3 


70GeV < M{jcjc) < 100 GeV 


\v{j)\ < 6.5 


M{jfjb) < 70 GeV ; M{jfjb) > 100 GeV 




M{jj) > 60 GeV 


^r]{jfjb) > 4 


Prii^i^rec) > 100 GeV 




M{jf/h£^urec) > 250 GeV 



Table 13.20: List of kinematical cuts applied in all results on the ^uf^ channel, j denotes 
any final-state quark or gluon, while £^ is the charged lepton. The subfixes c,f,b mean 
central, forward, backward respectively. Urec is the neutrino reconstructed following the 



prescription of Eq.( 13.40) 





no Higgs 


Mh = 200 GeV 




a 


£=100 fb"' 


a 


£=100 fb"' 


ratio 


all events 


12.46 fb 


1246 ± 35 


13.57 fb 


1357 ± 37 


0.918 


Mcut = 0.8 TeV 


3.19 fb 


319 ± 18 


1.45 fb 


145 ± 12 


2.200 


Mcut = 1.2 TeV 


1.28 fb 


128 ± 11 


0.41 fb 


41 ± 6 


3.122 


Mcut = 1.6 TeV 


0.60 fb 


60 ± 8 


0.14 fb 


14 ± 4 


4.286 



Table 13.21: Integrated 0{a%j^j) cross section for M{jcjcl'^) > Mcut and number of ex- 



pected events after one year at high luminosity having applied the cuts listed in Tab 13.20 



0(a|^) + 


no Higgs 


Mh = 200 GeV 




a 


£=100 fb"i 


a 


£=100 fb"' 


ratio 


all events 


40.70 fb 


4070 ± 64 


40.73 fb 


4073 ± 64 


0.999 


M,„i = 0.8 TeV 


7.61 fb 


761 ± 28 


5.14 fb 


514 ± 23 


1.481 


Mcut = 1.2 TeV 


2.53 fb 


253 ± 16 


1.73 fb 


173 ± 13 


1.462 


Mcut = 1.6 TeV 


1.00 fb 


100 ± 10 


0.55 fb 


55 ± 7 


1.818 



Table 13.22: Integrated ©(al^J + 0{aij,,al.) cross section for M{jcjcli^) > Mcut and 
number of expected events after one year at high luminosity having applied the cuts listed 



in Tab, 13.20 Interferences between the two perturbative orders are neglected. 



Fig 13.6T|l3.62| illustrate the final results of our analysis, showing that the top back- 
ground is basically under control. VV + 2 jets still provides a non-negligible contribution 
over the whole invariant mass spectrum, nevertheless differences between light-Higgs and 
no-Higgs can be appreciated. In Tab, 13.21|13.22 we show the integrated cross section at 
high energies as a function of the minimum invariant mass, comparing results of the pure 
EW and EW-I-QCD cases. Despite reducing the ratio between no-Higgs and light-Higgs 
cross sections, the inclusion of QCD background seems not to compromise the possibility 
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of finding signals of EWSB at LHC. We find that about 500 events are expected above 800 
GeV after one year of high luminosity running (£ = 100 fb~^) in case of a Higgs boson 
with mass 200 GeV. The Higgsless model predicts about 250 more events in accordance 
with the enhancement of the VV differential cross section at high energies. These num- 
bers refer to the muon channel only, and are obviously improved by summing up the muon 
and electron channels. It should nevertheless be noticed that imposing a minimum AR 
separation among coloured particles could degrade these preliminary results and requires 
further investigations. 

I Mh = 200 GeV I I NO Higgs | 
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No Higgs 



0.18 
0.16 
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I top OCD backgrojnd 




60 180 200 220 240 260 280 300 320 340 360 380 



xlO"* 



0.18- 
0.16 
0.14 
0.12 
0.1 
0.08 
0.06 
0.04 
0.02 
0, 



EW|no Hggs) 
I non-tcp QCD background 
I tcp OCDbarkgrcjnd 




160 180 200 220 240 260 280 300 320 340 360 380 



Figure 13.61: Invariant mass distribution of the two leptons and the two most central jets 
in the Standard Model with a light Higgs (on the left) and in the no-Higgs scenario (on 
the right). The cuts applied are hsted in Tab jl3.20[ 0{a%j,j) (EW) and (r'(a|^,,a|) (QCD) 
contributions to the differential cross section have been isolated and are shown separately. 
The QCD contributions are further split into top background (in blue) and VV + 2 jets 
(in green). 



13.52 Conclusions 

The results presented in this note should be intended as the starting point of a larger 
complete study which aims at providing an up-to-date and hopefully more reliable answer 
on the sensitivity of Boson-Boson Scattering as a probe of the symmetry breaking mech- 
anism at LHC. In spite of theoretical uncertaintes tipically estimated of C(10)%, these 
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Figure 13.62: Invariant-mass distribution of the two leptons plus the two most central 
jets in the pure EW (©(al^J) and EW+QCD {0{a%j^j) + 0{aj^j^al)) case. Comparison 
between light-Higgs {Mh = 200 GeV) and no-Higgs results. Interferences between the 



two different perturbative orders are neglected. The cuts applied are listed in Tab, 13.20 



results put a benchmark at parton level to the possibility of extracting the VV scattering 
signal from its irreducible background. There are clear indications that, with a complete 
calculation and appropriate kinematical cuts, the Higgs and no-Higgs cases show appre- 
ciable differences even after including the most relevant QCD contributions. We find in 
our analysis that the number of events expected in the Higgsless scenario is about 50% 
larger than the one predicted by the Standard Model with a 200 GeV Higgs boson. More 
refined selection criteria are under study in order to enhance further the Higgs signal as 
well as any possible deviation from the Standard Model expectations in the high invariant 
mass region. As an evolution for the near future, the analysis will be extended to cover 
V + 4:jets contributions so that a full estimate of the QCD background in the semileptonic 
channel will be available. 
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Searching for Extra Neutral 
Interactions at the LHC 



Roberta Armillis, Claudio Coriand, Alon E. Faraggi, Marco Guzzi, Nikos Irges 

Abstract 

We present a brief overview of some aspects of the theory and phenomenology of models 
containing extra neutral currents in abelian extensions of the Standard Model. We illus- 
trate the mechanism of anomaly cancellation as a way to infer the charge assignments of 
the fermion spectrum and then briefly discuss a variant of this approach where an effective 
action is rendered anomaly-free by a mechanism involving an axion via a Wess-Zumino 
term. Measurements at the LHC on the Z resonance in leptoproduction will be able to 
exclude a class of these models for variations of the cross section at the level of 4%, which 
is obtained at larger values of the anomalous coupling {g^ ~ 1). The anomalous nature is 
unlikely to be resolved with an inclusive NNLO analysis. 

14.53 Introduction 

The search for extra abelian gauge interactions will surely be an important component 
of the experimental program at the LHC. In fact, extra neutral interactions, or "extra 
Z'" , as they are commonly known, are predicted by several theoretical constructions such 
as Grand Unified Theories, Superstring Theory, and, more recently, by a class of models 
characterized by intersecting branes, just to mention a few (see [1] , [2] , |3] ) . One of the most 
relevant channels useful for the search of extra neutral currents is lepton pair production, 
at an invariant mass of the lepton pair not too large compared to the Z mass, since 
the invariant mass distribution da/dM"^ is rapidly falling, so to allow to gather enough 
statistics and separate the signal from the Standard Model background. Being the mass 
of the extra gauge boson (or bosons) undetermined, as is the coupling constant of the 
extra neutral current, a scan of the entire high energy tail of the Drell-Yan distribution is 
required. The width of these resonances is also undetermined by the theoretical models 
and may change considerably depending on the underlying assumptions of each theoretical 
construction. We will summarize below some of the main motivations of these searches 
discussing both anomaly-free models and a class of models where the anomaly is cancelled 
via the inclusion of a pseudoscalar, the Axi-Higgs [5j , which is part of the Higgs sector. 
In models of this type the phenomenology of the extra Z' shows some very distinctive 
features at the level of trilinear gauge interactions. A more detailed analysis of the topics 
addressed in this contribution can be found in [H El El |9l [TOl [TTl [12] . The numerical study 
that we present is performed on the abelian extension of [13], whose analysis is extended 
here, in part, to next-to-next-to-leading order in QCD. Other recent analysis addressing 
a general family-blind U{\) is in [T3] and briefly highlighted in the next section. A more 
detailed analysis of some of the issues addressed here are presented in [TU] . 
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14.54 Non Anomalous ?7(l)'s 



The interaction of an extra neutral gauge boson U{1) with the fermions of the Standard 
Model requires a suitable definition of their charges with respect to the new additional 
neutral interactions. Some of the most useful and powerful constraints in fixing the charges 
of these models is the requirement that they are free of anomalies. We will bring in a simple 
example to illustrate how the cancellation works for the case of one extra gauge factor, 
U{l)z, assuming a chiral spectrum that includes also a set of extra fermions z^^, with 
k = 1,2,..., n indicating a certain number of right-handed neutrinos. For simplicity we 
will assume that the interaction is family universal, which means that the cancellations 
take place generation by generation. We consider a model with the gauge symmetry 
5^7(3) X SU{2) X ^7(1) x C/(l)^. Breaking the gauge symmetry down to SU{3) x U{l)em 
requires an extra scalar sector compared to the SM. This can be achieved in different ways, 
for instance by including, beside the Higgs doublet, a SU{2)\y singlet (p, whose vacuum 
expectation value is v^p. In a more general framework, the constraints on the interactions 
of the U{l)z gauge boson with the fermions of the SM are relaxed if two Higgs doublets. 
Hi and H2, together with an extra scalar (p, are introduced (see |13]). 



14.55 Phenomenological Models: An Example 

Following [13] we consider three generations of quarks, g]^,-u^,(i^, and leptons, l]^,e^j^, 
i = 1,2,3, and a number n of right-handed neutrinos, u^, k = 1, ...,n, which are singlets 
under SU{3)c x SU{2)\y. We label the U{l)z charges as Zg, z^, Zd, zi, z^, for the standard 
model fermions, and are the charges of the right-handed neutrinos. The Higgs sector 
of the Standard Model is enlarged with an extra singlet (p, as shown in Table 1. is 
the electroweak Higgs. The cancellation of the anomalous interactions, for instance those 
involving the [5'C/(2)vv']^f^(l)z and \S\J{Z)c^\J{X)z gauge currents requires that 

zi = -3zg 

Zd = 2zq - Zu. (14.44) 
Similarly, the [C/(l)y]^C/(l)^ anomaly cancellation then implies that 

Ze = -2Zq - Zu. (14.45) 



Eqs. (14.44) and (14.45) together lead to the conclusion that only two independent real 
parameters, Zq and z^, describe the allowed U{l)z charges of the quarks and f7(l)y-charged 
leptons. Equivalently, the U{\)z charges may be expressed as a linear combination of Y 
and B-L: {zu - Zq)Y + {Azq - Zu){B - L). 

Additional restrictions on the U{\)z charges are imposed by the mixed gravitational- 
U{l)z and [^7(1)2]^ anomaly cancellation conditions 

1 " 

ly^Zk = -Azq + Zu , (14.46) 



Y,zk] =9^zl. (14.47) 



\fc=l / k=l 

For n = 1 or 2, the charge assignments compatible with these equations constrain 
[7(1)2 to be proportional to the hypercharge, giving a solution termed "sequential", while 
more general solutions can be found already for n = 3, which can be expressed in terms 
of a free parameter. The analysis of the possible neutrino mass terms compatible with a 
given charge assignment can also be used as a way to select the most interesting solutions 
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Table 14.23: Charge assignement for the Appelquist model 



of these equations. For instance, in the case n = 3 both Dirac and Majorana mass terms 
are possible in this minimal model (see [2] for more details). It is understood that a 
modification of the Higgs sector renders the study more involved [13] . 

14.56 Heterotic— string inspired Z' 

The heterotic-string gives rise to effective field theories that descend from the Eg x Eg 
or 5*0 (32) groups of the ten dimensional theories. The first case gives rise to additional 
Z's that arise in the 50(10) and Eq extensions of the Standard Model. A basis for the 
extra Z' arising in these models is formed by the two groups U{1)^ and U{1)^ via the 
decomposition Eq 50(10) x U{1)^ and 50(10) ^ 5C/(5) x U{1)^ [2^. Additional, flavor 
non-universal C/(l)'s, may arise in heterotic Eg x E^ string models from the U{1) currents 
in the Cartan subalgebra of the four dimensional gauge group, that are external to Eq. 
Non-universal Z's typically must be beyond the LHC reach, to avoid conflict with Flavor 
Changing Neutral Currents (FCNC) constraints. Recently |6] a novel Z' in quasi-realistic 
string models that do not descend from the Eg x Eg has been identified. Under the new 
C/(l) symmetry left-handed components and right-handed components in the 16 spinorial 
50(10) representation, of each Standard Model generation, have charge —1/2 and +1/2, 
respectively. As a consequence, the extra ?7(1) is family universal and anomaly free. 
It arises in left-right symmetric string models, in which the 50(10) symmetry is broken 
directly at the string level to 5C/(3) xU{1)b~lX SU{2) l x SU{2)rX [/(l)^' x [/(I)" x hidden 
|15j . The C/(l)" are flavor dependent U{l)s that are broken near the string scale. The 
Standard Model matter states are neutral under the hidden sector gauge group, which in 
these string models is typically a rank eight group. It is important to note that the fact 
that the spectrum is derived from a string vacuum that satisfies the modular invariance 
constraints, establishes that the model is free from gauge and gravitational anomalies. The 
pattern of U{l)z' charges in the quasi-realistic string models of ref. [15] does not arise in 
related string models in which the 50(10) symmetry is broken to a different subgroup. 

The important function of this Z' is that it forbids dimension four, five and six proton 
decay mediating operators [16]. The extra U{1) is anomaly free and family universal. 
It allows the fermions Yukawa couplings to the Higgs field and the generation of small 
neutrino masses via a seesaw mechanism. The existence of an extra Z' at low energies 
is motivated by proton longevity, and the suppression of the proton decay mediating 
operators. String models contain several U{1) symmetries that suppress the proton decay 
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mediating operators. However, these are typically non-family universal. They constrain 
the fermion mass terms and hence must be broken at a high scale. Thus, the existence 
of a U{1) symmetry that can remain unbroken down to low energies is highly nontrivial. 
The U{1) symmetry in ref. \15\ [6] satisfies all of these requirements. Furthermore, as 
the generation of small neutrino masses in the string models arises from the breaking 
of the B — L current, the extra U{1) allows lepton number violating terms, but forbids 
the baryon number violating terms. Hence, it predicts that i?-parity is violated and its 
phenomenological implications for SUSY collider searches differ substantially from models 
in which i?-parity is preserved. 

14.57 Anomalous U (l)'s: Cancelling the anomalies via higher 
dimensional operators 

The simplest effective theory that allows a cancellation of the gauge anomalies includes 
a Wess-Zumino term with a shifting axion. The role of the axion has been a matter of 
debate in the past [17^, since with a suitable gauge choice one can set the axion to vanish. 
Undoubtly, the axion can also be interpreted as the phase of an additional (second) Higgs 
when either large Yukawa couplings or a large vev of this additional Higgs is responsible 
for the decoupling of one or more chiral fermions [7]. According to this picture, an effec- 
tive anomalous theory is then the result of the partial decoupling of part of the fermion 
spectrum, leaving the left-over fermions in a reducible representation. The induced Wess- 
Zumino term is of the form b/MiF A F, a dimension-5 operator, where b is an axion and 
F is the field strength of a gauge interaction. The 1/Mi suppression of the axion-gauge 
interaction is directly related to the vev of the Higgs and b is its pseudoscalar phase [7] . 
The scale Mi is completely unrelated to the mechanism held responsible for the generation 
of a mass for the axion and remains a free parameter (see [18] ) . The mechanism of partial 
decoupling is a generic feature of effective anomalous theories that brings into the effective 
action a gauged axion. 

In general, models based on intersecting branes predict similar structures, though the 
axion of the effective theory is not a consequence of partial decoupling. The Standard 
Model, in this case, is enlarged with several extra anomalous C/(l)'s, one combination of 
which is anomaly free and is identified with the hypercharge, while the remaining U(l)'s 
are accompanied by axions. These models predict a single physical axion (the Axi-Higgs) 
[5] and contain generalized Chern-Simons terms in the effective action [5] , [19] , [7] . These 
terms are necessary in order to render the abelian interactions anomaly-free. On the 
contrary of non-anomalous U{1) models, which have similar properties and difficult to 
discern at the LHC, anomalous U(l)'s show some interesting features both in Drell-Yan 
at NNLO [12] , due to non-cancelling anomalous contribution and to the possibility of an 
axion exchange. We refer to [3 |^ for a discussion of the detection of modified trilinear 
interaction in the neutral sector ( Z77 vertex). 

14.58 The detection of extra Z-primes at the LHC 

We quantify below the rates for the invariant mass distributions in the case of the model 
of [13] . The approach follows closely |2U| EU [22] • Our convention for the couplings of the 
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Z = Z^Z' to the quarks and leptons are written below 
9v^'' = 2 



= 2 



9A 



z',j 



(14.48) 



where j is an index that runs over the quarks and the lepton, and sin Oy/ = Sw, cos 6w = Cw 
for brevity. The mixing in the neutral sector appears through the parameter e, defined as 



4 cos 2 9w 



K,+4.)[l + 0(6^)] 



ggz 



4cos^H."^'^^"^^+"^^'^^^ 



(14.49) 



where we have chosen [13] zh^ = = 0, = 246 GeV and tan /? = 1 and have defined 
g = e/s\D.9w,9Y = e / cos 9w- Precision studies at LEP constrain e to be smaller than 
10-3. 

The decay rates into leptons for the Z and the Z' are universal and are given by 



T{Z II) 



9 



1927rc2 



1 + + 1.409^^ - 12.77^^ 



(14.50) 



where i = u,d,c,s and Z = Z, Z' . For the Z' and Z decays into heavy quarks we have 
used (see Ref.[^ 



r{z ^ bb) 



V{Z tt) 



Ncae 



-Mzlig^f 



Co + Ci 



as{Mz) 



+ ■ 



+ (Qa ) 



-MzPiig^ 



^,t^2(3-/?') 



Ai + 



1 If A 



(14.51) 



where the coefficients cq, ci, are defined in Eq. (4.1) of [23] and where as in [23] we 

have defined 



1-4 



TT _3 + P fn _ 3 
2/3 4^ V2 ^ 4^ 



A2 = <i 1+ gQ.(A^2) 



TT /19 22/3 7l3^\ f n 3 
2/3 ^ VIO ^ ~5~ ^ ^ 



(14.52) 



2 47r^ 

Then, the total decay rate for the Z and Z' is obtained by summing over each fermionic 
contribution, for instance 

Tz' = Yl ^ ^*'^*) + ^ + ^ + 3r(^' ^ + 3r(Z' ^ zy,i7;). 

i=u,ci,c,s 

(14.53) 
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14.58.1 Calculation of the point-like cross section in the Z' case 

We come now briefly to discuss the quantiflcation of the invariant mass distributions 
around the resonances in leptoproduction. 

The colour-averaged inclusive differential cross section for the reaction li\ + H2 
h + h + X, is given by |25] 

Tav{Q^M^)Wv{T,Q^) r=^, (14.54) 



where V refers to the generic vector bosons Z, 7, Z' and where all the hadronic initial state 
information is contained in the hadronic structure function which is deflned as 

Wv/(t, Q^) = / dxi / dx2 I dx5{T - xxiX2)PDij{xi,X2, fJ''j?)Aij{x,Q'^, fi%) , 
Jo Jo Jo 

(14.55) 

where the quantity PDij{xi, X2, ^ip) contains all the information about the parton distri- 
bution functions and their evolution to the //^ scale and the A's are the flavour-dependent 
hard scatterings. 

The point-like cross sections for the case of the Z' boson are given by 
ctz'KQ ) ~ 



4Mz' sin^ 9w cos2 OwNc {Q^ - M|,)2 + M|,r|, 



0'Z','y{Q ) 



2\ _ '^^em 



QNc sin2 Ow cos2 Ow Q^{Q^ - M^Y + Ml,Tl, ' 



(yz',z[Q ) ^ 



Z',l z.l , Z',l z,i 
9v 9v +9a 9a 



96 sin^ Ow cos4 OwNc 

{Q^ - Ml){Q^ - Ml,) + MzTzMz'Tz' 
' [(Q2 _ MlY + MlTl] [(Q2 _ M|)2 + M|r|] 



(14.56) 



In Table 114.241 we have shown the numerical results for the total cross sections on the Z 
peak in the different models considered in our analysis. In the flrst line of each column 
we show the results for the total cross section in [fb], in the second line the total width 
Tz', expressed in GeV, and in the third line the observable atot x BR{Z' — > //), where 
BR{Z' — > //) = T zi^ii/^ z' ■ AH these quantities refer to the value of the coupling constant 
listed in the flrst column. 



In Table 14.25| we have shown a comparison between the Drell-Yan NNLO invariant 



mass distribution [21 for the Free fermionic model (anomaly free) in the TeV region 
and the SM background for different values of the coupling (see Ref. [TO]). We have 
performed the PDF evolution with CANDIA [22], and we have chosen the MRST [26] 
set as input distributions. The mass of the extra Z-prime has been taken Mz' = 2.5 
TeV, while tan (3 = 40. We observe that in correspondence of the peak value of the extra 
resonance, the cross section is enhanced of about 2 order of magnitude with respect to 
the SM, while moving away from the peak, the value of the Free fermionic cross section 
decreases rapidly and the difference with respect to the background is around 1-2%. 

14.58.2 Results for anomaly-free and anomalous models 

An analysis of the anomalous effects in Drell-Yan and in double prompt photon can be 
found in where several comparisons between anomaly- free and anomalous models 
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(the MLSOM) are reported. In Fig. 14.63 we show a comparison between the MLSOM 
and the anomaly-free extensions. We have included the ^r/ scale dependence, which 
appears as a band, and the variations with respect to Qz- As shown in this figure, the 
red lines correspond to the MLSOM, the green ones to the free fermionic model, the blue 
ones to the U{\)b-l model and the purple ones to the U{\)qj^u model. The first peak 
corresponds to Qz = 0.1, the second to Qz = Qy, the third to Qz = §2 and the fourth to 
gz = I- The width of each peak gets larger as gz grows, but the peak- value of the cross 
section decreases. Different choices of gz correspond to slightly different values of the mass 
of the extra Z' because of the relation between the Stiickelberg mass Mi and Mz' given 
in [7j. For a fixed value of the coupling, the effects due to the variations of the scales fiR 
and fip become visible only for gz = 1 and in this case they are around 2-3%. In the case 
gz = I (red line), the uppermost lines correspond to the choice = 2(5, = l/2fiF 
and fiR = 2^F, while the lowermost lines correspond to the choice fip = Q, fj,R = 1/2 fip 
and fiR = 2iJip. The peak- value obtained for the anomalous model is the largest one, with 
a cross section which is around 0.022 [pb/GeV], while the free fermionic model appears 
to be the smallest with a value around 0.006 [pb/GeV]. A sizeable coupling of the extra 



Z prime affects also the width and the height of the Z resonance. We show in Fig. 14.64 



that there is an overlap (red band) between the theoretical uncertainty for the invariant 
mass distribution in Drell-Yan on the Z resonance, due to the change of the perturbative 
order in the SM (green band) and in the MLSOM (blue band). From this figure it is 
evident that for larger values of the anomalous coupling {gz > 1) a class of anomalous 
models based on intersecting branes can be excluded, due to the small overlap for gz ^ 1 
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Table 14.24: Total cross sections, widths and atot x BR{Z' II), where BR{Z' II) = 
^ Z'^ii/^ z' , for the MLSOM and three anomaly-free extensions of the SM. 

We have presented a brief discussion of the search of extra neutral interactions at 
the LHC. While the anomaly- free construction are generated quite automatically in the 
context of GUT's and low energy string models, the possibility of detecting an anomalous 
gauge interaction can not be excluded. In this second case the anomaly cancellation 
procedure of the theory requires the introduction of a Wess-Zumino counterterm and the 
effective action contains a physical axion as a fingerprint of partial decoupling of part of 
the fermion spectrum. 

We have presented some numerical predictions for the NNLO invariant mass distribu- 
tions in leptoproduction using different anomaly free models as a case study, and shown 



153 




Figure 14.63: Comparisons among anomalous Drell-Yan in the MLSOM versus several 
anomaly- free models. 



that the cross sections are rather sensitive to the masses of the extra neutral gauge bosons 
and to the coupling g^. In this first analysis, the masses of the extra gauge bosons and 
the extra coupling can be considered in general, free parameters of the theory. In the 
anomalous case, anomaly effects in Drell-Yan have been found to be small. However, 
while the nature (whether anomalous or not) of the extra Z' will be difficult to resolve 
at the LHC in Drell-Yan or in double photon, precision study on the Z resonance can 
be used to set exclusion limits on many of these models especially at larger values of the 
anomalous coupling (g^ ~ 1). 



da""'^°/dQ [pb/GeV] for the FF model with Mz' = 2.5 TeV, tan/3 = 40, Candia evol. 
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Table 14.25: NNLO cross sections for the Free Fermionic model with a Mz' = 2.5 TeV for 
values of the coupling constant g^ larger than g^ = 0.1 
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Figure 14.64: Zoom on the Z resonance for anomalous Drell-Yan in tlie ~ fiji = Q at 
NLO/NNLO for all the models. 
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